exporter.hpp

Go to the documentation of this file.

#ifndef EXPORTED_HPP
#define EXPORTED_HPP
 
#include <Eigen/Dense>
#include <vector>
 
#include "exterior_dimension.hpp"
#include "mesh.hpp"
 
namespace Manicore {
 
 
  template<size_t dimension>
  class Exporter {
    public:
 
      Exporter(Mesh<dimension> const * mesh ,
               int r ,
               int acc 
              );
 
 
      int save(size_t k , 
                std::function<Eigen::VectorXd(size_t iT)> Fu_h ,
                const char *filename ,
                bool star=false ) const;
 
      int saveSq(size_t k , 
                    std::function<Eigen::VectorXd(size_t iT)> Fu_h ,
                    const char *filename ) const;
 
    private:
      std::vector<Eigen::Vector3d> _locations;
      std::vector<size_t> _cellIndices;
      std::vector<Eigen::Matrix<double,dimension,dimension>> _metricInv;
      std::vector<double> _volume;
      std::vector<Eigen::VectorXd> _polyEval;
      // TODO Make it works in arbitrary dimension. What to do with forms not match a scalar or a vector?
      std::vector<Eigen::Matrix<double,Dimension::ExtDim(1,3),Dimension::ExtDim(1,dimension)>> _pushforward;
      std::vector<Eigen::Matrix<double,Dimension::ExtDim(1,3),Dimension::ExtDim(dimension-1,dimension)>> _pushforwardStar;
      /* // Structure to store various size of pullback matrices
      template<size_t _k> struct PullbackHolder {
        std::vector<Eigen::Matrix<double,Dimension::ExtDim(_k,3),Dimension::ExtDim(_k,dimension)>> _pullback;
        PullbackHolder<_k-1> _next;
        template<size_t k> std::vector<Eigen::Matrix<double,Dimension::ExtDim(k,3),Dimension::ExtDim(k,dimension)>> & pullback() {
          if constexpr (k==_k) return _pullback;
          else return _next.template pullback<k>();
        }
        template<size_t k> std::vector<Eigen::Matrix<double,Dimension::ExtDim(k,3),Dimension::ExtDim(k,dimension)>> pullback() const {
          if constexpr (k==_k) return _pullback;
          else return _next.template pullback<k>();
        }
      };
      template<> struct PullbackHolder<0> {
        std::vector<Eigen::Matrix<double,1,1>> _pullback;
        template<size_t k> std::vector<Eigen::Matrix<double,1,1>> & pullback() {
          static_assert(k == 0);
          return _pullback;
        }
        template<size_t k> std::vector<Eigen::Matrix<double,1,1>> pullback() const {
          static_assert(k == 0);
          return _pullback;
        }
      };
      PullbackHolder<dimension> _pullbackEval;
      */ 
  };
}
#endif