  Different recently developed Krylov space methods for solving linear systems
are studied and compared for the solution of the Dirac equation on the lattice.
Stabilized Biconjugate Gradient (BiCGstab2) is shown to be a robust and
efficient solver for the calculation of the Wilson quark propagator in lattice
QCD.
