Another random thought: Do full tetrahedron integration in BCC space where metric is a general symmetric tensor. This leads 2nd moments of standard tetrahedron which can’t be too horrible. This may work because we use geometry to separate the variables.

This still turns out to take the integral of the −½ power of a linear combination of moments. Same nearly intractable problem as before.