quadraturerule.hpp

Go to the documentation of this file.

// Core data structures and methods required to implement the discrete de Rham sequence in 3D
//
// Provides:
//  - Generation of quadrature rules over cells and faces of the mesh
//
// Author: Daniele Di Pietro (daniele.di-pietro@umontpellier.fr)
// Author: Jérome Droniou (jerome.droniou@monash.edu)
//
 
/*
 *
 *      This library was developed around HHO methods, although some parts of it have a more
 * general purpose. If you use this code or part of it in a scientific publication, 
 * please mention the following book as a reference for the underlying principles
 * of HHO schemes:
 *
 * The Hybrid High-Order Method for Polytopal Meshes: Design, Analysis, and Applications. 
 *  D. A. Di Pietro and J. Droniou. Modeling, Simulation and Applications, vol. 19. 
 *  Springer International Publishing, 2020, xxxi + 525p. doi: 10.1007/978-3-030-37203-3. 
 *  url: https://hal.archives-ouvertes.fr/hal-02151813.
 *
 */
 
#ifndef QUADRATURERULE_HPP
#define QUADRATURERULE_HPP
 
#include <vector>
 
#include <mesh.hpp>
 
 
namespace HArDCore3D {
 
  struct QuadratureNode
  {
    double x, y, z, w;
    QuadratureNode(double x, double y, double z, double w)
      : x(x), y(y), z(z), w(w)
    {
      // Do nothing
    }
 
    inline Eigen::Vector3d vector() const {
      return Eigen::Vector3d(x,y,z);
    }
  };
 
  typedef std::vector<QuadratureNode> QuadratureRule;
 
  QuadratureRule generate_quadrature_rule(
                      const Cell & T,                
                      const int doe,              
                      const bool force_split = false 
                      );  
  QuadratureRule generate_quadrature_rule(
                      const Face & F,                
                      const int doe               
                      );  
  QuadratureRule generate_quadrature_rule(
                      const Edge & E,                
                      const int doe               
                      );  
}  // end of namespace HArDCore3D
#endif