[Numpy-discussion] linear algebra help

[Sebastian Walter]

2009/5/18 Stéfan van der Walt <stefan at sun.ac.za>:
> 2009/5/18 Sebastian Walter <sebastian.walter at gmail.com>:
>> B = numpy.dot(A.T, A)
>
> This multiplication should be avoided whenever possible -- you are
> effectively squaring your condition number.
Indeed.

>
> In the case where you have more rows than columns, use least squares.
> For square matrices use solve.  For large sparse matrices, use GMRES
> or any of the others available in scipy.sparse.linalg.

It is my impression that this is a linear algebra and not a numerics question.



>
> Regards
> Stéfan
> _______________________________________________
> Numpy-discussion mailing list
> Numpy-discussion at scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>