DRAFT: Use contiguous arrays for point containment

uxarray/grid/coordinates.py

@@ -757,4 +757,4 @@ def points_atleast_2d_xyz(points):
             "Points are neither Cartesian (shape N x 3) nor Spherical (shape N x 2)."
         )
 
-    return points_xyz
+    return np.ascontiguousarray(points_xyz)

uxarray/grid/point_in_face.py

@@ -205,6 +205,7 @@ def _point_in_face_query(
         Number of valid face‐hits per point.
     """
     pts = np.asarray(points, dtype=np.float64)
+    pts = np.ascontiguousarray(pts)
     if pts.ndim == 1:
         pts = pts[np.newaxis, :]
     # Cull with k-d tree

uxarray/plot/matplotlib.py

@@ -144,7 +144,10 @@ def _get_points_from_axis(ax: GeoAxes, *, pixel_ratio: float = 1):
     xs = np.cos(lats) * np.cos(lons)
     ys = np.cos(lats) * np.sin(lons)
     zs = np.sin(lats)
-    pts = np.vstack((xs, ys, zs)).T
+
+    # TODO:
+    pts = np.column_stack((xs, ys, zs))
+    # pts = np.vstack((xs, ys, zs)).T
 
     return pts, valid, nx, ny