Commit 6758aeee authored by Alberto Ramos's avatar Alberto Ramos

Added acc/rej (was missing!)

parent 041abadc
......@@ -148,8 +148,9 @@ if BCS == BC_SF_AFWB
end
println(flog, " # Thermalization")
HMC!(U,int,lp, gp, ymws, noacc=true)
for i in 1:nthm
dh[i], acc[i] = HMC!(U,int,lp, gp, ymws, noacc=true)
dh[i], acc[i] = HMC!(U,int,lp, gp, ymws)
pl[i] = plaquette(U,lp, gp, ymws)
if BCS == BC_PERIODIC
......@@ -175,7 +176,7 @@ println(flog, "\n # Measurements")
kfl = 0
for i in 1:nt
k = i + nthm
dh[k], acc[k] = HMC!(U,int,lp, gp, ymws, noacc=true)
dh[k], acc[k] = HMC!(U,int,lp, gp, ymws)
pl[k] = plaquette(U,lp, gp, ymws)
if BCS == BC_PERIODIC
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment