NumRecipes Class Reference

#include <NumRecipes.h>

List of all members.

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)
Static Private Member Functions
static void	recipesErr (const char *c)
static void	recipesErr (const char *c)

---

Detailed Description

Definition at line 30 of file NumRecipes.h.

---

Member Function Documentation

static double NumRecipes::gammln

(

double

x

)

[static]

double NumRecipes::gammln

(

double

x

)

[static]

Definition at line 40 of file NumRecipes.cxx.

References ganga-rec::j, and x.

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

{
   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);
}

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

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

Definition at line 57 of file NumRecipes.cxx.

References gcf(), and gser().

{
   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;
   }
}

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

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

Definition at line 72 of file NumRecipes.cxx.

References gcf(), gser(), and recipesErr().

Referenced by ChisqConsistency::ChisqConsistency().

{
   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;
   }
}

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

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

Definition at line 86 of file NumRecipes.cxx.

References a0, exp(), gammln(), NUMREC_EPS, NUMREC_ITMAX, and recipesErr().

Referenced by gammp(), and gammq().

   {
     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");
   }

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

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

Definition at line 115 of file NumRecipes.cxx.

References exp(), gammln(), NUMREC_EPS, NUMREC_ITMAX, and recipesErr().

Referenced by gammp(), and gammq().

{
   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;
   }
}

static void NumRecipes::recipesErr

(

const char *

c

)

[static, private]

void NumRecipes::recipesErr

(

const char *

c

)

[static, private]

Definition at line 142 of file NumRecipes.cxx.

Referenced by gammq(), gcf(), and gser().

{
   std::cout<<"ErrMsg(fatal)" << " Numerical Recipes run-time error...\n" << c 
                 << "\n ...now exiting to system..." << std::endl;
   ::abort();
}

---