Mathematics > Numerical Analysis
[Submitted on 27 Jan 2022 (v1), last revised 20 May 2022 (this version, v4)]
Title:GMRES using pseudoinverse for range symmetric singular systems
View PDFAbstract:Consider solving large sparse range symmetric singular linear systems $ A {\bf x}= {\bf b} $ which arise, for instance, in the discretization of convection diffusion equations with periodic boundary conditions, and partial differential equations for electromagnetic fields using the edge-based finite element method.
In theory, the Generalized Minimal Residual (GMRES) method converges to the least squares solution for inconsistent systems if the coefficient matrix $A$ is range symmetric, i.e. $ {\rm R}(A)= {\rm R}(A^{ \rm T } )$, where $ {\rm R}(A)$ is the range space of $A$.
We derived the necessary and sufficient conditions for GMRES to determine a least squares solution of inconsistent and consistent range symmetric systems assuming exact arithmetic except for the computation of the elements of the Hessenberg matrix.
In practice, GMRES may not converge due to numerical instability. In order to improve the convergence, we propose using the pseudoinverse for the solution of the severely ill-conditioned Hessenberg systems in GMRES. Numerical experiments on inconsistent systems indicate that the method is effective and robust. Finally, we further improve the convergence of the method by reorthogonalizing the Modified Gram-Schmidt procedure.
Submission history
From: Ken Hayami [view email][v1] Thu, 27 Jan 2022 10:32:42 UTC (1,846 KB)
[v2] Tue, 8 Feb 2022 00:54:40 UTC (2,340 KB)
[v3] Mon, 16 May 2022 09:30:36 UTC (2,340 KB)
[v4] Fri, 20 May 2022 13:05:19 UTC (2,340 KB)
Current browse context:
math.NA
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
alphaXiv (What is alphaXiv?)
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Hugging Face (What is Huggingface?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
Connected Papers (What is Connected Papers?)
CORE Recommender (What is CORE?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.