quarks then one could apply standard algorithms such as [Mincer]{}, [@31; @32], to achieve high loop orders. However, from a technical point of view this renormalization is carried out at a point of exceptional momentum since the operator is inserted at zero momentum. In this configuration one is effectively dealing with a reduced $2$-point function. It has recently been pointed out, [@33], that for matrix elements relating to quark masses then one could have infrared issues in extracting reliable results in the non-perturbative region as a consequence of this momentum configuration. To circumvent this technicality an alternative scheme has been developed which is called the RI${}^\prime$/SMOM scheme, [@33]. It differs from the RI${}^\prime$ scheme in that $3$-point functions are not subtracted at an exceptional point but instead at a symmetric point. This means that none of the external momenta are nullified, so that the potential infrared singularity embedded within the logarithms of the Green’s function are bypassed, [@33]. Whilst the one loop computations were performed in [@33], this has recently been extended to two loops in [@34; @35] for the scalar and tensor currents. However, given that there is recent interest in measuring twist-$2$ flavour non-singlet operator matrix elements on 