NumRecipes Class Reference

#include <NumRecipes.h>

Static Public Member Functions
static double	gammln (double x)
static double	gammq (double a, double x)
static double	gammp (double a, double x)
static void	gcf (double *gammcf, double a, double x, double *gln)
static void	gser (double *gamser, double a, double x, double *gln)
static double	gammln (double x)
static double	gammq (double a, double x)
static double	gammp (double a, double x)
static void	gcf (double *gammcf, double a, double x, double *gln)
static void	gser (double *gamser, double a, double x, double *gln)

Detailed Description

Definition at line 30 of file MdcRecoUtil/MdcRecoUtil-00-01-08/MdcRecoUtil/NumRecipes.h.

Member Function Documentation

gammln() [1/2]

double NumRecipes::gammln ( double x )

static

Definition at line 40 of file NumRecipes.cxx.

{
   double x,y,tmp,ser;
   static double cof[6]={76.18009172947146,-86.50532032941677,
       24.01409824083091,-1.231739572450155,
       0.1208650973866179e-2,-0.5395239384953e-5};
   int j;
 
   y=x=xx;
   tmp=x+5.5;
   tmp -= (x+0.5)*log(tmp);
   ser=1.000000000190015;
   for (j=0;j<=5;j++) ser += cof[j]/++y;
   return -tmp+log(2.5066282746310005*ser/x);
}

Referenced by ChisqConsistency::ChisqConsistency(), ChisqConsistency::ChisqConsistency(), gcf(), and gser().

gammln() [2/2]

static double NumRecipes::gammln ( double x )

static

gammp() [1/2]

double NumRecipes::gammp ( double a, double x )

static

Definition at line 57 of file NumRecipes.cxx.

{
   double gamser,gammcf,gln;
 
   if (x < 0.0 || a <= 0.0) std::cout<<"ErrMsg(error)" <<" Invalid arguments in routine gammp x=" << x << " a=" << a << endl;
   if (x < (a+1.0)) {
      gser(&gamser,a,x,&gln);
      return gamser;
   } else {
      gcf(&gammcf,a,x,&gln);
      return 1.0-gammcf;
   }
}

gammp() [2/2]

static double NumRecipes::gammp ( double a, double x )

static

gammq() [1/2]

double NumRecipes::gammq ( double a, double x )

static

Definition at line 72 of file NumRecipes.cxx.

{
   double gamser,gammcf,gln;
 
   if (x < 0.0 || a <= 0.0) recipesErr(" Invalid arguments in routine GAMMQ");
   if (x < (a+1.0)) {
      gser(&gamser,a,x,&gln);
      return 1.0-gamser;
   } else {
      gcf(&gammcf,a,x,&gln);
      return gammcf;
   }
}

Referenced by ChisqConsistency::ChisqConsistency().

gammq() [2/2]

static double NumRecipes::gammq ( double a, double x )

static

gcf() [1/2]

void NumRecipes::gcf ( double * gammcf, double a, double x, double * gln )

static

Definition at line 86 of file NumRecipes.cxx.

   {
     int n;
     double gold=0.0,g,fac=1.0,b1=1.0;
     double b0=0.0,anf,ana,an,a1,a0=1.0;
 
     *gln=gammln(a);
     a1=x;
     for (n=1;n<=NUMREC_ITMAX;n++) {
       an=(double) n;
       ana=an-a;
       a0=(a1+a0*ana)*fac;
       b0=(b1+b0*ana)*fac;
       anf=an*fac;
       a1=x*a0+anf*a1;
       b1=x*b0+anf*b1;
       if (a1) {
         fac=1.0/a1;
         g=b1*fac;
         if (fabs((g-gold)/g) < NUMREC_EPS) {
            *gammcf=exp(-x+a*log(x)-(*gln))*g;
            return;
         }
         gold=g;
       }
     }
     recipesErr(" a too large, NUMREC_ITMAX too small in routine GCF");
   }

Referenced by gammp(), and gammq().

gcf() [2/2]

static void NumRecipes::gcf ( double * gammcf, double a, double x, double * gln )

static

gser() [1/2]

void NumRecipes::gser ( double * gamser, double a, double x, double * gln )

static

Definition at line 115 of file NumRecipes.cxx.

{
   int n;
   double sum,del,ap;
 
   *gln=gammln(a);
   if (x <= 0.0) {
      if (x < 0.0) recipesErr(" x less than 0 in routine GSER");
      *gamser=0.0;
      return;
   } else {
      ap=a;
      del=sum=1.0/a;
      for (n=1;n<=NUMREC_ITMAX;n++) {
         ap += 1.0;
         del *= x/ap;
         sum += del;
         if (fabs(del) < fabs(sum)*NUMREC_EPS) {
            *gamser=sum*exp(-x+a*log(x)-(*gln));
            return;
         }
      }
      recipesErr(" a too large, NUMREC_ITMAX too small in routine GSER");
      return;
   }
}

Referenced by gammp(), and gammq().

gser() [2/2]

static void NumRecipes::gser ( double * gamser, double a, double x, double * gln )

static

---

The documentation for this class was generated from the following files:

• MdcRecoUtil/MdcRecoUtil-00-01-08/MdcRecoUtil/NumRecipes.h
• ProbTools/ProbTools-00-00-01/ProbTools/NumRecipes.h
• NumRecipes.cxx