1 #ifndef STAN_MCMC_HMC_NUTS_DIAG_E_NUTS_CLASSIC_HPP
2 #define STAN_MCMC_HMC_NUTS_DIAG_E_NUTS_CLASSIC_HPP
13 template <
class Model,
class BaseRNG>
16 expl_leapfrog, BaseRNG> {
26 Eigen::VectorXd& rho) {
27 return finish.
mInv.cwiseProduct(finish.
p).dot(rho - finish.
p) > 0
28 && finish.
mInv.cwiseProduct(start.
p).dot(rho - start.
p) > 0;
Probability, optimization and sampling library.
Point in a generic phase space.
bool compute_criterion(ps_point &start, diag_e_point &finish, Eigen::VectorXd &rho)
Point in a phase space with a base Euclidean manifold with diagonal metric.
diag_e_nuts_classic(const Model &model, BaseRNG &rng)