Should we "fix" exact zeros in PoissonianEnergy?
Several people seem to encounter problems when applying
PoissonianEnergy to a field that contains exact zeros.
Would it be acceptable to replace these zeros by
PoissonianEnergy.apply() in order to address the problem, or should this be done somewhere else?