“…with on-shell internal lines, that admit a one-fold integral representations [59,61]. The entire decomposition of Feynman integrals in terms of a complete set of MIs, including the ones corresponding to subdiagrams, as well as the derivation of contiguity relations for special functions admitting multi-fold integral representation [61,63,65], required the use of multivariate intersection numbers [75,76,77,78,79,80,81,46,82]. A recursive algorithm for computing intersection numbers for meromorphic n-forms was proposed in [60] and later refined and applied to a few paradigmatic cases of Feynman integral decomposition [63,65,66].…”