Merge pull request #59 from GComitini/integrate-early-exit

Add Nan/infinite guard to gauss_kronrod_quadrature (early exit)
Axect · Apr 25, 2024 · 6804677 · 6804677
2 parents 690c2e9 + 169ab6d
commit 6804677
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));