Merge pull request #1658 from OSGeo/backport-1655-to-6.2

[Backport 6.2] aeqd: for spherical forward path, go to higher precision ellipsoidal case when the point coordinates are super close to the origin (fixes #1654)

1 files changed, 1 insertions, 1 deletions

diff --git a/src/projections/aeqd.cpp b/src/projections/aeqd.cpp
index 04c3662e..eb38cbf3 100644
--- a/src/projections/aeqd.cpp
+++ b/src/projections/aeqd.cpp
@@ -151,7 +151,7 @@ oblcon:
                 return xy;
             }
             else
-                xy.x = xy.y = 0.;
+                return aeqd_e_forward(lp, P);
         else {
             xy.y = acos(xy.y);
             xy.y /= sin(xy.y);