We have extended the
HiRep code [1]},The code is publicly available, and can be
accessed at https://github.com/claudiopica/HiRep
for the main \(SU(N_c)\)  version, and at https://github.com/sa2c/HiRep for the \(Sp(2N)\)  fork. to
adapt it to treat \(Sp(2N)\)  (rather than \(SU(N_c)\) )
gauge theories and couple them to fermions in multiple representations of the group.
Ensembles with dynamical fermions can be produced by combining the hybrid Monte Carlo (HMC) algorithm,
and its extension with rational approximations for the Dirac matrix with fractional powers—the rational
hybrid Monte Carlo (RHMC).
The standard (R)HMC algorithm consists of the following three main steps.
