[Numpy-discussion] Solving Ax = b: inverse vs cholesky factorization
Mon Nov 8 14:17:46 CST 2010
On Mon, 08 Nov 2010 14:06:03 -0600, Bruce Southey wrote:
> Numpy uses SVD to get the (pseudo) inverse, which is usually very
> accurate at getting (pseudo) inverse.
solve(a, identity(a.shape, dtype=a.dtype))
It doesn't use xGETRI since that's not included in lapack_lite.
numpy.linalg.pinv OTOH does use SVD, but that's probably more costly.
More information about the NumPy-Discussion