From 76973c9e567e815053b6dc9e71b9d93e0e9b9964 Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Sat, 3 Apr 2021 20:17:15 +0200 Subject: Generic inverse: fix when initial guess if super close to result --- src/generic_inverse.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/generic_inverse.cpp b/src/generic_inverse.cpp index 41cf6c50..28d939fb 100644 --- a/src/generic_inverse.cpp +++ b/src/generic_inverse.cpp @@ -55,7 +55,7 @@ PJ_LP pj_generic_inverse_2d(PJ_XY xy, PJ *P, PJ_LP lpInitial) { return lp; } - if (fabs(deltaX) > 1e-6 || fabs(deltaY) > 1e-6) { + if (i == 0 || fabs(deltaX) > 1e-6 || fabs(deltaY) > 1e-6) { // Compute Jacobian matrix (only if we aren't close to the final // result to speed things a bit) PJ_LP lp2; -- cgit v1.2.3