From 2d0acd7d61374a167f532cd2851f46f63d5e09a4 Mon Sep 17 00:00:00 2001 From: Howard Butler Date: Sat, 13 Sep 2014 20:19:37 +0000 Subject: Apply patch in #244 along with test scenario git-svn-id: http://svn.osgeo.org/metacrs/proj/trunk@2520 4e78687f-474d-0410-85f9-8d5e500ac6b2 --- nad/testvarious | 18 ++++++++++++++++++ src/PJ_aeqd.c | 2 +- 2 files changed, 19 insertions(+), 1 deletion(-) diff --git a/nad/testvarious b/nad/testvarious index a1d11987..fded86ea 100755 --- a/nad/testvarious +++ b/nad/testvarious @@ -663,6 +663,24 @@ $EXE -f '%.13f' \ -E >>${OUT} <> ${OUT} +echo "Test bug 244 " >> ${OUT} +# +$EXE -f '%.12f' \ + +init=epsg:4326 \ + +to +proj=aeqd +lon_0=130.0 +lat_0=40.0 +a=6378137 +b=6378137 +units=m +no_defs \ + -E >>${OUT} <> ${OUT} +echo "Test bug 244 (2)" >> ${OUT} +# +$EXE -f '%.12f' \ + +proj=aeqd +lon_0=130.0 +lat_0=40.0 +a=6378137 +b=6378137 +units=m +no_defs \ + +to +init=epsg:4326 \ + -E >>${OUT} <sinph0 * sin(lp.phi)) * c_rh; xy.x *= sinc * P->cosph0; } - lp.lam = xy.y == 0. ? 0. : atan2(xy.x, xy.y); + lp.lam = atan2(xy.x, xy.y); } else if (P->mode == N_POLE) { lp.phi = HALFPI - c_rh; lp.lam = atan2(xy.x, -xy.y); -- cgit v1.2.3