uniformProposal.hpp

Go to the documentation of this file.

 
#ifndef MONTECARLO_1_UNIFORM_PROPOSAL_HPP
#define MONTECARLO_1_UNIFORM_PROPOSAL_HPP
 
#include "proposal.hpp"
#include "../domains/integration_domain.hpp"
 
#include <array>
#include <random>
#include <utility>
 
namespace mc {
namespace proposals {
 
template <size_t dim>
class UniformProposal : public Proposal<dim>
{
public:
    explicit UniformProposal(const mc::domains::IntegrationDomain<dim>& d);
 
    mc::geom::Point<dim> sample(std::mt19937& rng) const override;
 
    double pdf(const mc::geom::Point<dim>&) const override;
 
private:
    const mc::domains::IntegrationDomain<dim>& domain;
    mutable std::array<std::uniform_real_distribution<double>, dim> dist{};
    double vol_box;
};
 
} // namespace proposals
} // namespace mc
 
#include "uniformProposal.tpp"
 
#endif // MONTECARLO_1_UNIFORM_PROPOSAL_HPP