A simplified and flexible variant of GCROT for solving nonsymmetric linear systems

JE Hicken, DW Zingg - SIAM Journal on Scientific Computing, 2010 - SIAM
SIAM Journal on Scientific Computing, 2010SIAM
There is a need for flexible iterative solvers that can solve large-scale (>10^6 unknowns)
nonsymmetric sparse linear systems to a small tolerance. Among flexible solvers, flexible
GMRES (FGMRES) is attractive because it minimizes the residual norm over a particular
subspace. In practice, FGMRES is often restarted periodically to keep memory and work
requirements reasonable; however, like restarted GMRES, restarted FGMRES can suffer
from stagnation. This has led us to develop a flexible variant of the Krylov linear solver …
There is a need for flexible iterative solvers that can solve large-scale ( unknowns) nonsymmetric sparse linear systems to a small tolerance. Among flexible solvers, flexible GMRES (FGMRES) is attractive because it minimizes the residual norm over a particular subspace. In practice, FGMRES is often restarted periodically to keep memory and work requirements reasonable; however, like restarted GMRES, restarted FGMRES can suffer from stagnation. This has led us to develop a flexible variant of the Krylov linear solver GCROT (generalized conjugate residual with inner orthogonalization and outer truncation). Unlike the original GCROT algorithm, the proposed GCROT variant uses a simplified truncation strategy similar to loose GMRES (LGMRES). This modification is motivated by numerical experiments that suggest the specific subspace retained in the outer iteration of GCROT is less important than its size. The flexible GCROT variant appears to be well suited for advection-dominated problems. In particular, when applied to an adjoint problem from computational aerodynamics, the proposed GCROT variant is robust and efficient compared with several popular truncated Krylov subspace methods. Finally, a flexible version of LGMRES is easily constructed by recognizing algorithmic similarities to GCROT.
Society for Industrial and Applied Mathematics
この検索の最上位の結果を表示しています。 検索結果をすべて見る