From 13247554c9e245093da63145d28bdc965a1d52c4 Mon Sep 17 00:00:00 2001
From: Kristian Evers <kristianevers@gmail.com>
Date: Tue, 23 May 2017 23:31:41 +0200
Subject: Avoid zero division in PJ_eck3.c s_inverser().

Fixes https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=1789

Credit to OSS-Fuzz.
---
 src/PJ_eck3.c | 7 ++++++-
 1 file changed, 6 insertions(+), 1 deletion(-)

(limited to 'src')

diff --git a/src/PJ_eck3.c b/src/PJ_eck3.c
index d70838d2..3fe5c49f 100644
--- a/src/PJ_eck3.c
+++ b/src/PJ_eck3.c
@@ -24,9 +24,14 @@ static XY s_forward (LP lp, PJ *P) {           /* Spheroidal, forward */
 static LP s_inverse (XY xy, PJ *P) {           /* Spheroidal, inverse */
     LP lp = {0.0,0.0};
     struct pj_opaque *Q = P->opaque;
+    double denominator;
 
     lp.phi = xy.y / Q->C_y;
-    lp.lam = xy.x / (Q->C_x * (Q->A + asqrt(1. - Q->B * lp.phi * lp.phi)));
+    denominator = (Q->C_x * (Q->A + asqrt(1. - Q->B * lp.phi * lp.phi)));
+    if ( denominator == 0.0)
+        lp.lam = HUGE_VAL;
+    else
+        lp.lam = xy.x / denominator;
     return lp;
 }
 
-- 
cgit v1.2.3