From 80b9f98f2316b1d36f6409861cab6bcdddf23174 Mon Sep 17 00:00:00 2001
From: Anton Khirnov <anton@khirnov.net>
Date: Thu, 30 Jan 2020 13:22:46 +0100
Subject: 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.
---
 nlsolve.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

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;
-- 
cgit v1.2.3