implement shift-and-invert mode for non-symmetric matrices
add niter information in result
use LU decomposition instead of matrix inverse
finish Cholesky decomposition for symmetric matrices
cannot use Cholesky decomposition, use LU instead
using LU decomposition for nonsymmetric matrices