definitions.hpp

Go to the documentation of this file.

#ifndef DEFINITIONS_HPP_INCLUDED
#define DEFINITIONS_HPP_INCLUDED
 
#include <Eigen/Dense>
 
#include <array>
 
namespace Manicore {
  
  template<size_t d>
    using Simplex = std::array<Eigen::Vector<double,d>,d+1>;
  
 
  template<size_t dimension,size_t d>
    struct ParametrizedMap {
      virtual Eigen::Vector<double,dimension> I(Eigen::Vector<double,d> const &) const = 0;
      virtual Eigen::Vector<double,d> J(Eigen::Vector<double,dimension> const &) const = 0;
      std::vector<double> _extra;
      virtual ~ParametrizedMap() = default;
    };
 
 
  template<size_t d1,size_t d2>
    struct ParametrizedDerivedMap {
 
      virtual Eigen::Matrix<double,d1,d2> DI(Eigen::Vector<double,d2> const &) const = 0;
 
      virtual Eigen::Matrix<double,d2,d1> DJ(Eigen::Vector<double,d1> const &) const = 0;
      std::vector<double> _extra;
      virtual ~ParametrizedDerivedMap() = default;
    };
 
 
  template<size_t dimension>
    struct ParametrizedMetricMap {
      virtual Eigen::Matrix<double,dimension,dimension> metric(Eigen::Vector<double,dimension> const &) const = 0;
 
      virtual Eigen::Matrix<double,dimension,dimension> metric_inv(Eigen::Vector<double,dimension> const &) const = 0;
 
      virtual double volume(Eigen::Vector<double,dimension> const &) const = 0;
      std::vector<double> _extra;
      ParametrizedMetricMap(int o) : orientation(o) {;}
 
      const int orientation = 0;
      virtual ~ParametrizedMetricMap() = default;
    };
 
  template<size_t dimension,size_t d>
    ParametrizedMap<dimension,d> *DefaultMapping(size_t id) {return nullptr;}
  template<size_t dimension,size_t d>
    ParametrizedDerivedMap<dimension,d> *DefaultDerivedMapping(size_t id) {return nullptr;}
 
} // namespace
 
#endif