This is the function GC which is required for the computation of the finite sample variance for m and total multivariance
sums.of.products(a, b, c, type = "multi")