Typos varios

parent 2ba57494
...@@ -155,16 +155,16 @@ Returns the boundary to boundary propagator of the Schrodinger functional given ...@@ -155,16 +155,16 @@ Returns the boundary to boundary propagator of the Schrodinger functional given
""" """
function bndtobnd(bndpro::AbstractArray, U::AbstractArray, dpar::DiracParam, dws::DiracWorkspace, lp::SpaceParm{4,6,1,D}) where {D} function bndtobnd(bndpro::AbstractArray, U::AbstractArray, dpar::DiracParam, dws::DiracWorkspace, lp::SpaceParm{4,6,1,D}) where {D}
function krnl_bndtobnd!(psi::AbstractArray, bndp::AbstractArray, U::AbstractArray, lp) function krnl_bndtobnd!(psi::AbstractArray, bndp::AbstractArray, U::AbstractArray, lp::SpaceParm)
b=Int64(CUDA.threadIdx().x) b=Int64(CUDA.threadIdx().x)
r=Int64(CUDA.blockIdx().x) r=Int64(CUDA.blockIdx().x)
it = point_time((b, r), lp) it = point_time((b, r), lp)
if it == lp.il[end] if it == lp.iL[end]
psi[b,r] = gdagpmul(Pgamma{4,1},U[b,4,r],bndpro[b,r]) psi[b,r] = gdagpmul(Pgamma{4,1},U[b,4,r],bndpro[b,r])
else else
psi[b,r] = 0.0 * bndpro[b,r] psi[b,r] = 0.0*bndpro[b,r]
end end
return nothing return nothing
...@@ -174,7 +174,7 @@ function bndtobnd(bndpro::AbstractArray, U::AbstractArray, dpar::DiracParam, dws ...@@ -174,7 +174,7 @@ function bndtobnd(bndpro::AbstractArray, U::AbstractArray, dpar::DiracParam, dws
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_bndtobnd!(dws.sp, bndpro ,U, lp) CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_bndtobnd!(dws.sp, bndpro ,U, lp)
end end
res = -dpar.ct*CUDA.mapreduce(x -> x, +, dws.sp)./(lp.iL[1]*lp.iL[2]*lp.iL[3]) res = -dpar.ct*CUDA.mapreduce(x -> x, +, dws.sp) / (lp.iL[1]*lp.iL[2]*lp.iL[3])
return res return res
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