Eoft now is output as function of euclid time

src/io.jl

@@ -101,6 +101,7 @@ function write_log()
     println(log_file,"maxiter      = ", params["Measurements"]["maxiter"])
     println(log_file,"Flow times   = ", flow_times)
     println(log_file,"N noise      = ", N_noise)
+    println(log_file,"t_source     = ", params["Fermion"]["tsource"])
 
     flush(log_file)
     return nothing

src/meas.jl

@@ -21,7 +21,8 @@ function load_fields()
     global pp_corr_t = fill(zero(Float64),(lp.iL[4],N_noise,length(flow_times)));
     global ap_corr_t = fill(zero(ComplexF64),(lp.iL[4],N_noise,length(flow_times)));
 
-    global Eoft = Array{Complex{Float64}}(undef,1+length(flow_times));
+    global Eoft = Array{Complex{Float64}}(undef,1+length(flow_times),lp.iL[4]);
+    global Eofpla = Array{Complex{Float64}}(undef,lp.iL[4],lp.npls);
 
     global phat_t = Array{Float64}(undef,N_noise,length(flow_times));
     global Quark_cond = Array{Complex{Float64}}(undef,N_noise,length(flow_times));
@@ -40,7 +41,8 @@ function two_pt()
     println(log_file,"Measuring 2pt...")
     flush(log_file)
 
-    Eoft[1] = Eoft_plaq(U, gp, lp, ymws)
+    Eoft_plaq(Eofpla,U, gp, lp, ymws)
+    Eoft[1,:] .= sum(Eofpla,dims = 2)
 
     for noi in 1:N_noise
 
@@ -68,7 +70,10 @@ function two_pt()
             _,epslist = flw_adapt(U, psi, int, delta_t, neps, gp, dpar, lp, ymws, dws)
             neps = epslist[end]
 
-            noi == 1 ? Eoft[fl+1] = Eoft_plaq(U, gp, lp, ymws) : nothing
+            if noi ==1
+                Eoft_plaq(Eofpla,U, gp, lp, ymws)
+                Eoft[fl+1,:] .= sum(Eofpla,dims = 2)
+            end
 
             pp_density .= Array(norm2.(psi))
             ap_density .= Array(dot.(psi,dmul.(Gamma{4},psi)))

src/utils.jl

@@ -41,7 +41,7 @@ function read_ff(name::String)
     Sigma = Array{ComplexF64}(undef,1,N_noise,nflow);
     Phat_t = Array{Float64}(undef,1,N_noise,nflow);
     Sigma_cfl = Array{ComplexF64}(undef,1,N_noise,nflow);
-    Eoft = Array{Complex{Float64}}(undef,1,1+nflow);
+    Eoft = Array{Complex{Float64}}(undef,1,1+nflow,iL[4]);
 
     pp_corr_mc = Array{Float64}(undef,0,iL[4],N_noise);
     ap_corr_mc = Array{ComplexF64}(undef,0,iL[4],N_noise);
@@ -50,7 +50,7 @@ function read_ff(name::String)
     Sigma_mc = Array{ComplexF64}(undef,0,N_noise,nflow);
     Phat_t_mc = Array{Float64}(undef,0,N_noise,nflow);
     Sigma_cfl_mc = Array{ComplexF64}(undef,0,N_noise,nflow);
-    Eoft_mc = Array{Complex{Float64}}(undef,0,1+nflow);
+    Eoft_mc = Array{Complex{Float64}}(undef,0,1+nflow,iL[4]);
 
     while BDIO_get_uinfo(file) != 8
         BDIO_seek!(file)
@@ -116,7 +116,7 @@ function uwff(file::String)
     Ns = size(pp_corr_mc)[3]
     Nfl = size(pp_corr_t_mc)[4]
 
-    global Eoft = [uwreal(real.(Eoft_mc[:,i]),runame) for i in 1:Nfl+1]
+    global Eoft = [uwreal(real.(Eoft_mc[:,i,j]),runame) for j in 1:T, i in 1:Nfl+1]
 
     global pp_corr = [uwreal(sum(pp_corr_mc,dims = 3)[:,i,1]./Ns,runame) for i in 1:T]
     global ap_corr = [uwreal(real.(sum(ap_corr_mc,dims = 3)[:,i,1]./Ns),runame) for i in 1:T]