inv_gauss_dist.h

#ifndef _INV_GAUSS_DIST_H_
#define _INV_GAUSS_DIST_H_

#include <string>
#include <vector>
#include <boost/math/distributions/inverse_gaussian.hpp>
#include "distribution.h"

namespace stochastic {
class InverseGaussianDistribution : public Distribution {
 public:
  InverseGaussianDistribution() = delete;

  InverseGaussianDistribution(double mean, double scale);

  virtual ~InverseGaussianDistribution(){};

  std::string name() const override { return "InverseGaussianDist"; };

  std::vector<double> cumulative_dist_func(
      const std::vector<double>& locations) const override;

  std::vector<double> inv_cumulative_dist_func(
      const std::vector<double>& probabilities) const override;

 protected:
  double mean_;    
  double std_dev_; 
  boost::math::inverse_gaussian distribution_; 
};
}  // namespace stochastic

#endif  // _INV_GAUSS_DIST_H_