lehrfempp/make__quad__rule__nodal_8cc_source.html

#include "make_quad_rule_nodal.h"


#include <boost/core/ref.hpp>


namespace lf::quad {


QuadRule make_QuadRuleNodal(base::RefEl ref_el) {

  double ref_el_vol = 1.0;

  if (ref_el == base::RefEl::kTria()) {

    ref_el_vol = 0.5;

  }

  return QuadRule(ref_el, ref_el.NodeCoords(),

                  Eigen::VectorXd::Constant(ref_el.NumNodes(),

                                            ref_el_vol / ref_el.NumNodes()),

                  1);

}


}  // namespace lf::quad

lf::base::RefEl
Represents a reference element with all its properties.
Definition ref_el.h:109

lf::base::RefEl::NodeCoords
const Eigen::MatrixXd & NodeCoords() const
Get the coordinates of the nodes of this reference element.
Definition ref_el.h:241

lf::base::RefEl::kTria
static constexpr RefEl kTria()
Returns the reference triangle.
Definition ref_el.h:161

lf::base::RefEl::NumNodes
constexpr size_type NumNodes() const
The number of nodes of this reference element.
Definition ref_el.h:213

lf::quad::QuadRule
Represents a Quadrature Rule over one of the Reference Elements.
Definition quad_rule.h:58

lf::quad
Rules for numerical quadrature on reference entity shapes.
Definition gauss_quadrature.cc:15

lf::quad::make_QuadRuleNodal
QuadRule make_QuadRuleNodal(base::RefEl ref_el)
Create a quadrature rule that evaluates the quadrand only at the nodes of the reference element.
Definition make_quad_rule_nodal.cc:14