Scipy sparse invert or spsolve lead to UMFPACK_ERROR_OUT_OF_MEMORY

Question

I am trying to invert a large (150000,150000) sparse matrix as follows:

import scipy as sp
import scipy.sparse.linalg as splu

#Bs is a large sparse         


        

Answer 1

A sparse array only fits the non-zero entries of your matrix into memory. Now suppose you do an inversion. This means that almost all entries of the matrix become non-zero. Sparse matrices are memory optimized.

There are some operations which you can apply on sparse matrices without losing the "spare" property:

• Addition, just adding a constant can keep the sparse matrix sparse.