Tal-Ezer, Hillel; Turkel, Eli The iterative solver Risolv with application to the exterior Helmholtz problem. (English) Zbl 1209.65041 SIAM J. Sci. Comput. 32, No. 1, 463-475 (2010). Summary: The innermost computational kernel of many large-scale scientific applications is often a large set of linear equations of the form \(Ax=b\) which typically consumes a significant portion of the overall computational time required by the simulation. The traditional approach for solving this problem is to use direct methods. This approach is often preferred in industry because direct solvers are robust and effective for moderate size problems. However, direct methods can consume a huge amount of memory, and CPU time, in large-scale cases. In these cases, iterative techniques are the only viable alternative. Unfortunately, iterative methods lack the robustness of direct methods. The situation is especially difficult when the matrix is nonsymmetric. A lot of research has been devoted to trying to develop a robust iterative algorithm for nonsymmetric systems.The present paper describes a new robust and efficient algorithm aimed at solving iteratively nonsymmetric linear systems. It is based on looking for an approximation to the “optimal” polynomial \(P_m(z)\) which satisfies \(||P_m(z)||_{\infty}=\min_{Q\in\Pi_m}||Q(z)||_{\infty}, z\in D\), where \(\Pi_m\) is the set of all polynomials of degree \(m\) which satisfies \(Q_m(0)=1\) and \(D\) is a domain in the complex plane which includes all the eigenvalues of \(A\). The resulting algorithm is an efficient one, especially in the case where we have a set of linear systems which share the same matrix \(A\). We present several applications, including the exterior Helmholtz problem, which leads to a large indefinite, nonsymmetric, and complex system. Cited in 2 Documents MSC: 65F10 Iterative numerical methods for linear systems 65F08 Preconditioners for iterative methods 35J05 Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation 65N30 Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs 15A15 Determinants, permanents, traces, other special matrix functions 15A09 Theory of matrix inversion and generalized inverses 15A23 Factorization of matrices Keywords:Risolv algorithm; Krylov subspace method; numerical examples; GMRES method; finite element method; preconditioning; iterative methods; nonsymmetric linear systems; exterior Helmholtz problem Software:MatrixMarket; BiCGstab PDFBibTeX XMLCite \textit{H. Tal-Ezer} and \textit{E. Turkel}, SIAM J. Sci. Comput. 32, No. 1, 463--475 (2010; Zbl 1209.65041) Full Text: DOI