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

Bug in i/o

Browse source
This commit is contained in:
Alberto Ramos 2021-11-14 09:46:01 +01:00
parent 95f60c18b1
commit 7a438c69af
1 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/YM/YMio.jl
View file

@ -20,8 +20,8 @@ function import_lex64(fname, lp::SpaceParm)
dtr = [2,3,4,1] dtr = [2,3,4,1]
assign(i, V, i4) = SU3{Float64}(V[1,dtr[id],1],V[2,dtr[id],1],V[3,dtr[id],1], assign(id, V, i4) = SU3{Float64}(V[1,dtr[id],i3],V[2,dtr[id],i3],V[3,dtr[id],i3],
V[4,dtr[id],1],V[5,dtr[id],1],V[6,dtr[id],1]) V[4,dtr[id],i3],V[5,dtr[id],i3],V[6,dtr[id],i3])
Ucpu = Array{SU3{Float64}, 3}(undef, lp.bsz, lp.ndim, lp.rsz) Ucpu = Array{SU3{Float64}, 3}(undef, lp.bsz, lp.ndim, lp.rsz)
V = Array{ComplexF64, 3}(undef, 9, lp.ndim, lp.iL[3]) V = Array{ComplexF64, 3}(undef, 9, lp.ndim, lp.iL[3])
@ -32,7 +32,7 @@ function import_lex64(fname, lp::SpaceParm)
for i3 in 1:lp.iL[3] for i3 in 1:lp.iL[3]
b, r = point_index(CartesianIndex(i1,i2,i3,i4), lp) b, r = point_index(CartesianIndex(i1,i2,i3,i4), lp)
for id in 1:lp.ndim for id in 1:lp.ndim
Ucpu[b,id,r] = assign(i, V, i4) Ucpu[b,id,r] = assign(id, V, i3)
end end
end end
end end