[Numpy-discussion] What approach is used in linalg.solve ?
Fri Sep 21 09:37:42 CDT 2007
If you take a look at the source of numpy's linalg.py, you'll see that
solves uses dgesv /zgesv for real /complex solves. If you Google dgesv, you
DGESV computes the solution to a real system of linear equations
A * X = B,
where A is an N-by-N matrix and X and B are N-by-NRHS matrices.
The LU decomposition with partial pivoting and row interchanges is
used to factor A as
A = P * L * U,
where P is a permutation matrix, L is unit lower triangular, and U is
upper triangular. The factored form of A is then used to solve the
system of equations A * X = B.
Don't take my word for it though; that was just the first google hit I
found. Also, I don't know if scipy solve differs from numpy.solve here, nor
which you are using, so I recommend that you repeat the exercise on your
On 9/21/07, mark <firstname.lastname@example.org> wrote:
> Hello, anybody know what approach is used in linalg.solve?
> I used it in a paper and some reviewer wants to know.
> Some Gaussian elimination with pivoting or something more fancy?
> Numpy-discussion mailing list
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Numpy-discussion