error message if matrix is not square
no factorization if rhs is zero
default flag 'sum' is now 0 (--> unsymmetric matrix)