Set region for lagrange4pt

dschwoerer · dschwoerer · commit ebf2d07d071b · 2025-01-16T09:07:51.000+01:00
diff --git a/src/mesh/interpolation/lagrange_4pt_xz.cxx b/src/mesh/interpolation/lagrange_4pt_xz.cxx
@@ -132,7 +132,11 @@ Field3D XZLagrange4pt::interpolate(const Field3D& f, const std::string& region)
 
     // Then in X
     f_interp(x, y_next, z) = lagrange_4pt(xvals, t_x(x, y, z));
+    ASSERT2(std::isfinite(f_interp(x, y_next, z)));
+
   }
+  const auto region2 = y_offset != 0 ? fmt::format("RGN_YPAR_{:+d}", y_offset) : region;
+  f_interp.setRegion(region2);
   return f_interp;
 }
 

Original file line number	Diff line number	Diff line change
`@@ -132,7 +132,11 @@ Field3D XZLagrange4pt::interpolate(const Field3D& f, const std::string& region)`
`132`	`132`
`133`	`133`	`// Then in X`
`134`	`134`	`f_interp(x, y_next, z) = lagrange_4pt(xvals, t_x(x, y, z));`
	`135`	`+ ASSERT2(std::isfinite(f_interp(x, y_next, z)));`
	`136`	`+`
`135`	`137`	`}`
	`138`	`+ const auto region2 = y_offset != 0 ? fmt::format("RGN_YPAR_{:+d}", y_offset) : region;`
	`139`	`+ f_interp.setRegion(region2);`
`136`	`140`	`return f_interp;`
`137`	`141`	`}`
`138`	`142`