diff options
| author | Kristian Evers <kristianevers@gmail.com> | 2017-05-23 23:31:41 +0200 |
|---|---|---|
| committer | Kristian Evers <kristianevers@gmail.com> | 2017-05-24 11:20:09 +0200 |
| commit | 13247554c9e245093da63145d28bdc965a1d52c4 (patch) | |
| tree | 0a5bfa3dbc701fb8889dc52c1341dda35a05e1c5 | |
| parent | 9e402c2577e702dd0d45df35502dea08b774e072 (diff) | |
| download | PROJ-13247554c9e245093da63145d28bdc965a1d52c4.tar.gz PROJ-13247554c9e245093da63145d28bdc965a1d52c4.zip | |
Avoid zero division in PJ_eck3.c s_inverser().
Fixes https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=1789
Credit to OSS-Fuzz.
| -rw-r--r-- | src/PJ_eck3.c | 7 |
1 files changed, 6 insertions, 1 deletions
diff --git a/src/PJ_eck3.c b/src/PJ_eck3.c index d70838d2..3fe5c49f 100644 --- a/src/PJ_eck3.c +++ b/src/PJ_eck3.c @@ -24,9 +24,14 @@ static XY s_forward (LP lp, PJ *P) { /* Spheroidal, forward */ static LP s_inverse (XY xy, PJ *P) { /* Spheroidal, inverse */ LP lp = {0.0,0.0}; struct pj_opaque *Q = P->opaque; + double denominator; lp.phi = xy.y / Q->C_y; - lp.lam = xy.x / (Q->C_x * (Q->A + asqrt(1. - Q->B * lp.phi * lp.phi))); + denominator = (Q->C_x * (Q->A + asqrt(1. - Q->B * lp.phi * lp.phi))); + if ( denominator == 0.0) + lp.lam = HUGE_VAL; + else + lp.lam = xy.x / denominator; return lp; } |