From 9e90e82012fe67c8310318654a65a145acc967aa Mon Sep 17 00:00:00 2001
From: Kristian Evers <kristianevers@gmail.com>
Date: Wed, 5 Dec 2018 17:48:36 +0100
Subject: Allow ellipsoid flattening to be set to zero

Issue raised in https://github.com/OSGeo/proj.4/issues/1191.
---
 src/pj_ell_set.c | 5 ++---
 1 file changed, 2 insertions(+), 3 deletions(-)

(limited to 'src')

diff --git a/src/pj_ell_set.c b/src/pj_ell_set.c
index 6c89ea0d..e2d2750c 100644
--- a/src/pj_ell_set.c
+++ b/src/pj_ell_set.c
@@ -258,9 +258,8 @@ static int ellps_shape (PJ *P) {
         P->f = pj_atof (pj_param_value (par));
         if (HUGE_VAL==P->f)
             return proj_errno_set (P, PJD_ERR_INVALID_ARG);
-        if (0==P->f)
-            return proj_errno_set (P, PJD_ERR_INVALID_ARG);
-        P->rf = 1 / P->f;
+
+        P->rf = P->f != 0.0 ? 1.0/P->f: HUGE_VAL;
         P->es = 2*P->f - P->f*P->f;
         break;
 
-- 
cgit v1.2.3