From d85d78fc4ff716868db835b7d920637064c84573 Mon Sep 17 00:00:00 2001
From: Julian Litz <91479202+julianlitz@users.noreply.github.com>
Date: Sun, 17 May 2026 19:16:09 +0200
Subject: [PATCH 1/3] Test Issue #31

---
 src/PolarGrid/polargrid.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/PolarGrid/polargrid.cpp b/src/PolarGrid/polargrid.cpp
index a54aa3645..70b6f2aa9 100644
--- a/src/PolarGrid/polargrid.cpp
+++ b/src/PolarGrid/polargrid.cpp
@@ -211,7 +211,7 @@ void PolarGrid::initializeLineSplitting(std::optional<double> splitting_radius)
              i_r++) { /* We assume lengthRadialSmoother_ >= 3 in the further implementation */
             double uniform_theta_k = (2 * M_PI) / ntheta();
             double radius_r        = radius(i_r);
-            double radial_dist_h   = radius(i_r) - radius(i_r - 1);
+            double radial_dist_h   = radius(i_r + 1) - radius(i_r);
             ;
             double q = uniform_theta_k / radial_dist_h;
             if (q * radius_r > 1.0) {

From f7f353dc9bdb9651685b9d75f770c5f9fd11efad Mon Sep 17 00:00:00 2001
From: Julian Litz <91479202+julianlitz@users.noreply.github.com>
Date: Sun, 17 May 2026 19:45:32 +0200
Subject: [PATCH 2/3] Update solve_tests.cpp

---
 tests/GMGPolar/solve_tests.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/tests/GMGPolar/solve_tests.cpp b/tests/GMGPolar/solve_tests.cpp
index bc2472977..5e6e91502 100644
--- a/tests/GMGPolar/solve_tests.cpp
+++ b/tests/GMGPolar/solve_tests.cpp
@@ -301,7 +301,7 @@ using gmgpolar_test_suite = testing::Types<
         std::integral_constant<ResidualNormType, ResidualNormType::INFINITY_NORM>, // residualNormType
         std::integral_constant<double, 1e-12>, // absoluteTolerance
         std::integral_constant<double, 1e-10>, // relativeTolerance
-        std::integral_constant<int, 13>, // expected_iterations
+        std::integral_constant<int, 7>, // expected_iterations
         std::integral_constant<double, 6e-6>, // expected_l2_error
         std::integral_constant<double, 2e-5>, // expected_inf_error
         std::integral_constant<double, 0.3> // expected_residual_reduction
@@ -508,7 +508,7 @@ using gmgpolar_test_suite = testing::Types<
         std::integral_constant<ResidualNormType, ResidualNormType::EUCLIDEAN>, // residualNormType
         std::integral_constant<double, 1e-9>, // absoluteTolerance
         std::integral_constant<double, 1e-8>, // relativeTolerance
-        std::integral_constant<int, 8>, // expected_iterations
+        std::integral_constant<int, 6>, // expected_iterations
         std::integral_constant<double, 5e-6>, // expected_l2_error
         std::integral_constant<double, 2e-5>, // expected_inf_error
         std::integral_constant<double, 0.2> // expected_residual_reduction

From 4dc3ec0bbc813ae6f5f5bcbe3fad629ec5492c58 Mon Sep 17 00:00:00 2001
From: Julian Litz <91479202+julianlitz@users.noreply.github.com>
Date: Sun, 17 May 2026 19:46:14 +0200
Subject: [PATCH 3/3] Update pcg_tests.cpp

---
 tests/GMGPolar/pcg_tests.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/tests/GMGPolar/pcg_tests.cpp b/tests/GMGPolar/pcg_tests.cpp
index d03cdbe4a..336c31ff1 100644
--- a/tests/GMGPolar/pcg_tests.cpp
+++ b/tests/GMGPolar/pcg_tests.cpp
@@ -222,7 +222,7 @@ using gmgpolar_test_suite = testing::Types<
         std::integral_constant<ResidualNormType, ResidualNormType::EUCLIDEAN>, // residualNormType
         std::integral_constant<double, 1e-8>, // absoluteTolerance
         std::integral_constant<double, 1e-8>, // relativeTolerance
-        std::integral_constant<int, 8>, // expected_iterations
+        std::integral_constant<int, 7>, // expected_iterations
         std::integral_constant<double, 2e-6>, // expected_l2_error
         std::integral_constant<double, 9e-6>, // expected_inf_error
         std::integral_constant<double, 0.7> // expected_residual_reduction
@@ -349,7 +349,7 @@ using gmgpolar_test_suite = testing::Types<
         std::integral_constant<ResidualNormType, ResidualNormType::INFINITY_NORM>, // residualNormType
         std::integral_constant<double, 1e-12>, // absoluteTolerance
         std::integral_constant<double, 1e-10>, // relativeTolerance
-        std::integral_constant<int, 4>, // expected_iterations
+        std::integral_constant<int, 3>, // expected_iterations
         std::integral_constant<double, 6e-6>, // expected_l2_error
         std::integral_constant<double, 2e-5>, // expected_inf_error
         std::integral_constant<double, 0.3> // expected_residual_reduction