InverseGamma.h

Go to the documentation of this file.

#ifndef INVERSEGAMMA_H_
#define INVERSEGAMMA_H_
 
#include "MUQ/Modeling/Distributions/Distribution.h"
 
namespace muq {
  namespace Modeling {
 
    class InverseGamma : public Distribution {
    public:
 
      InverseGamma(double       alphaIn,
                   double       betaIn);
 
      InverseGamma(Eigen::VectorXd const& alphaIn,
                   Eigen::VectorXd const& betaIn);
 
 
      virtual ~InverseGamma() = default;
 
 
      const Eigen::VectorXd alpha;
      const Eigen::VectorXd beta;
 
    private:
 
      static double ComputeConstant(Eigen::VectorXd const& alphaIn,
                                    Eigen::VectorXd const& betaIn);
 
      const double logConst; // sum( log( beta^alpha / Gamma(alpha) ) )
 
      virtual double LogDensityImpl(ref_vector<Eigen::VectorXd> const& inputs) override;
 
      virtual Eigen::VectorXd SampleImpl(ref_vector<Eigen::VectorXd> const& inputs) override;
 
    };
  } // namespace Modeling
} // namespace muq
 
#endif