From 58597f501717812f64b5c7ac35942d3d0a267432 Mon Sep 17 00:00:00 2001 From: vbertone Date: Tue, 3 Sep 2024 09:42:21 +0200 Subject: [PATCH] Fixing bug affecting Q-type interpolation at the boundaries of the interpolation grid. Thanks to S. Rodini for spotting it.) --- src/kernel/qgrid.cc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/kernel/qgrid.cc b/src/kernel/qgrid.cc index 836e5e156..930effc4d 100644 --- 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