  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.
