lattice/latticegpu.jl
1
0
Fork 0
mirror of https://igit.ific.uv.es/alramos/latticegpu.jl.git synced 2025-05-14 19:23:42 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

Minor changes to test in src/main/times.jl

Browse source
This commit is contained in:
Alberto Ramos 2023-01-24 18:41:39 +01:00
parent cd8a9d4f7c
commit 4944740ad2
1 changed files with 7 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

14
src/main/times.jl
View file

@ -2,13 +2,13 @@ using CUDA, Logging, Random, TimerOutputs
CUDA.allowscalar(false) CUDA.allowscalar(false)
import Pkg import Pkg
Pkg.activate("/lhome/ific/a/alramos/s.images/julia/workspace/LatticeGPU") Pkg.activate("../..")
#Pkg.activate("/home/alberto/code/julia/LatticeGPU") #Pkg.activate("/home/alberto/code/julia/LatticeGPU")
using LatticeGPU using LatticeGPU
# Set lattice/block size # Set lattice/block size
ntwist = (0,0,0,0,0,0) ntwist = (0,0,0,0,0,0)
lp = SpaceParm{4}((8,8,8,8), (4,4,4,4), BC_PERIODIC, ntwist) lp = SpaceParm{4}((32,32,32,32), (4,4,4,4), BC_PERIODIC, ntwist)
println("Space Parameters: ", lp) println("Space Parameters: ", lp)
# Seed RNG # Seed RNG
@ -23,7 +23,7 @@ PREC = Float64
println("Precision: ", PREC) println("Precision: ", PREC)
# MD integrator # MD integrator
int = omf4(PREC, 0.1, 10) int = omf4(PREC, 1.0/40.0, 40)
println(int) println(int)
println("Allocating YM workspace") println("Allocating YM workspace")
@ -41,10 +41,10 @@ println("Time to take the configuration to memory: ")
# Set gauge parameters # Set gauge parameters
# FIRST SET: Wilson action/flow # FIRST SET: Wilson action/flow
println("\n## WILSON ACTION/FLOW TIMES") println("\n## WILSON ACTION/FLOW TIMES")
gp = GaugeParm{PREC}(GRP{PREC}, 6.0, 1.0) gp = GaugeParm{PREC}(GRP{PREC}, 6.0, 1.66666666)
println("Gauge Parameters: ", gp) println("Gauge Parameters: ", gp)
flwint = wfl_rk3(PREC, 0.005, 1.0E-6) flwint = zfl_rk3(PREC, 0.005, 1.0E-6)
println(flwint) println(flwint)
println("Initial Action: ") println("Initial Action: ")
@ -54,8 +54,8 @@ println("Initial Action: ")
HMC!(U,int.eps,1,lp, gp, ymws, noacc=true) HMC!(U,int.eps,1,lp, gp, ymws, noacc=true)
pl = Vector{Float64}() pl = Vector{Float64}()
for i in 1:4 for i in 1:100
@time dh, acc = HMC!(U,int,lp, gp, ymws, noacc=true) @time dh, acc = HMC!(U,int,lp, gp, ymws)
println("# HMC: ", acc, " ", dh) println("# HMC: ", acc, " ", dh)
push!(pl, plaquette(U,lp, gp, ymws)) push!(pl, plaquette(U,lp, gp, ymws))
println("# Plaquette: ", pl[end], "\n") println("# Plaquette: ", pl[end], "\n")