Anasazi Version of the Day
Orthogonalization fault. More...
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.