I have a sparse matrix, symmetric, real. Dimension 101680 , number of non zeros ~6M.
Calling cholesky
takes around 11 seconds.
@time K = SuiteSparse.CHOLMOD.cholesky(K)
Calling ldlt
results in several minutes of “nothing happening”, by the end of which there is still no result.
@time K = SuiteSparse.CHOLMOD.ldlt(K)
I wonder if anyone has been successful using this solver?
Edit: ldlt took almost half an hour.