Remove in-place squeezes

Author: lshaw8317

src/blosc2/linalg.py

@@ -127,9 +127,9 @@ def matmul(x1: blosc2.Array, x2: blosc2.NDArray, **kwargs: Any) -> blosc2.NDArra
                         result[chunk + (slice(row, row_end), slice(col, col_end))] += np.matmul(bx1, bx2)
 
     if x1_is_vector:
-        result.squeeze(axis=-2)
+        result = result.squeeze(axis=-2)
     if x2_is_vector:
-        result.squeeze(axis=-1)
+        result = result.squeeze(axis=-1)
 
     return result
 

src/blosc2/ndarray.py

@@ -4681,8 +4681,7 @@ def slice(self, key: int | slice | Sequence[slice], **kwargs: Any) -> NDArray:
                 for order, nchunk in enumerate(aligned_chunks):
                     chunk = self.schunk.get_chunk(nchunk)
                     newarr.schunk.update_chunk(order, chunk)
-                newarr.squeeze(axis=np.where(mask)[0])  # remove any dummy dims introduced
-                return newarr
+                return newarr.squeeze(axis=np.where(mask)[0])  # remove any dummy dims introduced
 
         key = (start, stop)
         ndslice = super().get_slice(key, mask, **kwargs)
@@ -4724,16 +4723,7 @@ def squeeze(self, axis: int | Sequence[int]) -> NDArray:
         >>> a.shape
         (23, 11)
         """
-        axis = [axis] if isinstance(axis, int) else axis
-        mask = [False for i in range(self.ndim)]
-        for a in axis:
-            if a < 0:
-                a += self.ndim  # Adjust axis to be within the array's dimensions
-            if mask[a]:
-                raise ValueError("Axis values must be unique.")
-            mask[a] = True
-        super().squeeze(axis_mask=mask)
-        return self
+        return blosc2.squeeze(self, axis=axis)
 
     def indices(self, order: str | list[str] | None = None, **kwargs: Any) -> NDArray:
         """