  The short-range modes of the fermionic determinant can be absorbed in the
gauge action using the loop expansion. The coefficients of this expansion and
the zeroes of the polynomial approximating the remainder can be optimized by a
simple, practical method. When the multiboson approach is used, this
optimization results in a faster simulation with fewer auxiliary fields.
