I was doing a physical problem, and then it comes to this Gaussian integral. The dimension of the integral is very large (dimension = 300~600), and it is difficult to find the maximum of the integrand. Is it possible to use Monte Carlo technique to do this integration?

To make the problem more specific, I write it in the following. I need to calculate C = A/B, where $B = \int d\bar{\phi_i}d\phi_i e^{-\phi_i^**\phi_i}f(\phi_i^*,\phi_i)$, $A = \int d\bar{\phi_i}d\phi_i e^{-\phi_i^**\phi_i}f(\phi_i^*,\phi_i)*H(\phi_i^*,\phi_i)$. The number of $\phi_i$, i.e., the dimension of the integration, can change from $10^2$ to $10^3$. The function $f(\phi_i^*,\phi_i)$ is positive definite, however it is difficult to find the maximum of $e^{-\phi_i^**\phi_i}f(\phi_i^*,\phi_i)$ and $f$ itself is very complex.

The problem is much like the original problem considered by Metropolis, who invented importance sampling with otheres. Is it really possible to calculate by Monte Carlo directly?