From 0529b07f81d3c027e101c6e1eddb4685e957934d Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Sun, 24 Mar 2019 17:05:11 +0100 Subject: tmerc inverse: avoid division by zero Fixes https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=13893 Credit to OSS Fuzz --- src/projections/tmerc.cpp | 4 ++++ test/gie/builtins.gie | 8 ++++++++ 2 files changed, 12 insertions(+) diff --git a/src/projections/tmerc.cpp b/src/projections/tmerc.cpp index c91c5174..bb56f8ae 100644 --- a/src/projections/tmerc.cpp +++ b/src/projections/tmerc.cpp @@ -188,6 +188,10 @@ static PJ_LP approx_s_inv (PJ_XY xy, PJ *P) { double h, g; h = exp(xy.x / static_cast(P->opaque)->esp); + if( h == 0 ) { + proj_errno_set(P, PJD_ERR_INVALID_X_OR_Y); + return proj_coord_error().lp; + } g = .5 * (h - 1. / h); h = cos (P->phi0 + xy.y / static_cast(P->opaque)->esp); lp.phi = asin(sqrt((1. - h * h) / (1. + g * g))); diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index a28ff45c..f0290fa2 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -5058,6 +5058,14 @@ expect -0.001790493 0.000895247 accept -200 -100 expect -0.001790493 -0.000895247 + +------------------------------------------------------------------------------- +operation +proj=tmerc +R=1 +------------------------------------------------------------------------------- +direction inverse +accept -1e200 0 +expect failure errno invalid_x_or_y + =============================================================================== Tobler-Mercator Cyl, Sph -- cgit v1.2.3