A quadrature-based eigensolver with a Krylov subspace method for shifted linear systems for Hermitian eigenproblems in lattice QCD

DOI10.14495/JSIAML.2.115MaRDI QIDQ2843152zbMATH OpenOpenAlexFDO

Authors Hiroshi Ohno, Yoshinobu Kuramashi, Tetsuya Sakurai, Hiroto Tadano

Publication date 9 August 2013

Published in JSIAM Letters (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1004.0292, https://www.jstage.jst.go.jp/A_PRedirectJournalInit?sryCd=jsiaml&kijiCd=2_0_115&screenID=AF06S010&noVol=2&noIssue=0

zbMATH Keywords

numerical experimentseigenvalue problemHermitian matrixquadrature methodlattice QCDshifted linear systemsquantum chromodynamics (QCD)shifted conjugate gradient method

Mathematics Subject Classification ID

Hermitian, skew-Hermitian, and related matrices (15B57) Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Iterative numerical methods for linear systems (65F10) Strong interaction, including quantum chromodynamics (81V05)

Abstract: We consider a quadrature-based eigensolver to find eigenpairs of Hermitian matrices arising in lattice quantum chromodynamics. To reduce the computational cost for finding eigenpairs of such Hermitian matrices, we propose a new technique for solving shifted linear systems with complex shifts by means of the shifted CG method. Furthermore using integration paths along horizontal lines corresponding to the real axis of the complex plane, the number of iterations for the shifted CG method is also reduced. Some numerical experiments illustrate the accuracy and efficiency of the proposed method by comparison with a conventional method.

Recommendations

• Computing and deflating eigenvalues while solving multiple right-hand side linear systems with an application to quantum chromodynamics
• A deflated conjugate gradient method for multiple right hand sides and multiple shifts
• Numerical determination of partial spectrum of Hermitian matrices using a Lánczos method with selective reorthogonalization
• An accelerated conjugate gradient algorithm to compute low-lying eigenvalues -- a study for the Dirac operator in SU(2) lattice QCD
• BiCGStab() for families of shifted linear systems

Cited in

(8)

• Designing rational filter functions for solving eigenvalue problems by contour integration
• Kω -- open-source library for the shifted Krylov subspace method of the form \((z I - H)x = b\)
• Complex moment-based eigensolver coupled with two Krylov subspaces
• Computing partial spectra with least-squares rational filters
• Periodic band structure calculation by the Sakurai-Sugiura method with a fast direct solver for the boundary element method with the fast multipole representation
• Nonlinear eigenvalue problems and contour integrals
• Numerical determination of partial spectrum of Hermitian matrices using a Lánczos method with selective reorthogonalization
• Computing and deflating eigenvalues while solving multiple right-hand side linear systems with an application to quantum chromodynamics

Describes a project that uses

Uses Software

• CIRR

This page was built for publication: A quadrature-based eigensolver with a Krylov subspace method for shifted linear systems for Hermitian eigenproblems in lattice QCD

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2843152)