1 #ifndef STAN_MATH_FWD_SCAL_FUN_GAMMA_Q_HPP
2 #define STAN_MATH_FWD_SCAL_FUN_GAMMA_Q_HPP
33 while (
fabs(delta) > 1
e-6) {
37 delta = s / ((k + x1.
val_) * (k + x1.
val_));
40 T der1 = (1.0 - u) * (dig - l) +
exp(x1.
val_ * l) * S / g;
68 while (
fabs(delta) > 1
e-6) {
72 delta = s / ((k + x1.
val_) * (k + x1.
val_));
75 T der1 = (1.0 - u) * (dig - l) +
exp(x1.
val_ * l) * S / g;
fvar< T > fabs(const fvar< T > &x)
fvar< T > log(const fvar< T > &x)
fvar< T > exp(const fvar< T > &x)
double e()
Return the base of the natural logarithm.
fvar< T > pow(const fvar< T > &x1, const fvar< T > &x2)
fvar< T > tgamma(const fvar< T > &x)
fvar< T > gamma_q(const fvar< T > &x1, const fvar< T > &x2)
fvar< T > digamma(const fvar< T > &x)