nlsolve: take the absolute value of delta in the divergence check

Also, reduce the bound from 1e18 to 1e12 which should be sufficient for
most cases.

1 files changed, 1 insertions, 1 deletions

diff --git a/nlsolve.c b/nlsolve.c
index f81e701..d1ab0ff 100644
--- a/nlsolve.c
+++ b/nlsolve.c
@@ -211,7 +211,7 @@ int tdi_nlsolve_solve(NLSolveContext *ctx, NLEqCallback eq_eval,
             ret = 0;
             goto finish;
         } else if ((fast_abort && fabs(s->delta[max_idx]) > 1e6) ||
-                   s->delta[max_idx] > 1e18) {
+                   fabs(s->delta[max_idx]) > 1e12) {
             tdi_log(&ctx->logger, 2, "max(delta) %g, aborting\n", s->delta[max_idx]);
             ret = -EDOM;
             goto finish;