Quadrature.h

Go to the documentation of this file.

#ifndef QUADRATURE_H_
#define QUADRATURE_H_
 
#include <memory>
 
#include <Eigen/Core>
 
namespace muq {
  namespace Approximation {
 
    class Quadrature {
    public:
 
      Quadrature(unsigned int dimIn) : dim(dimIn){};
 
 
      virtual ~Quadrature() = default;
 
      virtual void Compute(unsigned int quadOrder) = 0;
 
      virtual unsigned int Exactness(unsigned int quadOrder) const;
 
      virtual void Compute(Eigen::RowVectorXi const& orders) { assert(orders.size()==1); Compute(orders(0));}
 
      virtual unsigned int Dim() const{return dim;};
 
      virtual Eigen::MatrixXd const& Points() const{return pts;};
 
      virtual Eigen::VectorXd const& Weights() const{return wts;};
 
 
    protected:
      unsigned int dim;
      Eigen::MatrixXd pts;
      Eigen::VectorXd wts;
    };
  }
}
 
#endif /* QUADRATURE_H_ */