diff options
| -rw-r--r-- | src/gauss.cpp | 9 | ||||
| -rw-r--r-- | test/gie/builtins.gie | 3 |
2 files changed, 10 insertions, 2 deletions
diff --git a/src/gauss.cpp b/src/gauss.cpp index 49ccfa1c..a34a8f5b 100644 --- a/src/gauss.cpp +++ b/src/gauss.cpp @@ -65,13 +65,18 @@ void *pj_gauss_ini(double e, double phi0, double *chi, double *rc) { } *chi = asin(sphi / en->C); en->ratexp = 0.5 * en->C * e; + double srat_val = srat(en->e * sphi, en->ratexp); + if (srat_val == 0.0) { + free(en); + return nullptr; + } if( .5 * phi0 + M_FORTPI < 1e-10 ) { - en->K = 1.0 / srat(en->e * sphi, en->ratexp); + en->K = 1.0 / srat_val; } else { en->K = tan(.5 * *chi + M_FORTPI) / ( pow(tan(.5 * phi0 + M_FORTPI), en->C) * - srat(en->e * sphi, en->ratexp) ); + srat_val ); } return ((void *)en); } diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index 2cc44632..91f0cb22 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -4946,6 +4946,9 @@ accept 0 0 expect 0.000000000000 12493602.143489977345 +operation +proj=sterea +a=9999 +b=.9 +lat_0=73 +expect failure + =============================================================================== Gauss-Schreiber Transverse Mercator (aka Gauss-Laborde Reunion) Cyl, Sph&Ell |