Add check for not positive definite (in factorize)