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

Small bug corrections

Browse source
This commit is contained in:
Alberto Ramos 2021-10-27 15:24:35 +02:00
parent 563eb55e48
commit 0b8869dc89
4 changed files with 5 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

2
src/LatticeGPU.jl
View file

@ -33,7 +33,7 @@ export vector_field, scalar_field, nscalar_field, scalar_field_point
include("MD/MD.jl") include("MD/MD.jl")
using .MD using .MD
export IntrScheme export IntrScheme
export omf4, leapfrog export omf4, leapfrog, omf2
include("YM/YM.jl") include("YM/YM.jl")

2
src/MD/MD.jl
View file

@ -33,7 +33,7 @@ end
omf2(::Type{T}, eps, ns) where T = IntrScheme{3,T}((r1omf2,r2omf2,r3omf2), eps, ns) omf2(::Type{T}, eps, ns) where T = IntrScheme{3,T}((r1omf2,r2omf2,r3omf2), eps, ns)
omf4(::Type{T}, eps, ns) where T = IntrScheme{6,T}((r1omf4,r2omf4,r3omf4,r4omf4,r5omf4,r6omf4), eps, ns) omf4(::Type{T}, eps, ns) where T = IntrScheme{6,T}((r1omf4,r2omf4,r3omf4,r4omf4,r5omf4,r6omf4), eps, ns)
leapfrog(::Type{T}, eps, ns) where T = IntrScheme{2,T}((0.5,1.0,0.5), eps, ns) leapfrog(::Type{T}, eps, ns) where T = IntrScheme{2,T}((0.5,1.0), eps, ns)
import Base.show import Base.show

2
src/YM/YM.jl
View file

@ -30,7 +30,7 @@ struct GaugeParm{T,G}
function GaugeParm{T}(::Type{G}, bt, c0, cG) where {T,G} function GaugeParm{T}(::Type{G}, bt, c0, cG) where {T,G}
function degree(::SU2{T}) where T function degree(::Type{SU2{T}}) where T
return 2 return 2
end end
function degree(::Type{SU3{T}}) where T <: AbstractFloat function degree(::Type{SU3{T}}) where T <: AbstractFloat

4
src/YM/YMhmc.jl
View file

@ -56,11 +56,10 @@ function hamiltonian(mom, U, lp, gp, ymws)
return K+V return K+V
end end
function HMC!(U, eps, ns, lp::SpaceParm, gp::GaugeParm, ymws::YMworkspace{T}; noacc=false) where T function HMC!(U, int::IntrScheme, lp::SpaceParm, gp::GaugeParm, ymws::YMworkspace{T}; noacc=false) where T
@timeit "HMC trayectory" begin @timeit "HMC trayectory" begin
int = omf4(T, eps, ns)
ymws.U1 .= U ymws.U1 .= U
randomize!(ymws.mom, lp, ymws) randomize!(ymws.mom, lp, ymws)
@ -87,6 +86,7 @@ function HMC!(U, eps, ns, lp::SpaceParm, gp::GaugeParm, ymws::YMworkspace{T}; no
end end
return dh, acc return dh, acc
end end
HMC!(U, eps, ns, lp::SpaceParm, gp::GaugeParm, ymws::YMworkspace{T}; noacc=false) where T = HMC!(U, omf4(T, int.eps, int.ns), lp, gp, ymws; noacc=noacc)
function MD!(mom, U, int::IntrScheme{NI, T}, lp::SpaceParm, gp::GaugeParm{T}, ymws::YMworkspace{T}) where {NI, T <: AbstractFloat} function MD!(mom, U, int::IntrScheme{NI, T}, lp::SpaceParm, gp::GaugeParm{T}, ymws::YMworkspace{T}) where {NI, T <: AbstractFloat}