From 0333785b7d89b3af154793b4f4bb939e1bdaf731 Mon Sep 17 00:00:00 2001
From: Thomas Knudsen <thokn@sdfe.dk>
Date: Tue, 10 Oct 2017 00:35:55 +0200
Subject: Resolve OSS Fuzz issue 3595 division by zero by replacing atan(y/x)
 with atan2(y,x)

---
 src/geocent.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'src')

diff --git a/src/geocent.c b/src/geocent.c
index e340e73d..a8cec043 100644
--- a/src/geocent.c
+++ b/src/geocent.c
@@ -430,7 +430,7 @@ void pj_Convert_Geocentric_To_Geodetic (GeocentricInfo *gi,
     while (SDPHI*SDPHI > genau2 && iter < maxiter);
 
 /*	ellipsoidal (geodetic) latitude */
-    *Latitude=atan(SPHI/fabs(CPHI));
+    *Latitude=atan2(SPHI, fabs(CPHI));
 
     return;
 #endif /* defined(USE_ITERATIVE_METHOD) */
-- 
cgit v1.2.3