U_sampler.hpp

Go to the documentation of this file.

 
#pragma once
 
#include "../../utils/Data.hpp"
#include "../../utils/Params.hpp"
#include <random>

class U_sampler {
 
protected:
  Params &params;

  Data &data;

  int total_iterations = 0;

  int accepted_U = 0;

  std::random_device rd;

  mutable std::mt19937 gen;


  const double a_over_sigma = params.a / params.sigma;

  const double tau_power_sigma = std::pow(params.tau, params.sigma);

  const int n = data.get_n();

  double U = 1.0;



  double log_conditional_density_U(double u) const;

  double log_conditional_density_V(double v) const;

 
public:
  U_sampler(Params &p, Data &d) : params(p), data(d) {};

  virtual void update_U() = 0;

  double get_U() const { return U; }

  double get_acceptance_rate() const {
    if(total_iterations == 0)
      return 0.0;
    return static_cast<double>(accepted_U) / static_cast<double>(total_iterations);
  }

  virtual ~U_sampler() = default;
};