I implemented a Python algorithm using numpy.linalg.solve and wanted to optimize it with scipy.sparse.linalg.spsolve. To my surprise, results are numerically significantly different. I made a ...
Describe your issue. When one tries to solve a sparse linear system AX = B where A is of type scipy.sparse.csc.csc_matrix of size 46679 x 46679 and B is of type numpy.ndarray of size 46679 x 46680 ...