  We introduce a constrained energy functional to describe dielectric response.
We demonstrate that the local functional is a generalization of the long ranged
Marcus energy. Our re-formulation is used to implement a cluster Monte Carlo
algorithm for the simulation of dielectric media. The algorithm avoids solving
the Poisson equation and remains efficient in the presence of spatial
heterogeneity, nonlinearity and scale dependent dielectric properties.
