22 double v00,
double v01,
double v02,
double v10,
double v11,
double v20);
31 double invertFunc(
double (*func)(
double),
double y,
double xL,
double xR,
double dx,
double x0,
double relative_accuracy=1
e-10);
34 double invertTableDirect(vector<double>*
x, vector<double>*
y,
double y0,
double x0,
double relative_accuracy=1
e-10);
44 double adaptiveSimpsonsAux(
double (*
f)(
double),
double a,
double b,
double epsilon,
double S,
double fa,
double fb,
double fc,
int bottom);
45 double adaptiveSimpsons(
double (*
f)(
double),
double a,
double b,
double epsilon=1
e-15,
int maxRecursionDepth=50);
47 double qiu_simpsons(
double (*
f)(
double),
double a,
double b,
double epsilon=1
e-15,
int maxRecursionDepth=50);
48 double qiu_simpsonsRel(
double (*
f)(
double),
double a,
double b,
double epsilon=1
e-15,
int maxRecursionDepth=50);
60 inline long irand(
long LB,
long RB)
63 return LB + (long)((RB-LB+1)*(drand48()-1
e-25));
66 inline double drand(
double LB,
double RB)
71 double dw = width*1
e-30;
72 return LB+dw+(width-2*dw)*drand48();
84 void get_bin_average_and_count(istream& is, ostream&
os, vector<double>*
bins,
long col_to_bin=0,
void (*func)(vector<double>*)=NULL,
long wanted_data_columns=-1,
bool silence=
false);