Accelerating Wilson Fermion Matrix Inversions by Means of the Stabilized   Biconjugate Gradient Algorithm
  The stabilized biconjugate gradient algorithm BiCGStab recently presented by van der Vorst is applied to the inversion of the lattice fermion operator in the Wilson formulation of lattice Quantum Chromodynamics. Its computational efficiency is tested in a comparative study against the conjugate gradient and minimal residual methods. Both for quenched gauge configurations at beta= 6.0 and gauge configurations with dynamical fermions at beta=5.4, we find BiCGStab to be superior to the other methods. BiCGStab turns out to be particularly useful in the chiral regime of small quark masses. 
A. Frommer, V. Hannemann, Th. Lippert, B. Noeckel and K. Schilling