ncluded PhatP correlator.

main.jl

@@ -29,7 +29,10 @@ while RUN_ON
 
     @timeit "HMC" Gauge_update()
 
-    @timeit "Two point function" two_pt()
+    @timeit "Two point function" begin
+        two_pt()
+        two_pt_rnd()
+    end
 
     @timeit "One point function" one_pt()
 

src/meas.jl

@@ -49,6 +49,7 @@ function two_pt()
         niter = propagator!(psi,U, dpar, dws, lp,params["Measurements"]["maxiter"], params["Measurements"]["tolerance"],params["Fermion"]["tsource"])
 
         println(log_file,"CG converged after ",niter," iterations with residue ",CUDA.mapreduce(x -> norm2(x), +, dws.sr))
+        flush(log_file)
 
         pp_density .= Array(norm2.(psi))
         ap_density .= Array(dot.(psi,dmul.(Gamma{4},psi)))
@@ -134,6 +135,7 @@ function one_pt()
             end
 
             println(log_file,"CG converged after ",niter," iterations with residue ",CUDA.mapreduce(x -> norm2(x), +, dws.sr))
+            flush(log_file)
 
             @timeit "CPU to GPU" copyto!(dws.sp,psi_CPU)
             Quark_cond[noi,fl] = sum(dot.(dws.sp,psi))