Fix np.unique() calls for numpy 2.0

Numpy now returns a 2D array for reverse indexing, but we want
the flattened version.

Author: xylar

jigsawpy/bisect.py

@@ -208,6 +208,9 @@ def bisect(mesh):
             return_index=True,
             return_inverse=True, axis=0)
 
+        # for numpy 2.x compatibility
+        erev = erev.ravel()
+
         edge = edge[efwd, :]
 
     if (quad.size != 0):
@@ -216,6 +219,9 @@ def bisect(mesh):
             return_index=True,
             return_inverse=True, axis=0)
 
+        # for numpy 2.x compatibility
+        qrev = qrev.ravel()
+
         quad = quad[qfwd, :]
 
     if (hexa.size != 0):
@@ -224,6 +230,9 @@ def bisect(mesh):
             return_index=True,
             return_inverse=True, axis=0)
 
+        # for numpy 2.x compatibility
+        hrev = hrev.ravel()
+
         hexa = hexa[hfwd, :]
 
 #------------------------------ map unique joins to elements