As anticipated, the implementation of HMC/RHMC algorithms for fermions in arbitrary representations
of \(SU(N_c)\)  gauge groups is extensively discussed in Ref. [1]},
and its generalisation to the fundamental representation of \(Sp(2N)\)  in Ref. [2]}.
We pause here to discuss in further depth the case of multiple representations, given
the limited extent of the literature on the subject [3]}, [4]}, [5]}. In the rest of this subsection,
we follow closely the discussion in Ref. [1]}, and
refer the reader to this publication for details, while we highlight the differences
required in our implementation.
