  A scheme for separating the high- and low-frequency molecular dynamics modes
in Hybrid Monte Carlo (HMC) simulations of gauge theories with dynamical
fermions is presented. The algorithm is tested in the Schwinger model with
Wilson fermions.
