|
Anasazi Version of the Day
|
Orthogonalization fault. More...
#include <AnasaziTsqrOrthoManager.hpp>
Orthogonalization fault.
Stewart (SISC 2008) gives a Block Gram-Schmidt (BGS) with reorthogonalization algorithm. An "orthogonalization fault" is what happens when the second BGS pass does not succeed (which is possible in BGS, but not possible in (non-block) Gram-Schmidt if you use Stewart's randomization procedure). Stewart gives an algorithm for recovering from an orthogonalization fault, but the algorithm is expensive: it involves careful reorthogonalization with non-block Gram-Schmidt. We choose instead to report the orthogonalization fault and let users recover from it.
Definition at line 72 of file AnasaziTsqrOrthoManager.hpp.
1.7.4