GMRES/CR and Arnoldi/Lanczos as Matrix Approximation Problems

نویسندگان

  • Anne Greenbaum
  • Lloyd N. Trefethen
چکیده

The GMRES and Arnoldi algorithms, which reduce to the CR and Lanczos algorithms in the symmetric case, both minimize p(A)b over polynomials p of degree n. The difference is that p is normalized at z 0 for GMRES and at z x for Arnoldi. Analogous "ideal GMRES" and "ideal Arnoldi" problems are obtained if one removes b from the discussion and minimizes p(/l)II instead. Investigation of these true and ideal approximation problems gives insight into how fast GMRES converges and how the Arnoldi iteration locates eigenvalues.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Tikhonov regularization via flexible Arnoldi reduction

Flexible GMRES, introduced by Saad, is a generalization of the standard GMRES method for the solution of large linear systems of equations. It is based on the flexible Arnoldi process for reducing a large square matrix to a small matrix. We describe how the flexible Arnoldi process can be applied to implement one-parameter and multi-parameter Tikhonov regularization of linear discrete ill-posed...

متن کامل

Investigating the Effects of MINRES with Local Reorthogonalization

The minimum residual method, MINRES, is an iterative method for solving a n x n linear system, Ax = b, where is A is a symmetric matrix. It searches for a vector, xk, in the k th Krylov subspace that minimizes the residual, rk = b − Axk. Due to the symmetric nature, the vectors are computed to be the Lanczos vectors, which simplifies to a three-term recurrence, where vk+1 is computed as a linea...

متن کامل

Krylov Subspace Approximation for TEM Simulation in the Time Domain

Forward transient electromagnetic modeling requires the numerical solution of a linear constant-coefficient initial-value problem for the quasi-static Maxwell equations. After discretization in space this problem reduces to a large system of ordinary differential equations, which is typically solved using finite-difference time-stepping. We compare standard time-stepping schemes such as the exp...

متن کامل

On Best Approximations of Polynomials in Matrices in the Matrix 2-Norm

We show that certain matrix approximation problems in the matrix 2-norm have uniquely defined solutions, despite the lack of strict convexity of the matrix 2-norm. The problems we consider are generalizations of the ideal Arnoldi and ideal GMRES approximation problems introduced by Greenbaum and Trefethen [SIAM J. Sci. Comput., 15 (1994), pp. 359–368]. We also discuss general characterizations ...

متن کامل

CSSP 13:733-758, 1994 Krylov Space Methods on State-Space Control Models

We give an overview of various Lanczos/Krylov space methods and how they are being used for solving certain problems in Control Systems Theory based on state-space models. The matrix methods used are based on Krylov sequences and are closely related to modern iterative methods for standard matrix problems such as sets of linear equations and eigenvalue calculations. We show how these methods ca...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • SIAM J. Scientific Computing

دوره 15  شماره 

صفحات  -

تاریخ انتشار 1994