[Numpy-discussion] np.dot and scipy sparse matrices

[Peter Prettenhofer]

Hi everybody,

I recently got the latest numpy version (2.0.0.dev-7297785) from the
git repo and realized that `np.dot` causes a segfault if its operands
are scipy sparse matrices. Here's some code to reproduce the problem::

   import numpy as np
   from scipy import sparse as sp
   A = np.random.rand(10, 10)
   S = sp.csr_matrix(A)
   _ = np.dot(A, A)   # this works OK
   _ = np.dot(S, S)   # this segfaults!

thanks,
 Peter
-- 
Peter Prettenhofer