From ff1e495525969352f83507950dd4e1fd3cec69c2 Mon Sep 17 00:00:00 2001
From: Even Rouault <even.rouault@spatialys.com>
Date: Fri, 29 Mar 2019 23:09:16 +0100
Subject: tpeqd: avoid division by zero

Fixes https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=13948
Credit to OSS Fuzz
---
 src/projections/tpeqd.cpp | 4 ++++
 1 file changed, 4 insertions(+)

(limited to 'src')

diff --git a/src/projections/tpeqd.cpp b/src/projections/tpeqd.cpp
index 20921de4..9904bb8c 100644
--- a/src/projections/tpeqd.cpp
+++ b/src/projections/tpeqd.cpp
@@ -87,6 +87,10 @@ PJ *PROJECTION(tpeqd) {
     Q->sc = Q->sp1 * Q->cp2;
     Q->ccs = Q->cp1 * Q->cp2 * sin(Q->dlam2);
     Q->z02 = aacos(P->ctx, Q->sp1 * Q->sp2 + Q->cp1 * Q->cp2 * cos (Q->dlam2));
+    if( Q->z02 == 0.0 ) {
+        // Actually happens when both lat_1 = lat_2 and |lat_1| = 90
+        return pj_default_destructor(P, PJD_ERR_LAT_1_OR_2_ZERO_OR_90);
+    }
     Q->hz0 = .5 * Q->z02;
     A12 = atan2(Q->cp2 * sin (Q->dlam2),
         Q->cp1 * Q->sp2 - Q->sp1 * Q->cp2 * cos (Q->dlam2));
-- 
cgit v1.2.3