lehrfempp/quad__rule__cache_8cc_source.html

#include "quad_rule_cache.h"


#include "make_quad_rule.h"


namespace lf::quad {


const QuadRule& QuadRuleCache::Get(base::RefEl ref_el, unsigned degree) const {

  LF_ASSERT_MSG(degree >= 0, "degree must be non-negative.");

  auto& vector = cache_[ref_el.Id()];

  if (vector.size() < degree + 1) {

    vector.resize(degree + 1);

  }

  auto& qr = vector[degree];

  if (qr.NumPoints() == 0) {

    qr = make_QuadRule(ref_el, degree);

  }

  LF_ASSERT_MSG(qr.Degree() >= degree, "something is horribly wrong.");

  return qr;

}


}  // namespace lf::quad

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

lf::base::RefEl::Id
constexpr unsigned int Id() const
Return a unique id for this reference element.
Definition ref_el.h:494

lf::quad::QuadRuleCache::Get
const QuadRule & Get(base::RefEl ref_el, unsigned degree) const
Retrieve a quadrature rule for reference element ref_el with degree d. Is identical to calling make_Q...
Definition quad_rule_cache.cc:14

lf::quad::QuadRuleCache::cache_
std::array< std::deque< QuadRule >, 5 > cache_
Definition quad_rule_cache.h:70

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_QuadRule
QuadRule make_QuadRule(base::RefEl ref_el, unsigned degree)
Returns a QuadRule object for the given Reference Element and Degree.
Definition make_quad_rule.cc:21