Fixing bug affecting Q-type interpolation at the boundaries of the in…

…terpolation grid. Thanks to S. Rodini for spotting it.)
vbertone · Sep 3, 2024 · 58597f5 · 58597f5
1 parent 8cfd62f
commit 58597f5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/kernel/qgrid.cc b/src/kernel/qgrid.cc
@@ -256,7 +256,7 @@ namespace apfel
 
     // Return if "Q" is outside the grid range (no sum will be
     // performed).
-    if (Q < _Qg.front() - eps12 || Q > _Qg.back() + eps12)
+    if (Q < _Qg.front() || Q > _Qg.back())
       return bounds;
 
     // If Q falls in the tiny gap at the thresholds assume the node