FlatMarginal.h

Go to the documentation of this file.

#ifndef STANDARD_FLAT_H
#define STANDARD_FLAT_H
 
#include <random>
#include <cmath>
#include <limits>
#include <stdexcept>
 
#include "IMarginalDistribution.h"
#include "Include.h"
#include "AbstractConfig.h"
 
 using gsl_rng_sptr = std::unique_ptr<gsl_rng, decltype(&gsl_rng_free)>;
 
struct FlatMarginalCfg : public AbstractConfig {
    double a;   
    double b;   
 
    FlatMarginalCfg() = default;
 
    FlatMarginalCfg(double a, double b) : a(a), b(b) {}
};
 
class FlatMarginal final : public IMarginalDistribution {
public:
    explicit FlatMarginal(double a, double b, unsigned int seed = std::random_device{}());
 
    std::vector<double> rvs(std::size_t n) override;
 
    double logpdf(double x) override;
 
    PDFDiff f_df_ddf(double x) override;
 
    double cdf(double x) override;
 
    double ppf(double p) override;
 
    double mean() override;
 
    double std() override;
 
private:
    double a, b;                                    
    const gsl_rng_type* rng_tp {gsl_rng_mt19937};   
    gsl_rng_sptr eng_{gsl_rng_alloc(rng_tp), &gsl_rng_free};      
};
 
#endif