GEOSHasZ now checks both DoubleNotANumber and isFinite

git-svn-id: http://svn.osgeo.org/postgis/trunk@1099 b70326c6-7e19-0410-871a-916f4a2858ee

diff --git a/lwgeom/lwgeom_geos_wrapper.cpp b/lwgeom/lwgeom_geos_wrapper.cpp
index f5d110d57b3ca326b771279036501ee8f5eba187..c3d6d94646ab814a97ba01718c94923b45b2c2ab 100644 (file)
--- a/lwgeom/lwgeom_geos_wrapper.cpp
+++ b/lwgeom/lwgeom_geos_wrapper.cpp
@@ -1719,5 +1719,5 @@ GEOSHasZ(Geometry *g)
        double az = g->getCoordinate()->z;
        //sprintf(msg, "ZCoord: %g", az);
        //NOTICE_MESSAGE(msg);
-       return (isfinite(az));
+       return (isfinite(az) && ! az == DoubleNotANumber);
 }