SparseM_solve function

Helper function to slm_fit_csr