wedge_felippa_rule

wedge_felippa_rule, a Fortran77 code which generates Felippa quadrature rules 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

f77_rule, a Fortran77 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.f, the source code.
• wedge_felippa_rule.sh, commands to compile the source code.