Bug fix

CliMA · Jul 19, 2023 · 849f18b · 849f18b
1 parent b15fe7d
commit 849f18b
Showing 1 changed file with 5 additions and 5 deletions.
diff --git a/src/Operators/finitedifference.jl b/src/Operators/finitedifference.jl
@@ -3478,14 +3478,14 @@ function copyto_stencil_kernel!(out, bc, space, bds)
     #TODO pass in (Nq, Nh)!
     Nq = 1
     Nh = 1
-    space = axes(out)
-    if space isa Spaces.ExtrudedFiniteDifferenceSpace
-        QS = Spaces.quadrature_style(space)
+    spaceout = axes(out)
+    if spaceout isa Spaces.ExtrudedFiniteDifferenceSpace
+        QS = Spaces.quadrature_style(spaceout)
         Nq = Quadratures.degrees_of_freedom(QS)
-        Nh = Topologies.nlocalelems(Spaces.topology(space))
+        Nh = Topologies.nlocalelems(Spaces.topology(spaceout))
     end
     nitems = Nq * Nq * Nh
-    if gid < nitems
+    if gid ≤ nitems
         h = cld(gid, Nq * Nq)
         j = cld(gid - (h - 1) * Nq * Nq, Nq)
         i = gid - (h - 1) * Nq * Nq - (j - 1) * Nq