GP-CMRH: An inner product free iterative method for block two-by-two nonsymmetric linear systems

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Du, Kui, Fan, Jia-Jun
Natura: Preprint
Pubblicazione: 2025
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866916949124972544
author Du, Kui
Fan, Jia-Jun
author_facet Du, Kui
Fan, Jia-Jun
contents We propose an inner product free iterative method called GP-CMRH for solving block two-by-two nonsymmetric linear systems. GP-CMRH relies on a new simultaneous Hessenberg process that reduces two rectangular matrices to upper Hessenberg form simultaneously, without employing inner products. Compared with GPMR [SIAM J. Matrix Anal. Appl., 44 (2023), pp. 293--311], GP-CMRH requires less computational cost per iteration and may be more suitable for high performance computing and low or mixed precision arithmetic due to its inner product free property. Our numerical experiments demonstrate that GP-CMRH and GPMR exhibit comparable convergence behavior (with GP-CMRH requiring slightly more iterations), yet GP-CMRH consumes less computational time in most cases. GP-CMRH significantly outperforms GMRES and CMRH in terms of convergence rate and runtime efficiency.
format Preprint
id arxiv_https___arxiv_org_abs_2509_11272
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle GP-CMRH: An inner product free iterative method for block two-by-two nonsymmetric linear systems
Du, Kui
Fan, Jia-Jun
Numerical Analysis
15A06, 65F10, 65F50
We propose an inner product free iterative method called GP-CMRH for solving block two-by-two nonsymmetric linear systems. GP-CMRH relies on a new simultaneous Hessenberg process that reduces two rectangular matrices to upper Hessenberg form simultaneously, without employing inner products. Compared with GPMR [SIAM J. Matrix Anal. Appl., 44 (2023), pp. 293--311], GP-CMRH requires less computational cost per iteration and may be more suitable for high performance computing and low or mixed precision arithmetic due to its inner product free property. Our numerical experiments demonstrate that GP-CMRH and GPMR exhibit comparable convergence behavior (with GP-CMRH requiring slightly more iterations), yet GP-CMRH consumes less computational time in most cases. GP-CMRH significantly outperforms GMRES and CMRH in terms of convergence rate and runtime efficiency.
title GP-CMRH: An inner product free iterative method for block two-by-two nonsymmetric linear systems
topic Numerical Analysis
15A06, 65F10, 65F50
url https://arxiv.org/abs/2509.11272