Fractional Inversion in Krylov Space

Abstract

The fractional inverse $M^{-\gamma }$ (real $\gamma >0$) of a matrix $M$ is expanded in a series of Gegenbauer polynomials. If the spectrum of $M$ is confined to an ellipse not including the origin, convergence is exponential, with the same rate as for Chebyshev inversion. The approximants can be improved recursively and lead to an iterative solver for $M^{\gamma }x=b$ in Krylov space. In case of $\gamma =1/2$, the expansion is in terms of Legendre polynomials, and rigorous bounds for the truncation error are derived.

Contribution to LAT97 proceedings, 3 pages

Cited by in corpus (15)

• Exact chiral symmetry, topological charge and related topics
• Numerical techniques for lattice QCD in the $\epsilon$--regime
• Topological Charge and The Spectrum of Exactly Massless Fermions on the Lattice
• On the Neuberger overlap operator
• Locality with staggered fermions
• Numerical Methods for the QCD Overlap Operator:III. Nested Iterations
• Filtered overlap: speedup, locality, kernel non-normality and Z_A~1
• Topology and chiral symmetry breaking in SU(N) gauge theories
• Light hadron and diquark spectroscopy in quenched QCD with overlap quarks on a large lattice
• Lanczos Approach to the Inverse Square Root of a Large and Sparse Matrix
• A comparative study of numerical methods for the overlap Dirac operator--a status report
• An Exact Algorithm for Any-flavor Lattice QCD with Kogut-Susskind Fermion
• Improving the Dirac Operator in Lattice QCD
• Chirally improved Dirac operators: Studying the sensitivity to topological excitations for zero and finite temperature
• Recent results using the overlap Dirac operator