wflow.lua

require 'common'
require 'gaugeact'

--trace(doTrace)

local nx = nx or 4
local nt = nt or 8
local seed = seed or os.time()
latsize = { nx, nx, nx, nt }
local nd = #latsize
qopqdp.lattice(latsize)
qopqdp.profile(0)
qopqdp.verbosity(0)
qopqdp.seed(seed)

vol = 1
printf("latsize =")
for k,v in ipairs(latsize) do vol=vol*v; printf(" %i",v); end
printf("\n")
printf("seed = %i\n", seed)

U = qopqdp.gauge()
if inlat then
   U:load(inlat)
else
   U:unit()
end

coeffs = { plaq=1 }
--coeffs = { plaq=1, rect=-0.05 }
eps = 0.025
tmax = nx*nx/16
tmax = 1
nsteps = tmax/eps

function stats(u, i)
  local ss,st = plaq(u)
  local t = eps*i
  local e = 6*(2-ss-st)
  printf("%g\t%g\n", t, e*t*t)
end

t0 = clock()

stats(U, 0)
for i=1,nsteps do
  wflow(U, coeffs, eps, 1)
  stats(U, i)
  if(i%10==0) then collectgarbage() end
end

t1 = clock()

printf("# time: %g seconds\n", t1-t0)