[SciPy-user] Solve for symmetric matrix?

[Anne Archibald]

On 26/04/07, Anand Patil <anand at soe.ucsc.edu> wrote:
> If the matrix has a zero eigenvalue, it has a nullspace so it's singular so you can't really invert it. However, if you know in advance that the RHS vector is orthogonal to the nullspace you can solve for the remaining components using a low-rank submatrix.

Let me reiterate my support for the SVD. It allows you to find a
least-squares solution for matrices with nullspaces, real or
numerical, even if the result vector is not in the span of the matrix.
It also lets you check how far your answer is from being a real
solution, in both spaces.



Anne