GaussianMarginal.h

Go to the documentation of this file.

#ifndef STANDARD_NORMAL_H
#define STANDARD_NORMAL_H
 
#include <random>
 
#include "IMarginalDistribution.h"
#include "AbstractConfig.h"
#include "Math.h"
 
using gsl_rng_sptr = std::unique_ptr<gsl_rng, decltype(&gsl_rng_free)>;
 
struct GaussianMarginalCfg : public AbstractConfig {
    double mu {0.};     
    double sigma {1.};  
 
    GaussianMarginalCfg() = default;
 
    GaussianMarginalCfg(double mu, double sigma) : mu(mu), sigma(sigma) {}
};
 
class GaussianMarginal final : public IMarginalDistribution {
public:
    explicit GaussianMarginal(double mu, double sigma, 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 mu, sigma;                               
    const gsl_rng_type* rng_tp {gsl_rng_mt19937};   
    gsl_rng_sptr eng_{gsl_rng_alloc(rng_tp), &gsl_rng_free};      
};
 
#endif