Improvements in main file

main/sfwfl.jl

@@ -56,43 +56,50 @@ println(gp)
 
 wflw = wfl_rk3(Float64, 0.01, 1.0E-7)
 zflw = zfl_rk3(Float64, 0.01, 1.0E-7)
-println(zflw)
 println(" ## ")
 
 ymws = YMworkspace(SU3, Float64, lp)
 println(" # Plaquette: ", plaquette(U, lp, gp, ymws))
 
-flw_ns = round((0.3*l)^2/(8*wflw.eps)/10)
+flw_nc = 20
 
 Ecl = Array{Float64, 2}(undef, lp.iL[end], 6)
 Epl = Array{Float64, 2}(undef, lp.iL[end], 6)
 Qt  = Vector{Float64}(undef, lp.iL[end])
 
+ndtr = 10
 println(wflw)
 ymws.U1 .= U
-for j in 1:flw_ns
+tf = 0.0
+for j in 1:flw_nc+1
     Eoft_plaq(Epl, U, gp, lp, ymws)
     Eoft_clover(Ecl, U, gp, lp, ymws)
     Qtop(Qt, U, gp, lp, ymws)
     
-    t = j*wflw.eps
-    @printf("     FLOW   t=%8.4f:  %20.12e  %20.12e      %10.4e \n",
-            t, t^2*sum(Epl), t^2*sum(Ecl), sum(Qt))
-    flw(U, wflw, 10, gp, lp, ymws)
+    @printf("     FLOW   c=%8.4f:  %20.12e  %20.12e      %10.4e \n",
+            sqrt(8*tf)/l, l^3*sum(Epl), l^3*sum(Ecl), sum(Qt))
+
+    c = 0.6*j/flw_nc
+    dt = (cl)^2/8 - tf
+    flw_adapt(U, wflw, dt, gp, lp, ymws)
+    tf += dt
 end
 U = ymws.U1
 
 println(zflw)
-ymws.U1 .= U
-for j in 1:flw_ns
+tf = 0.0
+for j in 1:flw_nc+1
     Eoft_plaq(Epl, U, gp, lp, ymws)
     Eoft_clover(Ecl, U, gp, lp, ymws)
     Qtop(Qt, U, gp, lp, ymws)
     
-    t = j*zflw.eps
-    @printf("     FLOW   t=%8.4f:  %20.12e  %20.12e      %10.4e \n",
-            t, t^2*sum(Epl), t^2*sum(Ecl), sum(Qt))
-    flw(U, zflw, 10, gp, lp, ymws)
+    @printf("     FLOW   c=%8.4f:  %20.12e  %20.12e      %10.4e \n",
+            sqrt(8*tf)/l, l^3*sum(Epl), l^3*sum(Ecl), sum(Qt))
+
+    c = 0.6*j/flw_nc
+    dt = (cl)^2/8 - tf
+    flw_adapt(U, zflw, dt, gp, lp, ymws)
+    tf += dt
 end
 U = ymws.U1
 println(" # Plaquette: ", plaquette(U, lp, gp, ymws))