Some recursions on Arnoldi's method and IOM for large non-Hermitian linear systems

✍ Scribed by Z. Jia

• Publisher: Elsevier Science
• Year: 2000
• Tongue: English
• Weight: 282 KB
• Volume: 39
• Category: Article
• ISSN: 0898-1221
• DOI: 10.1016/s0898-1221(99)00338-7

✦ Synopsis

Arnoldi's method and the Incomplete Orthogonalization Method (IOM) for large non-Hermitian linear systems are studied. It is shown that the inverse of a general nonsingular j x j Hessenberg matrix can be updated in O(j 2) flops from that of its (j -1) x (j -1) principal snbma~ trix. The updating recursion of inverses of the Hessenberg matrices does not need any QR or LU decompostion as commonly used in the literature. Some updating recursions of the residual norms and the approximate solutions obtained by these two methods are derived. These results are appealing because they allow one to decide when the methods converge and show one how to compute approximate solutions very cheaply and easily. (~) 2000 Elsevier Science Ltd. All rights reserved.