an:00903501
Zbl 0856.65043
??tuller, J??lius
Ordered modified Gram-Schmidt orthogonalization revised
EN
J. Comput. Appl. Math. 63, No. 1-3, 221-227 (1995).
00034805
1995
j
65F25
ordered modified Gram-Schmidt orthogonalization; Gram-Schmidt algorithm
\textit{B. Nour-Omid}, \textit{W. S. Dunbar} and \textit{A. D. Woodbury} [Ordered modified Gram-Schmidt orthogonalization, unpublished paper (1991)] proposed an ordered modified Gram-Schmidt algorithm, which was supposed to improve the orthogonality state of the solution. Thorough analysis, even of a simple, planar case, shows [cf. \textit{J. ??tuller}, Some comments on the ordered modified Gram-Schmidt orthogonalization, Research Report No. 586, Institute of Computer Science, Academy of Sciences of the Czech Republic, Prague, (1994)] that, yet in the exact arithmetic, one cannot expect to obtain -- independently of the forward, reverse, or any other type, including ``ordered'', orthogonalization -- the desired solution: a vector in the orthogonal complement of the given vectors. For the planar case, some simple rules of the thumb can be given to minimize the error, but we are not sure whether they can directly be generalized to higher dimensions as has been done by Nour-Omid et al. [loc. cit.]. Naturally, in finite precision, the situation is even worse. The unique way out we see in the iterative Gram-Schmidt methods [cf. \textit{A. Bj??rck}, Linear Algebra Appl. 197-198, 297-316 (1994; Zbl 0801.65039)] for which we [loc. cit.] presented several algorithms with different criterions of efficiency applied to.
D.Braess (Bochum)
Zbl 0801.65039