There's a bug in the Apache code to calculate Poisson masses.  For large integer realizations, 0 gets returned as the probability.   Not a good approximation when one is working on the log scale.  I hacked in my own logProbability calculator using Alexei's logGamma function.
There's a bug in the Apache code to calculate Poisson masses.  For large integer realizations, 0 gets returned as the probability.   Not a good approximation when one is working on the log scale.  I hacked in my own logProbability calculator using Alexei's logGamma function.

git-svn-id: 67bc77c75b8364e4e9cdff0eb6560f5818674cd8@628 ca793f91-a31e-0410-b540-2769d408b6a1

armanbilge/BEAST_sandbox,armanbilge/BEAST_sandbox,svn2github/beast-mcmc,svn2github/beast-mcmc,armanbilge/BEAST_sandbox,armanbilge/BEAST_sandbox,svn2github/beast-mcmc,svn2github/beast-mcmc,svn2github/beast-mcmc,armanbilge/BEAST_sandbox