Block Algorithms for Quark Propagator Calculation

Abstract

Computing quark propagators in lattice QCD is equivalent to solving large, sparse linear systems with multiple right-hand sides. Block algorithms attempt to accelerate the convergence of iterative Krylov-subspace methods by solving the multiple systems simultaneously. This paper compares a block generalisation of the quasi-minimal residual method (QMR), Block Conjugate Gradient on the normal equation, Block Lanczos and (${\hat{I}}_5^3$-symmetric) Block BiConjugate Gradient.

3 pages, 1 figure, LaTeX2e, uses espcrc2 and epsf. Poster presented at Lattice '97