1 #ifndef STAN_MATH_PRIM_SCAL_PROB_VON_MISES_RNG_HPP 2 #define STAN_MATH_PRIM_SCAL_PROB_VON_MISES_RNG_HPP 33 using boost::variate_generator;
38 static const char*
function(
"von_mises_rng");
43 double r = 1 +
pow((1 + 4 * kappa * kappa), 0.5);
44 double rho = 0.5 * (r -
pow(2 * r, 0.5)) / kappa;
45 double s = 0.5 * (1 + rho * rho) / rho;
51 W = (1 + s * Z) / (s + Z);
52 double Y = kappa * (s - W);
54 done = Y * (2 - Y) - U2 > 0;
57 done =
log(Y / U2) + 1 - Y >= 0;
61 double sign = ((U3 >= 0) - (U3 <= 0));
fvar< T > cos(const fvar< T > &x)
void check_finite(const char *function, const char *name, const T_y &y)
Check if y is finite.
fvar< T > log(const fvar< T > &x)
fvar< T > fmod(const fvar< T > &x1, const fvar< T > &x2)
void check_positive_finite(const char *function, const char *name, const T_y &y)
Check if y is positive and finite.
double von_mises_rng(double mu, double kappa, RNG &rng)
fvar< T > acos(const fvar< T > &x)
double pi()
Return the value of pi.
fvar< T > pow(const fvar< T > &x1, const fvar< T > &x2)
double uniform_rng(double alpha, double beta, RNG &rng)