31 #define EPS_FD_JAC 1.0e-7 32 #define TOLF EPS_FD_JAC 33 #define TOLX EPS_FD_JAC 37 #define MAX_ROOT_EQNS 8 49 void LUBackSubst (
double **,
int,
int *,
double b[]);
52 void TridiagonalSolve(
double *,
double *,
double *,
double *,
double *,
int);
60 double BesselIO(
double);
62 double BesselKO(
double);
73 void ODE_Solve(
double *,
int,
double,
double,
double,
74 void (*rhs)(
double,
double *,
double *),
int method);
80 void QRUpdate(
double **,
double **,
int ,
double *,
double *);
81 void RSolve(
double **,
int ,
double *,
double *);
82 void QRSolve(
double **,
int ,
double *,
double *,
double *);
83 void QRDecompose (
double **,
int ,
double *,
double *,
int *);
84 void rotate(
double **,
double **,
int ,
int ,
double ,
double );
90 double GaussQuadrature(
double (*func)(
double,
void *),
void *,
double,
double,
97 #define PRNG_DEFAULT 0 102 #define PRNG PRNG_DEFAULT 109 unsigned int SeedGenerator(
long);
110 #if PRNG == PRNG_ECUYER 112 #elif PRNG == PRNG_MT 113 void init_genrand64(
unsigned long long);
114 unsigned long long genrand64_int64(
void);
115 double genrand64_real1(
void);
122 int CubicSolve (
double,
double,
double,
double *z);
123 int Brent(
double (*func)(
double,
void *),
void *,
double,
double,
124 double,
double,
double *);
125 void Broyden(
double *,
int ,
int *,
void (*vecfunc)(
int,
double *,
double *));
126 void FDJacobian(
int ,
double *,
double *,
double **,
127 void (*vecfunc)(
int ,
double *,
double *));
128 void LineSearch (
int ,
double *,
double ,
double *,
129 double *,
double *,
double *,
double *,
double ,
130 int *,
void ( *vecfunc)(
int ,
double *,
double *));
132 int QuarticSolve (
double,
double,
double,
double,
double *);
133 int Ridder(
double (*func)(
double,
void *),
void *,
134 double,
double,
double,
double,
double *);
141 void FinalizeTable2D (Table2D *);
151 double *a,
double *b,
double *c,
double *d);
154 void SplineCoeffs (
double *x,
double *f,
double dfL,
double dfR,
int n,
155 double *a,
double *b,
double *c,
double *d);
void MatrixMultiply(double **A, double **B, double **C, int n)
Definition: math_lu_decomp.c:163
double BesselK1(double x)
Definition: math_misc.c:159
void RandomSeed(long int seed, long int offset)
Definition: math_random.c:105
double BesselI1(double x)
Definition: math_misc.c:107
double BesselJ0(double x)
Definition: math_misc.c:13
void TridiagonalSolve(double *am, double *a0, double *ap, double *b, double *y, int n)
Definition: math_lu_decomp.c:182
double RandomNumber(double rmin, double rmax)
Definition: math_random.c:36
void MatrixInverse(double **A, double **Ainv, int n)
Definition: math_lu_decomp.c:137
double PowerLawRandomNumber(double xmin, double xmax, double n)
Definition: math_random.c:242
void SplineCoeffs(double *x, double *f, double dfL, double dfR, int n, double *a, double *b, double *c, double *d)
Definition: math_interp.c:95
void WriteBinaryTable2D(char *fname, Table2D *tab)
Definition: math_table2D.c:528
void LUBackSubst(double **a, int n, int *indx, double b[])
Definition: math_lu_decomp.c:92
int Brent(double(*Func)(double, void *), void *param, double x1, double x2, double abs_acc, double rel_acc, double *xroot)
Definition: math_root_finders.c:14
int InverseLookupTable2D(Table2D *tab, double y, double f, double *x)
Definition: math_table2D.c:196
void ODE_Solve(double *y0, int nvar, double xbeg, double xend, double dx, void(*rhs)(double, double *, double *), int method)
Definition: math_ode.c:25
void MonotoneSplineCoeffs(double *x, double *y, double *dydx, int n, double *a, double *b, double *c, double *d)
Definition: math_interp.c:13
double NR_ran2(long int *idum)
Definition: math_random.c:306
int Table2DInterpolate(Table2D *tab, double x, double y, double *f)
Definition: math_table2D.c:430
int CubicSolve(double b, double c, double d, double z[])
Definition: math_root_finders.c:591
int QuarticSolve(double b, double c, double d, double e, double *z)
Definition: math_root_finders.c:694
void QuickSort(int *x, int first, int last)
Definition: math_misc.c:210
double BesselJ1(double x)
Definition: math_misc.c:45
void VectorCartesianComponents(double *v, double x1, double x2, double x3)
Definition: math_misc.c:275
double GaussianRandomNumber(double mu, double sigma)
Definition: math_random.c:212
void SortArray(double *z, int n)
Definition: math_misc.c:249
void InitializeTable2D(Table2D *tab, double xmin, double xmax, int nx, double ymin, double ymax, int ny)
Definition: math_table2D.c:15
double BesselKn(int n, double x)
Definition: math_misc.c:184
int Ridder(double(*Func)(double, void *), void *param, double x1, double x2, double abs_acc, double rel_acc, double *xroot)
Definition: math_root_finders.c:120
double GaussQuadrature(double(*func)(double, void *), void *par, double xb, double xe, int nstep, int order)
Definition: math_quadrature.c:20
int LUDecompose(double **a, int n, int *indx, double *d)
Definition: math_lu_decomp.c:12
int QuadraticSolve(double a, double b, double c, double *x)
Definition: math_root_finders.c:545