diff options
| author | Even Rouault <even.rouault@spatialys.com> | 2021-04-03 20:17:15 +0200 |
|---|---|---|
| committer | Even Rouault <even.rouault@spatialys.com> | 2021-04-03 21:59:30 +0200 |
| commit | 76973c9e567e815053b6dc9e71b9d93e0e9b9964 (patch) | |
| tree | e8177c84a7af49aa1cd6d9d5af649d9d1203ecf4 | |
| parent | 734b69804386dda2181933bf8377f45f7105aeb1 (diff) | |
| download | PROJ-76973c9e567e815053b6dc9e71b9d93e0e9b9964.tar.gz PROJ-76973c9e567e815053b6dc9e71b9d93e0e9b9964.zip | |
Generic inverse: fix when initial guess if super close to result
| -rw-r--r-- | src/generic_inverse.cpp | 2 |
1 files changed, 1 insertions, 1 deletions
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; |