Arnoldi’s method and the incomplete orthogonalization method (IOM) for large non-Hermitian linear systerns are studied. It is shown that the inverse of a general nonsingular j × j Hessenberg matrir can be update...Arnoldi’s method and the incomplete orthogonalization method (IOM) for large non-Hermitian linear systerns are studied. It is shown that the inverse of a general nonsingular j × j Hessenberg matrir can be updated in O(j2) flops from that of its (j -1) × (j - 1) principal submatrir. 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.展开更多
文摘Arnoldi’s method and the incomplete orthogonalization method (IOM) for large non-Hermitian linear systerns are studied. It is shown that the inverse of a general nonsingular j × j Hessenberg matrir can be updated in O(j2) flops from that of its (j -1) × (j - 1) principal submatrir. 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.
