Wang-Landau algorithm for continuous models and joint density of states

We present a modified Wang-Landau algorithm for models with continuous degrees of freedom. We demonstrate this algorithm with the calculation of the joint density of states of ferromagnet Heisenberg models and a model polymer chain. The joint density of states contains more information than the density of states of a single variable-energy, but is also much more time consuming to calculate. We present strategies to significantly speed up this calculation for large systems over a large range of energy and order parameter.