Returns condition number of zero if matrix is not invertible.

It seems the rcond estimate algorithm implicitly assumes the matrix is invertible when it does solves of various e vectors. As is, it can fail to converge. After LU decomposition, checking the rank is only O(n), so this does not significantly impact overall computational complexity.

Fixes #2889 (closed).