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 np.finfo(x.dtype).tiny
in PoissonianEnergy.apply()
in order to address the problem, or should this be done somewhere else?