UniformBox.h

Go to the documentation of this file.

#ifndef UNIFORMBOX_H_
#define UNIFORMBOX_H_
 
#include "MUQ/Utilities/VariadicMacros.h"
#include "MUQ/Utilities/Exceptions.h"
 
#include "MUQ/Modeling/Distributions/Distribution.h"
 
 
namespace muq {
  namespace Modeling {
 
    class UniformBox : public Distribution {
    public:
 
      virtual ~UniformBox() = default;
 
 
      UniformBox(Eigen::MatrixXd const& bounds);
 
 
      template<typename... Args>
      inline UniformBox(std::pair<double,double> const& pair1, Args... args) : UniformBox(CreateBoundsPairs(pair1, args...)) {}
 
      inline UniformBox(std::vector<std::pair<double,double>> const& pairList) : UniformBox(CreateBoundsPairs(pairList)) {}
 
      template<typename... Args>
      inline UniformBox(double lb1, double ub1, Args... args) : UniformBox(CreateBoundsDouble(lb1, ub1, args...)) {}
 
      inline UniformBox(std::vector<double> const& bounds) : UniformBox(CreateBoundsDouble(bounds)) {}
 
    private:
 
 
      virtual double LogDensityImpl(ref_vector<Eigen::VectorXd> const& inputs) override;
 
      virtual Eigen::VectorXd SampleImpl(ref_vector<Eigen::VectorXd> const& inputs) override;
 
 
      //static Eigen::MatrixXd CreateBounds(std::vector<double>& bounds);
 
      static Eigen::MatrixXd CreateBoundsPairs(std::vector<std::pair<double,double>> const& bounds);
      STATIC_VARIADIC_TO_VECTOR(CreateBoundsPairs, (std::pair<double, double>), (Eigen::MatrixXd))
 
      static Eigen::MatrixXd CreateBoundsDouble(std::vector<double> const& bounds);
      STATIC_VARIADIC_TO_VECTOR(CreateBoundsDouble, (double), (Eigen::MatrixXd))
 
          //static Eigen::MatrixXd CreateBounds(std::vector<std::pair<double,double>>& bounds);
          //STATIC_VARIADIC_TO_VECTOR_PART2(CreateBounds, double, (Eigen::MatrixXd))
 
 
      
      const Eigen::MatrixXd bounds;
 
      // Use the values stored in bounds to compute the volume
      static double ComputeVolume(Eigen::MatrixXd const& boundsIn);
 
      // The volume of the region enclosed by the bounds
      const double volume;
 
    };
  } // namespace Modeling
} // namespace muq
 
#endif