lattice/latticegpu.jl
1
0
Fork 0
mirror of https://igit.ific.uv.es/alramos/latticegpu.jl.git synced 2025-11-04 16:36:03 +01:00
Code Issues Projects Releases Packages Wiki Activity Actions

Added basic flow routines

Browse source
This commit is contained in:
Alberto Ramos 2021-09-25 09:36:08 +02:00
commit b347d33dbf
1 changed files with 44 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

44
src/YM/YMflow.jl Normal file
View file

@ -0,0 +1,44 @@
###
### "THE BEER-WARE LICENSE":
### Alberto Ramos wrote this file. As long as you retain this
### notice you can do whatever you want with this stuff. If we meet some
### day, and you think this stuff is worth it, you can buy me a beer in
### return. <alberto.ramos@cern.ch>
###
### file: YMflow.jl
### created: Sat Sep 25 08:37:14 2021
###
function wfl_euler(U, ns, eps, lp::SpaceParm, ymws::YMworkspace)
for i in 1:ns
force_wilson(ymws, U, lp)
U .= expm(U, ymws.frc1, 2*eps)
end
return nothing
end
function wfl_rk3(U, ns, eps, lp::SpaceParm, ymws::YMworkspace)
for i in 1:ns
c0 = eps/2
force_wilson(ymws, U, lp)
ymws.mom .= ymws.frc1
U .= expm(U, ymws.mom, c0)
c0 = -34*eps/36
c1 = 16*eps/9
force_wilson(ymws, U, lp)
ymws.mom .= c0.*ymws.mom .+ c1.*ymws.frc1
U .= expm(U, ymws.mom)
c1 = 6*eps/4
force_wilson(ymws, U, lp)
ymws.mom .= c1.*ymws.frc1 .- ymws.mom
U .= expm(U, ymws.mom)
end
return nothing
end