2 GENERATE POISSON RANDOM NUMBER IF alambda IS LESS THAN xover
3 OR ELSE GENERATE GAUSSIAN RANDOM NUMBER
10 #include <DIGRand/DIGRand.h>
11 //#include "DIGGauss.h"
12 #include <DIGRand/DIGGauss.h>
13 //#include "DIGPoisson.h"
14 #include <DIGRand/DIGPoisson.h>
17 int Poisson(double pALambda)
19 static const double xover = 36.0;
26 if(pALambda >= xover) {
27 int g = (int) (Gauss(pALambda, (double) sqrt(pALambda)) + 0.5);
33 e = (double) exp(-pALambda);