From 9e402c2577e702dd0d45df35502dea08b774e072 Mon Sep 17 00:00:00 2001
From: Kristian Evers <kristianevers@gmail.com>
Date: Tue, 23 May 2017 23:21:09 +0200
Subject: Avoid zero division in PJ_aea.c setup().

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

Credit to OSS-Fuzz.
---
 src/PJ_aea.c | 2 ++
 1 file changed, 2 insertions(+)

(limited to 'src')

diff --git a/src/PJ_aea.c b/src/PJ_aea.c
index e6bb413c..524657b2 100644
--- a/src/PJ_aea.c
+++ b/src/PJ_aea.c
@@ -165,6 +165,8 @@ static PJ *setup(PJ *P) {
             cosphi = cos(Q->phi2);
             m2 = pj_msfn(sinphi, cosphi, P->es);
             ml2 = pj_qsfn(sinphi, P->e, P->one_es);
+            if (ml2 == ml1)
+                return NULL;
             Q->n = (m1 * m1 - m2 * m2) / (ml2 - ml1);
         }
         Q->ec = 1. - .5 * P->one_es * log((1. - P->e) /
-- 
cgit v1.2.3