Fix: include global offset in monotonic spline

boutproject · Jan 31, 2025 · 2d72bab · 2d72bab
1 parent cd383e2
commit 2d72bab
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 2 deletions.
diff --git a/src/mesh/interpolation/hermite_spline_xz.cxx b/src/mesh/interpolation/hermite_spline_xz.cxx
@@ -170,6 +170,8 @@ void XZHermiteSplineBase<monotonic>::calcWeights(const Field3D& delta_x,
 #ifdef HS_USE_PETSC
   IndConverter conv{localmesh};
 #endif
+  [[maybe_unused]] const int y_global_offset =
+      localmesh->getYProcIndex() * (localmesh->yend - localmesh->ystart + 1);
   BOUT_FOR(i, getRegion(region)) {
     const int x = i.x();
     const int y = i.y();
@@ -303,7 +305,8 @@ void XZHermiteSplineBase<monotonic>::calcWeights(const Field3D& delta_x,
 #endif
 #endif
     if constexpr (monotonic) {
-      const auto gind = gf3daccess->xyzg(i_corn, y + y_offset, k_corner(x, y, z));
+      const auto gind =
+          gf3daccess->xyzg(i_corn, y + y_offset + y_global_offset, k_corner(x, y, z));
       gf3daccess->get(gind);
       gf3daccess->get(gind.xp(1));
       gf3daccess->get(gind.zp(1));

diff --git a/tests/integrated/test_suite b/tests/integrated/test_suite
@@ -188,7 +188,7 @@ class Test(threading.Thread):
                 self.output += "\n(It is likely that a timeout occured)"
             else:
                 # ❌ Failed
-                print("\u274C", end="")  # No newline
+                print("\u274c", end="")  # No newline
             print(" %7.3f s" % (time.time() - self.local.start_time), flush=True)
 
     def _cost(self):