Sampling from High-dimensional Distributions

Suppose we have a function f : R^p → R, the sampling problem is to output samples X ∼ π, where π(x) = (1/Z) exp(−f(x)), where Z is the normalization constant. In most cases, the normalization constant is computationally intractable, and so we would only have access to π up to proportionality.