if (p <= pow(1 - dist.success_fraction(), dist.trials())) { p <= pdf(dist, 0) == cdf(dist, 0) return 0; So the only reasonable result is zero.
if (p <= pow(1 - dist.success_fraction(), dist.trials())) { p <= pdf(dist, 0) == cdf(dist, 0)
return 0; So the only reasonable result is zero.
} And root finder would fail otherwise. Added, and similarly for complement