[SciPy-user] null space of a matrix

[Nils Wagner]

On Wed, 29 Jun 2005 12:15:24 -0700
  Robert Kern <rkern at ucsd.edu> wrote:
> 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.
> 

Is it a great effort to add this as a built-in function
e.g. null(A) similar to Matlab's function ?

Nils

> -- 
> 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
> 
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-user