wedge_felippa_rule

wedge_felippa_rule, a Python code which returns 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:

python_rule, a Python 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.py, the source code.
• wedge_felippa_rule.sh, runs all the tests.
• wedge_felippa_rule.txt, the output file.