wedge_felippa_rule

wedge_felippa_rule, a C++ code which generates a Felippa quadrature rule to estimate integrals over the interior of the unit wedge in 3D.

The interior of the unit wedge in 3D is defined by the constraints:

        0 <= X
        0 <= Y
             X + Y <= 1
       -1 <= Z <= +1
      

Licensing:

The information on this web page is distributed under the MIT license.

Languages:

wedge_felippa_rule is available in a C version and a C++ version and a Fortran90 version and a MATLAB version and an Octave version and a Python version.

Related Data and Programs:

wedge_felippa_rule_test

cpp_rule, a C++ code which computes a quadrature rule which estimates the integral of a function f(x), which might be defined over a one dimensional region (a line) or more complex shapes such as a circle, a triangle, a quadrilateral, a polygon, or a higher dimensional region, and which might include an associated weight function w(x).

Reference:

1. Carlos Felippa,
   A compendium of FEM integration formulas for symbolic work,
   Engineering Computation,
   Volume 21, Number 8, 2004, pages 867-890.

Source Code:

• wedge_felippa_rule.cpp, the source code.
• wedge_felippa_rule.sh, compiles the source code.
• wedge_felippa_rule.hpp, the include file.