  Quark propagators with arbitrary sources and sinks can be obtained more
efficiently using a pseudofermion method with a mode-shifted action.
Mode-shifting solves the problem of critical slowing down (for light quarks)
induced by low eigenmodes of the Dirac operator. The method allows the full
physical content of every gauge configuration to be extracted, and should be
especially helpful for unquenched QCD calculations. The method can be applied
for all the conventional quark actions: Wilson, Sheikoleslami-Wohlert,
Kogut-Susskind, as well as Ginsparg-Wilson compliant overlap actions. The
statistical properties of the method are examined and examples of physical
processes under study are presented.
