quadraturerule.hpp

Go to the documentation of this file.

//
// 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>
//#include <cell.hpp>
//#include <edge.hpp>
 
namespace HArDCore2D {
 
  struct QuadratureNode
  {
    double x, y, w;
    QuadratureNode(double x, double y, double w)
      : x(x), y(y), w(w)
    {
      // Do nothing
    }
 
    inline Eigen::Vector2d vector() const {
      return Eigen::Vector2d(x,y);
    }
  };
 
  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 Edge & E,                
                      const int doe               
                      );  
}  // end of namespace HArDCore2D
#endif