void set_nominal_stepsize_and_T(const double e, const double t) {
if (e > 0 && t > 0) {
this->nom_epsilon_ = e;
T_ = t;
update_L_();
}
}
void set_nominal_stepsize(const double e) {
if (e > 0) {
this->nom_epsilon_ = e;
update_L_();
}
}
void set_T(const double t) {
if (t > 0) {
T_ = t;
update_L_();
}
}
base_static_hmc(M &m, BaseRNG& rng, std::ostream* o, std::ostream* e)
: base_hmc<M, P, H, I, BaseRNG>(m, rng, o, e), T_(1) {
update_L_();
}
base_static_hmc(const Model& model, BaseRNG& rng)
: base_hmc<Model, Hamiltonian, Integrator, BaseRNG>(model, rng),
T_(1), energy_(0) {
update_L_();
}
base_static_hmc(Model &model, BaseRNG& rng,
std::ostream* o, std::ostream* e)
: base_hmc<Model, Hamiltonian, Integrator, BaseRNG>(model, rng, o, e),
T_(1) {
update_L_();
}
