From d5c2e2b6e5b890f5548e4a79e60c38ecb44fce14 Mon Sep 17 00:00:00 2001
From: Thomas Ceulemans <thomas.ceulemans@kuleuven.be>
Date: Wed, 16 Oct 2024 15:58:06 +0200
Subject: [PATCH] Added boundary treatment fix for 1D imager.

---
 src/solver/solver.tpp | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/src/solver/solver.tpp b/src/solver/solver.tpp
index b1c0d644..f4a60b0f 100644
--- a/src/solver/solver.tpp
+++ b/src/solver/solver.tpp
@@ -888,7 +888,9 @@ accel inline Size Solver ::trace_ray_imaging(const Model& model, const Vector3D&
         }
     }
 
-    if (!accessed_non_boundary) {
+    // Spherically symmetric models have a correct boundary treatment already, thus no need to
+    // reduce their accuracy in a region between the outer boundary and the nearest point.
+    if (!accessed_non_boundary && !model.parameters->spherical_symmetry()) {
         return original_id1;
     }
     return id1;