From 9b561406e1ddcde02e868bef908e5cd5140745c7 Mon Sep 17 00:00:00 2001
From: Even Rouault <even.rouault@spatialys.com>
Date: Thu, 3 Oct 2019 20:48:06 +0200
Subject: aeqd: for spherical forward path, go to higher precision ellipsoidal
 case when the point coordinates are super close to the origin (fixes #1654)

---
 test/gie/builtins.gie | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)

(limited to 'test')

diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie
index f86b7e79..676086e0 100644
--- a/test/gie/builtins.gie
+++ b/test/gie/builtins.gie
@@ -136,6 +136,29 @@ accept  45  45
 expect  3_860_398.3783  5_430_089.0490
 roundtrip   100
 
+
+# Test oblique aeqd with point very close lon_0, lat_0, on a perfect sphere
+operation +proj=aeqd +a=6371008.771415 +b=6371008.771415 +lat_0=30.2345 +lon_0=-120.2345
+tolerance 1 mm
+accept    -120.234501 30.234501
+expect    -0.096      0.111
+roundtrip 1
+
+accept    -120.2345   30.2345
+expect    0.000       0.000
+roundtrip 1
+
+# Same on an ellipsoid very close to the sphere
+operation +proj=aeqd +a=6371008.771415 +b=6371008.771414 +lat_0=30.2345 +lon_0=-120.2345
+tolerance 1 mm
+accept    -120.234501 30.234501
+expect    -0.096      0.111
+roundtrip 1
+
+accept    -120.2345   30.2345
+expect    0.000       0.000
+roundtrip 1
+
 -------------------------------------------------------------------------------
 # Test the Modified Azimuthal Equidistant / EPSG 9832. Test data from the EPSG
 Guidance Note 7 part 2, April 2018,  p. 85
-- 
cgit v1.2.3