From dea2d656c4ca658239cb3d4d300c47b2477eb2ba Mon Sep 17 00:00:00 2001
From: Frank Warmerdam <warmerdam@pobox.com>
Date: Mon, 21 Oct 2013 20:56:49 +0000
Subject: Change handling of values nearly 90degrees away from the origin
 (#114)

git-svn-id: http://svn.osgeo.org/metacrs/proj/trunk@2400 4e78687f-474d-0410-85f9-8d5e500ac6b2
---
 src/PJ_omerc.c | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

(limited to 'src')

diff --git a/src/PJ_omerc.c b/src/PJ_omerc.c
index c8888a83..3529b73c 100644
--- a/src/PJ_omerc.c
+++ b/src/PJ_omerc.c
@@ -48,8 +48,11 @@ FORWARD(e_forward); /* ellipsoid */
 			F_ERROR;
 		v = 0.5 * P->ArB * log((1. - U)/(1. + U));
 		temp = cos(P->B * lp.lam);
-		u = (fabs(temp) < TOL) ? P->AB * lp.lam :
-			P->ArB * atan2((S * P->cosgam + V * P->singam) , temp); 
+                if(fabs(temp) < TOL) {
+                    u = P->A * lp.lam;
+                } else {
+                    u = P->ArB * atan2((S * P->cosgam + V * P->singam), temp);
+                }
 	} else {
 		v = lp.phi > 0 ? P->v_pole_n : P->v_pole_s;
 		u = P->ArB * lp.phi;
-- 
cgit v1.2.3