sampler_tools.hpp

Go to the documentation of this file.

#pragma once
#include <chrono>
 
namespace ppl {
namespace mcmc {
 
inline size_t random_seed()
{
    return std::chrono::duration_cast<std::chrono::milliseconds>
        (std::chrono::system_clock::now().time_since_epoch()).count();
}
 
template <class UniformDistType, class GenType>
inline bool accept_or_reject(double p, 
                             UniformDistType&& unif_sampler,
                             GenType&& gen)
{
    return (unif_sampler(gen) <= p);
}
 
} // namespace mcmc
} // namespace ppl