Switch to SVD instead of LUP

[rreusser]

I think a singular matrix means that there may just be a choice of schemes given by the nullspace of the matrix. I'd have to sit down and think about it, but if that's the case, then SVD would give a more complete solution.