IntegrationRules.hpp

Go to the documentation of this file.

/**
 * \file IntegrationRules.hpp
 *
 * Integration rules
 *
 */
 
#pragma once
 
 
namespace MoFEM {
 
namespace IntRules {
 
namespace NCC {
 
#include <triangle_ncc_rule.h>
#include <tetrahedron_ncc_rule.h>
 
} // namespace NCC
 
namespace NCO {
 
#include <triangle_nco_rule.h>
#include <tetrahedron_nco_rule.h>
 
} // namespace NCO
 
/** Xiao--Gimbutas quadrature rules from quadraturerules.org Q000002. */
namespace XiaoGimbutas {
 
/**
 * Raw quadrature rule data.
 *
 * Points are stored point-by-point in barycentric coordinates. The weights are
 * normalised for a domain of volume one; multiply them by the physical domain
 * volume when evaluating an integral.
 */
struct Rule {
  unsigned int numPoints;
  unsigned int numBarycentricCoordinates;
  const double *points;
  const double *weights;
};
 
enum {
  triangleRuleCount = 30,
  tetrahedronRuleCount = 15,
};
 
/// Rule of polynomial order n is stored at index n - 1.
extern const Rule triangleRules[triangleRuleCount];
 
/// Rule of polynomial order n is stored at index n - 1.
extern const Rule tetrahedronRules[tetrahedronRuleCount];
 
/**
 * Get a triangle rule for the requested polynomial order.
 *
 * Order zero uses the order-one rule. A null pointer is returned for a
 * negative order or an order greater than triangleRuleCount.
 */
const Rule *getTriangleRule(int order);
 
/**
 * Get a tetrahedron rule for the requested polynomial order.
 *
 * Order zero uses the order-one rule. A null pointer is returned for a
 * negative order or an order greater than tetrahedronRuleCount.
 */
const Rule *getTetrahedronRule(int order);
 
} // namespace XiaoGimbutas
 
} // namespace IntRules
 
}