Fix for old scipy

Author: flying-sheep

src/fast_array_utils/stats/_generic_ops.py

@@ -90,7 +90,8 @@ def _generic_op_cs(
     # convert to array so dimensions collapse as expected
     x = (sp.csr_array if x.format == "csr" else sp.csc_array)(x, **_dtype_kw(dtype, op))  # type: ignore[arg-type]
     rv = cast("NDArray[Any] | types.coo_array | np.number[Any]", getattr(x, op)(axis=axis))
-    return rv.toarray() if isinstance(rv, types.coo_array) else rv
+    # old scipy versions’ sparray.{max,min}() return a 1×n/n×1 sparray here, so we squeeze
+    return rv.toarray().squeeze() if isinstance(rv, types.coo_array) else rv
 
 
 @generic_op.register(types.DaskArray)

src/fast_array_utils/stats/_utils.py

@@ -110,7 +110,7 @@ def _get_shape(a: NDArray[Any] | np.number[Any] | types.CupyArray, *, axis: Lite
             assert axis is not None
             return (1, a.size) if axis == 0 else (a.size, 1)
         case _:  # pragma: no cover
-            msg = f"{keepdims=}, {type(a)}"
+            msg = f"{keepdims=}, {a.ndim=}, {type(a)=}"
             raise AssertionError(msg)