square_felippa_rule

square_felippa_rule, a Fortran77 code which generates the points and weights of a Felippa quadrature rule over the interior of a square in 2D.

Actually, the word "square" is meant to designate any quadrature region defined by:

        A(1) <= X <= B(1)
        A(2) <= Y <= B(2)
      

Licensing:

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

Languages:

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

Related Data and Programs:

square_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:

• square_felippa_rule.f, the source code.
• square_felippa_rule.sh, BASH commands to compile the source code.