Fix in volume sum. Should be improved for GPU

parent 6b8b4a32
...@@ -156,8 +156,8 @@ function one_pt() ...@@ -156,8 +156,8 @@ function one_pt()
for t in 1:lp.iL[4] for t in 1:lp.iL[4]
for i in 1:lp.iL[1] for j in 1:lp.iL[2] for k in 1:lp.iL[3] for i in 1:lp.iL[1] for j in 1:lp.iL[2] for k in 1:lp.iL[3]
b,r = point_index(CartesianIndex{lp.ndim}((i,j,k,t)),lp) b,r = point_index(CartesianIndex{lp.ndim}((i,j,k,t)),lp)
Quark_cond[t,noi,fl] += dot(dws.sp[b,r],psi[b,r]) CUDA.@allowscalar Quark_cond[t,noi,fl] += dot(dws.sp[b,r],psi[b,r])
Quark_cond_cfl[t,noi,fl] += -sum(norm2.(dws.sp)) CUDA.@allowscalar Quark_cond_cfl[t,noi,fl] += -norm2.(dws.sp[b,r])
end end end end end end
end end
####### #######
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment