Как рассчитать число Пуассона в Scala

Question

Для этого нет стандартной библиотечной функции. Как я могу сделать это эффективно?

Answer 1

Я перевел this answer на Scala:

def recursive_poisson_helper(m:Long, r:Long, p:Double, i:Long):Double = { 
    if (r == i) { 
    p 
    } else { 
    recursive_poisson_helper(m, r, (p * m)/(i + 1), i + 1) 
    } 
} 

def efficient_poisson(m:Long, r:Long): Double = { 
    val p = math.exp(-m) 
    recursive_poisson_helper(m, r, p, 0) 
} 

Answer 2

Scala Breeze, http://www.scalanlp.org/ имеет класс Пуассона в его пакете stats.distributions.

case class Poisson(mean: Double)(implicit rand: RandBasis = Rand)