[SciPy-user] null space of a matrix

[Robert Kern]

Ryan Krauss wrote:
> This is probably a Numeric question, but is there an easy way to find 
> the null space of a matrix in SciPy?

The SVD is one way.

u, s, vh = linalg.svd(A)
null_mask = (s <= eps)
null_space = compress(null_mask, vh, axis=0)

The rows in null_space span the null space of A.

-- 
Robert Kern
rkern at ucsd.edu

"In the fields of hell where the grass grows high
  Are the graves of dreams allowed to die."
   -- Richard Harter