[SciPy-dev] implement reorthogonalization method in gmres

[Darcoux Christine]

It seems that the vectors produced by the Gram-Schmidt process may not be
orthogonal (probably due to the presence of roundoff errors or when the
matrix is ill conditionned). I would suggest to add an option to test for
loss of orthogonality and reorthogonalize if needed. See
http://www4.ncsu.edu/~ctk/newton/SOLVERS/nsoli.m for an example.

regards,

Christine
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20100223/fb26a0a2/attachment.html>