From 17f2f7cf8bcaa5a4edc9e94d2bd6d8e633455c03 Mon Sep 17 00:00:00 2001
From: Even Rouault <even.rouault@spatialys.com>
Date: Wed, 20 Mar 2019 22:22:38 +0100
Subject: lcc: avoid division by zero

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

(limited to 'src')

diff --git a/src/projections/lcc.cpp b/src/projections/lcc.cpp
index a1fe79a9..5eee0d14 100644
--- a/src/projections/lcc.cpp
+++ b/src/projections/lcc.cpp
@@ -108,6 +108,10 @@ PJ *PROJECTION(lcc) {
         if (secant) { /* secant cone */
             sinphi = sin(Q->phi2);
             Q->n = log(m1 / pj_msfn(sinphi, cos(Q->phi2), P->es));
+            if (Q->n == 0) {
+                // Not quite, but es is very close to 1...
+                return pj_default_destructor(P, PJD_ERR_ECCENTRICITY_IS_ONE_OR_GREATER);
+            }
             Q->n /= log(ml1 / pj_tsfn(Q->phi2, sinphi, P->e));
         }
         Q->c = (Q->rho0 = m1 * pow(ml1, -Q->n) / Q->n);
-- 
cgit v1.2.3