[SciPy-dev] implement reorthogonalization method in gmres

Darcoux Christine bouloumag@gmail....
Tue Feb 23 19:26:02 CST 2010


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.scipy.org/pipermail/scipy-dev/attachments/20100223/fb26a0a2/attachment.html 


More information about the SciPy-Dev mailing list