Merge pull request #11421 from livace/fix_point_in_polygon_with_nan

Fix pointInPolygon with nan

(cherry picked from commit 42534959)

src/Functions/GeoUtils.h

@@ -301,6 +301,9 @@ bool PointInPolygonWithGrid<CoordinateType>::contains(CoordinateType x, Coordina
     if (has_empty_bound)
         return false;
 
+    if (std::isnan(x) || std::isnan(y))
+        return false;
+
     CoordinateType float_row = (y + y_shift) * y_scale;
     CoordinateType float_col = (x + x_shift) * x_scale;
 

tests/queries/0_stateless/00500_point_in_polygon_nan.reference

+0

tests/queries/0_stateless/00500_point_in_polygon_nan.sql

+SELECT pointInPolygon((nan, 10.000100135803223), [(39.83154, 21.41527), (2., 1000.0001220703125), (39.90033, 21.37195), (1.000100016593933, 10.000100135803223), (39.83051, 21.42553), (39.82898, 21.41382), (39.83043, 21.41432), (39.83154, 21.41527)]);