fix point_in_polygon method when the asum was a multiple of 2*phi (#4495

)

Co-authored-by: maxcapodi78 <Shark78>
Co-authored-by: adimaria <[email protected]>

pyaedt/modeler/geometry_operators.py

@@ -1558,6 +1558,8 @@ def point_in_polygon(point, polygon, tolerance=1e-8):
 
         The method implements the radial algorithm (https://es.wikipedia.org/wiki/Algoritmo_radial)
 
+        This version supports also self-intersecting polygons.
+
         point : List
             List of ``[x, y]`` coordinates.
         polygon : List
@@ -1591,9 +1593,10 @@ def point_in_polygon(point, polygon, tolerance=1e-8):
             if abs(abs(a) - math.pi) < tol:
                 return 0
             asum += a
+        r = asum % (2*math.pi)
         if abs(asum) < tol:
             return -1
-        elif abs(abs(asum) - 2*math.pi) < tol:
+        elif r < tol or (2*math.pi - r) < tol:
             return 1
         else:  # pragma: no cover
             raise Exception("Unexpected error!")