IMPL: add Nan/infinite guard to gauss_kronrod_quadrature (early exit)

Axect · Apr 23, 2024 · 169ab6d · 169ab6d
1 parent dea021a
commit 169ab6d
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/numerical/integral.rs b/src/numerical/integral.rs
@@ -240,6 +240,9 @@ where
                     tol
                 };
                 if (G - K).abs() < tol_curr || a == b || max_iter == 0 {
+                    if ! G.is_finite() {
+                        return G;
+                    }
                     I += G;
                 } else {
                     S.push((a, c, tol / 2f64, max_iter - 1));