Upon Fourier transformation the convolution becomes a product of the Fourier transform ${\cal F}[f]$ of the function $f$ and the Fourier transformed Gaussian measure, which is again a Gaussian with covariance matrix $C^{-1}$, $${\cal F}[\mu_{C}*f](k) = \exp\left(-\tfrac{1}{2}\sum_{n,m}k_n C_{nm} k_m\right){\cal F}[f](k).$$ Inverse Fourier transformation, with $k_n=i\partial/\partial x_n$, then gives Equation (1). This should hold irrespective of whether $f$ is polynomial or not.
Carlo Beenakker
- 188.2k
- 18
- 448
- 651