dcell_integral.hpp

Go to the documentation of this file.

#ifndef DCELL_INTEGRAL_HPP_INCLUDED
#define DCELL_INTEGRAL_HPP_INCLUDED
 
#include "integral.hpp"
 
namespace Manicore {
 
 
  template<size_t dimension,size_t d> requires(d > 0 && d <= dimension)
  struct dCell_mass {
      dCell_mass(size_t i_cell , int r , const QuadratureRule<d> & quad , const Integral<dimension,d> &integral );
      dCell_mass() {;}
 
 
      std::array<Eigen::MatrixXd,dimension+1> masses;
  };
 
 
  template<size_t dimension,size_t d> requires(d > 0 && d <= dimension)
  struct dCell_traces {
      dCell_traces(size_t i_cell , int r , int dqr , 
                   const std::vector<dCell_mass<dimension,d-1>> & b_masses ,
                   const Integral<dimension,d> &integral , 
                   const Integral<dimension,d-1> &integral_b );
      dCell_traces() {;}
 
      // The array index is the form degree (0 to d-1), and the vector index is the boundary
 
      std::array<std::vector<Eigen::MatrixXd>,dimension> traces;
      std::array<std::vector<Eigen::MatrixXd>,dimension> starTraces; // TODO: only used for the stabilization term in the L2-product, check if this is really necessary
  };
  template<size_t dimension> 
  struct dCell_traces<dimension,1> {
      dCell_traces(size_t i_cell , int r , const Mesh<dimension>* mesh );
      dCell_traces() {;}
 
      // The array index is the form degree (0 to d-1), and the vector index is the boundary
 
      std::array<std::vector<Eigen::MatrixXd>,dimension> traces;
      std::array<std::vector<Eigen::MatrixXd>,dimension> starTraces;
  };
}
 
#endif