I need to numerically evaluate/approximate high-order moments of high-dimensional Gaussian measures/distributions with given mathematical expectations and covariance matrices. The Gaussian measures have dimension $ d$ with say $ d>1000$ and the moments have order $ k$ with say $ k>>1000$ . Hence the classical Wick-Isserlis theorem/formula looks useless because it would involve an astronomically largeRead more