I'm trying to convert a huge scipy csc_matrix (27998 x 1306127, int32) to a persistent zarr array. Here is the code I'm using: However, this fails because scipy cannot convert sparse matrices with ...

It appears that adding or subtracting numpy.ndarrays with scipy.sparse matrices returns a numpy.matrix. Is the inconsistency in returned array type (see code below) a bug or is it intentional? It is ...