MultiIndexFactory.h

Go to the documentation of this file.

#ifndef MULTIINDEXFACTORY_H_
#define MULTIINDEXFACTORY_H_
 
#include "MUQ/config.h"
 
#include "MUQ/Utilities/MultiIndices/MultiIndexSet.h"
#include "MUQ/Utilities/MultiIndices/MultiIndexLimiter.h"
 
namespace muq {
  namespace Utilities{
 
    class  MultiIndexFactory{
 
    public:
 
      static std::shared_ptr<MultiIndexSet> CreateTotalOrder(unsigned int const length,
                                                             unsigned int const maxOrder,
                                                             unsigned int const minOrder = 0,
                                                             std::shared_ptr<MultiIndexLimiter> limiter = std::make_shared<NoLimiter>());
 
 
      static std::vector<std::shared_ptr<MultiIndexSet>> CreateTriTotalOrder(unsigned int const length,
                                                                             unsigned int const maxOrder,
                                                                             unsigned int const minOrder = 0,
                                                                             std::shared_ptr<MultiIndexLimiter> limiter = std::make_shared<NoLimiter>());
 
      static std::shared_ptr<MultiIndexSet> CreateHyperbolic(unsigned int const length,
                                                             unsigned int const maxOrder,
                                                             const double q = 0.5,
                                                             std::shared_ptr<MultiIndexLimiter> limiter = std::make_shared<NoLimiter>());
 
 
      static std::vector<std::shared_ptr<MultiIndexSet>> CreateTriHyperbolic(unsigned int const length,
                                                                             unsigned int const maxOrder,
                                                                             const double q = 0.5,
                                                                             std::shared_ptr<MultiIndexLimiter> limiter = std::make_shared<NoLimiter>());
 
      static std::shared_ptr<MultiIndexSet> CreateFullTensor(unsigned int const length,
                                                             unsigned int const order,
                                                             std::shared_ptr<MultiIndexLimiter> limiter = std::make_shared<NoLimiter>());
 
      static std::shared_ptr<MultiIndexSet> CreateFullTensor(const Eigen::RowVectorXi& orders,
                                                             std::shared_ptr<MultiIndexLimiter> limiter = std::make_shared<NoLimiter>());
 
      static std::shared_ptr<MultiIndexSet> CreateAnisotropic(const Eigen::RowVectorXf& weights, const double epsilon);
 
 
      static std::shared_ptr<MultiIndex> CreateSingleTerm(int totalDim, int nonzeroDim, int order);
 
 
    private:
 
      static void RecursiveHyperbolicFill(const double maxOrderPow,
                                          std::shared_ptr<MultiIndexSet> output,
                                          unsigned int const currDim,
                                          Eigen::RowVectorXi &base,
                                          const double q,
                                          std::shared_ptr<MultiIndexLimiter> limiter);
 
 
      static void RecursiveTotalOrderFill(unsigned int const maxOrder,
                                          unsigned int const minOrder,
                                          std::shared_ptr<MultiIndexSet> output,
                                          unsigned int const currDim,
                                          Eigen::RowVectorXi &base,
                                          std::shared_ptr<MultiIndexLimiter> limiter);
 
 
      static void RecursiveTensor(const Eigen::RowVectorXi& orders,
                                  std::shared_ptr<MultiIndexSet> output,
                                  unsigned int const currDim,
                                  Eigen::RowVectorXi &base,
                                  std::shared_ptr<MultiIndexLimiter> limiter,
                                  bool allInactive);
 
    };// class MultiIndexFactory
  }// namespace Utilities
}// namespace muq
 
 
#endif