We hence have four contributions to the MD evolution:
the gauge contribution is supplemented by those coming from the
HMC treatment of
the \((f)\)  pseudofermion, from the HMC treatment of one \((as)\)  pseudofermion,
and from the RHMC application to the third, \((as)\)  pseudofermion.
We illustrate the size of each, by showing in Fig. REF 
their contribution to the force as it enters the Hamiltonian evolution—see the
Hamilton equations governing the MD evolution, in Eqs. (19) and (20) of Ref. [1]}—averaged
over one of the ensembles.
The acceptance rate is in the range of \(75-85\%\) . To accelerate the
(computationally demanding) inversion of the
Dirac operator \(Q^2_m\) ,
we use the second order Omelyan integrator [2]}
in the MD evolution
and the even-odd preconditioning of the fermion matrix [3]},
applied to the (R)HMC algorithm as discussed in Ref. [1]}.
