Reset of poly_indexs due to an unknown segmentation fault

Author: AntSimi

src/py_eddy_tracker/poly.py

@@ -815,7 +815,7 @@ def box_indexes(x, y, step):
 
 
 @njit(cache=True)
-def poly_indexs(x_p, y_p, x_c, y_c):
+def poly_indexs_(x_p, y_p, x_c, y_c):
     """
     Index of contour for each postion inside a contour, -1 in case of no contour
 
@@ -873,6 +873,34 @@ def poly_indexs(x_p, y_p, x_c, y_c):
     return indexs
 
 
+@njit(cache=True)
+def poly_indexs(x_p, y_p, x_c, y_c):
+    """
+    index of contour for each postion inside a contour, -1 in case of no contour
+
+    :param array x_p: longitude to test
+    :param array y_p: latitude to test
+    :param array x_c: longitude of contours
+    :param array y_c: latitude of contours
+    """
+    nb_p = x_p.shape[0]
+    nb_c = x_c.shape[0]
+    indexs = -ones(nb_p, dtype=numba_types.int32)
+    for i in range(nb_c):
+        x_, y_ = reduce_size(x_c[i], y_c[i])
+        x_c_min, y_c_min = x_.min(), y_.min()
+        x_c_max, y_c_max = x_.max(), y_.max()
+        v = create_vertice(x_, y_)
+        for j in range(nb_p):
+            if indexs[j] != -1:
+                continue
+            x, y = x_p[j], y_p[j]
+            if x > x_c_min and x < x_c_max and y > y_c_min and y < y_c_max:
+                if winding_number_poly(x, y, v) != 0:
+                    indexs[j] = i
+    return indexs
+
+
 @njit(cache=True)
 def insidepoly(x_p, y_p, x_c, y_c):
     """