IndexedScalarBasis.h

Go to the documentation of this file.

#ifndef INDEXEDSCALARBASIS_H_
#define INDEXEDSCALARBASIS_H_
 
#include <functional>
#include <string>
 
#include "MUQ/Modeling/WorkPiece.h"
#include "MUQ/Utilities/RegisterClassName.h"
 
namespace muq {
  namespace Approximation {
 
    class IndexedScalarBasis : public muq::Modeling::WorkPiece {
    public:
 
 
      typedef std::function<std::shared_ptr<IndexedScalarBasis>()> ScalarBasisConstructorType;
      typedef std::map<std::string, ScalarBasisConstructorType> ScalarBasisMapType;
      static std::shared_ptr<ScalarBasisMapType> GetScalarBasisMap();
 
      // Factory method for consttructing polynomial from family
      static std::shared_ptr<IndexedScalarBasis> Construct(std::string const& polyName);
 
 
      IndexedScalarBasis();
 
      virtual ~IndexedScalarBasis() = default;
 
 
      virtual double BasisEvaluate(int    const order,
                                   double const x) const = 0;
 
 
      virtual Eigen::VectorXd EvaluateAllTerms(int    const maxOrder,
                                               double const x) const
      {
        Eigen::VectorXd output(maxOrder+1);
        for(int i=0; i<=maxOrder; ++i)
          output(i) = BasisEvaluate(i,x);
        return output;
      };
 
      virtual double DerivativeEvaluate(int    const polyOrder,
                                        int    const derivOrder,
                                        double const x) const = 0;
 
 
    private:
 
      virtual void EvaluateImpl(muq::Modeling::ref_vector<boost::any> const& inputs);
 
 
    };
  } // namespace Approximation
} // namespace muq
 
 
#define REGISTER_SCALARBASIS_FAMILY(NAME) static auto regScalarBasis ##NAME     \
    = muq::Approximation::IndexedScalarBasis::GetScalarBasisMap()->insert(std::make_pair(#NAME, muq::Utilities::shared_factory<NAME>()));
 
 
 
#endif