Integrating products of quadratic forms

We prove that if $q1, \ldots, qm: {\Bbb R}ⁿ \longrightarrow {\Bbb R}$ are quadratic forms in variables $x1, \ldots, xn$ such that each $qk$ depends on at most $r$ variables and each $qk$ has common variables with at most $r$ other forms, then the average value of the product $\left(1+ q1\right) \cdots \left(1+qm\right)$ with respect to the standard Gaussian measure in ${\Bbb R}^n$ can be approximated within relative error $\epsilon >0$ in quasi-polynomial $n^{O(1)} m^{O(\ln m -\ln \epsilon)}$ time, provided $|qk(x)| \leq \gamma |x|² /r$ for some absolute constant $\gamma > 0$ and $k=1, \ldots, m$. When $qk$ are interpreted as pairwise squared distances for configurations of points in Euclidean space, the average can be interpreted as the partition function of systems of particles with mollified logarithmic potentials. We sketch a possible application to testing the feasibility of systems of real quadratic equations.