Tolerance as input in backflow

input/FerFlow.in

@@ -38,3 +38,4 @@ N_noise = 6
 Flow_times = [2.0,3.0]
 Nsaves = 25
 tsource = 4
+tol = 1.0e-6

src/io.jl

@@ -60,7 +60,7 @@ function load_structs()
     global dpar = DiracParam{Float64}(SU3fund,(1/(2*params["Fermion"]["kappa"])) - 4,params["Fermion"]["csw"],ntuple(i -> exp(((i!=4)*im*params["Fermion"]["theta"]/lp.iL[i]) + ((i==4)*im*params["Fermion"]["theta_t"]/lp.iL[i])),4),0.0,params["Fermion"]["ct"]);
     global dws = DiracWorkspace(SU3fund,Float64,lp);
     global ymws = YMworkspace(SU3,Float64,lp);
-    global int = wfl_rk3(Float64, params["Frontflow"]["epsilon"], 1.0)
+    global int = wfl_rk3(Float64, params["Frontflow"]["epsilon"], params["Backflow"]["tol"])
     global intsch = omf4(Float64,params["HMC"]["eps"], params["HMC"]["ns"]);
 
     global flow_times = params["Backflow"]["Flow_times"]