I am trying to obtain the Cholesky decomposition of a huge $150,000 \times 150,000$ sparse matrix with randomly distributed non-zero elements. I have only the entries for which the values are non-zero.
I tried using the sparse Cholesky decomposition from SciPy library but it gave me an error of not enough memory. Any suggestions for efficient ways to proceed with the Cholesky decomposition?