Poisson distribution with rate `lambda` (> 0, ≤ 1e9). `function` selects the mass (pmf), cumulative distribution (cdf) or quantile / inverse-CDF.