get neig/own only for internal faces

exasim-project · Feb 15, 2025 · 3e5215b · 3e5215b
1 parent 139ea2c
commit 3e5215b
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/finiteVolume/cellCentred/interpolation/linear.cpp b/src/finiteVolume/cellCentred/interpolation/linear.cpp
@@ -41,10 +41,10 @@ void computeLinearInterpolation(
         exec,
         {0, dstS.size()},
         KOKKOS_LAMBDA(const size_t facei) {
-            size_t own = static_cast<size_t>(ownerS[facei]);
-            size_t nei = static_cast<size_t>(neighS[facei]);
             if (facei < nInternalFaces)
             {
+                size_t own = static_cast<size_t>(ownerS[facei]);
+                size_t nei = static_cast<size_t>(neighS[facei]);
                 dstS[facei] = weightS[facei] * srcS[own] + (1 - weightS[facei]) * srcS[nei];
             }
             else