Go to the documentation of this file.
66 static double hyp2f0(
double a,
double b,
double x,
int type,
double *err)
68 double a0, alast,
t, tlast, maxt;
69 double n, an, bn, u, sum, temp;
87 u = an * (bn *
x /
n);
91 if ((temp > 1.0) && (maxt > (DBL_MAX / temp)))
139 (0.5 + (0.125 + 0.25 *
b - 0.5 *
a + 0.25 *
x - 0.25 *
n) /
x);
143 alast *= 2.0 / 3.0 -
b + 2.0 *
a +
x -
n;
181 static double hy1f1a(
double a,
double b,
double x,
double *err)
183 double h1, h2,
t, u, temp, acanc,
asum, err1, err2;
191 t =
x + temp * (
a -
b);
200 h1 =
hyp2f0(
a,
a -
b + 1, -1.0 /
x, 1, &err1);
206 h2 =
hyp2f0(
b -
a, 1.0 -
a, 1.0 /
x, 2, &err2);
237 if (
asum == INFINITY ||
asum == -INFINITY)
251 static double hy1f1p(
double a,
double b,
double x,
double *err)
253 double n, a0, sum,
t, u, temp, maxn;
283 u =
x * (an / (bn *
n));
287 if ((temp > 1.0) && (maxt > (DBL_MAX / temp))) {
296 c = (sumc - sum) -
y;
310 *err =
fabs(
c / sum);
328 double asum, psum, acanc, pcanc, temp;
Array< double, 1, 3 > e(1./3., 0.5, 2.)
RealScalar EIGEN_BLAS_FUNC() asum(int *n, RealScalar *px, int *incx)
set noclip points set clip one set noclip two set bar set border lt lw set xdata set ydata set zdata set x2data set y2data set boxwidth set dummy x
const EIGEN_DEVICE_FUNC LogReturnType log() const
const EIGEN_DEVICE_FUNC ExpReturnType exp() const
static double hy1f1a(double a, double b, double x, double *err)
static double hy1f1p(double a, double b, double x, double *err)
static double hyp2f0(double a, double b, double x, int type, double *err)
void sf_error(const char *func_name, sf_error_t code, const char *fmt,...)
double hyperg(double a, double b, double x)
gtsam
Author(s):
autogenerated on Sat Nov 16 2024 04:02:26