PoissonDistribution.h

Go to the documentation of this file.

#ifndef _utl_PoissonDistribution_h_
#define _utl_PoissonDistribution_h_


#include <utl/Math.h>


namespace utl {

  typedef unsigned long ulong;

  inline double LogPoissonPDF(const ulong k, const double x)
  { return -x - LogGamma(k+1) + k * log(x); }

  inline double PoissonPDF(const ulong k, const double x)
  { return exp(LogPoissonPDF(k,x)); }

}


#endif