From 3e2b57c676680bbf299cb7bef4dd2eab49148905 Mon Sep 17 00:00:00 2001
From: Peter Hill <peter.hill@york.ac.uk>
Date: Thu, 30 Jul 2026 17:22:51 +0100
Subject: [PATCH] Fix typos in `LaplacePetsc` after conversion

Fixes #3442
---
 .../laplace/impls/petsc/petsc_laplace.cxx      | 18 +++++++++---------
 1 file changed, 9 insertions(+), 9 deletions(-)

diff --git a/src/invert/laplace/impls/petsc/petsc_laplace.cxx b/src/invert/laplace/impls/petsc/petsc_laplace.cxx
index 33481f9bde..bc1cd91d51 100644
--- a/src/invert/laplace/impls/petsc/petsc_laplace.cxx
+++ b/src/invert/laplace/impls/petsc/petsc_laplace.cxx
@@ -336,7 +336,7 @@ FieldPerp LaplacePetsc::solve(const FieldPerp& b, const FieldPerp& x0) {
   }
   if (!isOuterBoundaryFlagSet(INVERT_RHS)) {
     BOUT_FOR_SERIAL(index, indexer->getRegionOuterX()) {
-      rhs(index) = isInnerBoundaryFlagSet(INVERT_SET) ? x0[index] : 0.0;
+      rhs(index) = isOuterBoundaryFlagSet(INVERT_SET) ? x0[index] : 0.0;
     }
   }
 
@@ -488,10 +488,10 @@ void LaplacePetsc::setSecondOrderMatrix(int y, bool inner_X_neumann,
   // Set the boundaries
   if (inner_X_neumann) {
     const auto dx = sliceXZ(coords->dx, y);
-    const auto g11 = sliceXZ(coords->g11, y);
+    const auto g_11 = sliceXZ(coords->g_11, y);
 
     BOUT_FOR_SERIAL(i, indexer->getRegionInnerX()) {
-      const auto factor = 1. / dx[i] / std::sqrt(g11[i]);
+      const auto factor = 1. / dx[i] / std::sqrt(g_11[i]);
       operator2D(i, i) = -factor;
       operator2D(i, i.xp()) = factor;
     }
@@ -503,10 +503,10 @@ void LaplacePetsc::setSecondOrderMatrix(int y, bool inner_X_neumann,
   }
   if (outer_X_neumann) {
     const auto dx = sliceXZ(coords->dx, y);
-    const auto g11 = sliceXZ(coords->g11, y);
+    const auto g_11 = sliceXZ(coords->g_11, y);
 
     BOUT_FOR_SERIAL(i, indexer->getRegionOuterX()) {
-      const auto factor = 1. / dx[i] / std::sqrt(g11[i]);
+      const auto factor = 1. / dx[i] / std::sqrt(g_11[i]);
       operator2D(i, i) = factor;
       operator2D(i, i.xm()) = -factor;
     }
@@ -558,10 +558,10 @@ void LaplacePetsc::setFourthOrderMatrix(int y, bool inner_X_neumann,
   // Set boundaries
   if (inner_X_neumann) {
     const auto dx = sliceXZ(coords->dx, y);
-    const auto g11 = sliceXZ(coords->g11, y);
+    const auto g_11 = sliceXZ(coords->g_11, y);
 
     BOUT_FOR_SERIAL(i, indexer->getRegionInnerX()) {
-      const auto factor = 1. / dx[i] / std::sqrt(g11[i]);
+      const auto factor = 1. / dx[i] / std::sqrt(g_11[i]);
       operator2D(i, i) = (-25.0 / 12.0) * factor;
       operator2D(i, i.xp(1)) = 4.0 * factor;
       operator2D(i, i.xp(2)) = -3.0 * factor;
@@ -580,10 +580,10 @@ void LaplacePetsc::setFourthOrderMatrix(int y, bool inner_X_neumann,
 
   if (outer_X_neumann) {
     const auto dx = sliceXZ(coords->dx, y);
-    const auto g11 = sliceXZ(coords->g11, y);
+    const auto g_11 = sliceXZ(coords->g_11, y);
 
     BOUT_FOR_SERIAL(i, indexer->getRegionOuterX()) {
-      const auto factor = 1. / dx[i] / std::sqrt(g11[i]);
+      const auto factor = 1. / dx[i] / std::sqrt(g_11[i]);
       operator2D(i, i) = (25.0 / 12.0) * factor;
       operator2D(i, i.xm(1)) = -4.0 * factor;
       operator2D(i, i.xm(2)) = 3.0 * factor;