Algorithms 3 and 4 depict how to compute the long and short permutation contributions to R abc ijk , respectively. The computation of the R abc ijk residual has been adapted from the work of Lesiuk et al. . Using the conjugate permutation operator, from the permutation operator Algorithm 2 Linear-scaling computation of triples contribution to R ab ij (Equation )