- Added a test in solve to check whether the LU factorization or the
cholesky decomposition succeeded, otherwise a LinAlgError is raised
(similar to scipy's solve) + tests