Increase MAX_ITER so Mollweide forward projection works near the pole…

…s. (#3082)
OSGeo · Mar 5, 2022 · 8a88101 · 8a88101
1 parent c34b0ac
commit 8a88101
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 1 deletion.
diff --git a/src/projections/moll.cpp b/src/projections/moll.cpp
@@ -10,7 +10,7 @@ PROJ_HEAD(moll, "Mollweide") "\n\tPCyl, Sph";
 PROJ_HEAD(wag4, "Wagner IV") "\n\tPCyl, Sph";
 PROJ_HEAD(wag5, "Wagner V") "\n\tPCyl, Sph";
 
-#define MAX_ITER    10
+#define MAX_ITER    30
 #define LOOP_TOL    1e-7
 
 namespace { // anonymous namespace

diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie
@@ -3599,6 +3599,18 @@ roundtrip 1
 accept  -2 -1
 expect  -201113.698641813 -124066.283433860
 roundtrip 1
+accept  0.0 89.99
+expect  0.00 9050917.562466157600
+roundtrip 1
+accept  0.0 89.999
+expect  0.00 9050964.513822982088
+roundtrip 1
+accept  0.0 -89.99
+expect  0.0 -9050917.562466157600
+roundtrip 1
+accept  0.0 -89.999
+expect  0.00 -9050964.513822982088
+roundtrip 1
 
 direction inverse
 accept  200 100