TRIANGLE_SYMQ_RULE
Symmetric Quadrature Rules for Triangles.


TRIANGLE_SYMQ_RULE is a FORTRAN90 library which returns symmetric quadrature rules, with exactness up to total degree 50, over the interior of an arbitrary triangle in 2D, by Hong Xiao and Zydrunas Gimbutas.

The original source code, from which this library was developed, is available from the Courant Mathematics and Computing Laboratory, at http://www.cims.nyu.edu/cmcl/quadratures/quadratures.html ,

Licensing:

The computer code and data files made available on this web page are distributed under the GNU LGPL license.

Languages:

TRIANGLE_SYMQ_RULE is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version.

Related Data and Programs:

CUBE_FELIPPA_RULE, a FORTRAN90 library which returns the points and weights of a Felippa quadrature rule over the interior of a cube in 3D.

GNUPLOT, FORTRAN90 programs which illustrate how a program can write data and command files so that gnuplot can create plots of the program results.

PYRAMID_FELIPPA_RULE, a FORTRAN90 library which returns Felippa's quadratures rules for approximating integrals over the interior of a pyramid in 3D.

SIMPLEX_GM_RULE, a FORTRAN90 library which defines Grundmann-Moeller quadrature rules over the interior of a simplex in M dimensions.

SQUARE_FELIPPA_RULE, a FORTRAN90 library which returns the points and weights of a Felippa quadrature rule over the interior of a square in 2D.

STROUD, a FORTRAN90 library which defines quadrature rules for a variety of M-dimensional regions, including the interior of the square, cube and hypercube, the pyramid, cone and ellipse, the hexagon, the M-dimensional octahedron, the circle, sphere and hypersphere, the triangle, tetrahedron and simplex, and the surface of the circle, sphere and hypersphere.

TETRAHEDRON_ARBQ_RULE, a FORTRAN90 library which returns quadrature rules, with exactness up to total degree 15, over the interior of a tetrahedron in 3D, by Hong Xiao and Zydrunas Gimbutas.

TETRAHEDRON_FELIPPA_RULE, a FORTRAN90 library which returns Felippa's quadratures rules for approximating integrals over the interior of a tetrahedron in 3D.

TOMS886, a FORTRAN90 library which defines the Padua points for interpolation in a 2D region, including the rectangle, triangle, and ellipse, by Marco Caliari, Stefano de Marchi, Marco Vianello. This is a version of ACM TOMS algorithm 886.

TRIANGLE_DUNAVANT_RULE, a FORTRAN90 library which sets up a Dunavant quadrature rule over the interior of a triangle in 2D.

TRIANGLE_FELIPPA_RULE, a FORTRAN90 library which returns Felippa's quadratures rules for approximating integrals over the interior of a triangle in 2D.

TRIANGLE_INTEGRALS, a FORTRAN90 library which returns the exact value of the integral of any monomial over the interior of the unit triangle in 2D.

TRIANGLE_LYNESS_RULE, a FORTRAN90 library which returns Lyness-Jespersen quadrature rules over the interior of a triangle in 2D.

TRIANGLE_MONTE_CARLO, a FORTRAN90 library which uses the Monte Carlo method to estimate the integral of a function over the interior of the unit triangle in 2D.

TRIANGLE_NCC_RULE, a FORTRAN90 library which defines Newton-Cotes Closed (NCC) quadrature rules over the interior of a triangle in 2D.

TRIANGLE_NCO_RULE, a FORTRAN90 library which defines Newton-Cotes Open (NCO) quadrature rules over the interior of a triangle in 2D.

TRIANGLE_WANDZURA_RULE, a FORTRAN90 library which sets up a quadrature rule of exactness 5, 10, 15, 20, 25 or 30 over the interior of a triangle in 2D.

WEDGE_FELIPPA_RULE, a FORTRAN90 library which returns quadratures rules for approximating integrals over the interior of the unit wedge in 3D.

Reference:

  1. Hong Xiao, Zydrunas Gimbutas,
    A numerical algorithm for the construction of efficient quadrature rules in two and higher dimensions,
    Computers and Mathematics with Applications,
    Volume 59, 2010, pages 663-676.

Source Code:

Examples and Tests:

EQUI08 is a degree 8 rule in an equilateral triangle.

SIMP08 is a degree 8 rule in a simplex.

USER08 is a degree 8 rule in a user specified triangle at (1,0), (4,4), (0,3).

List of Routines:

You can go up one level to the FORTRAN90 source codes.


Last revised on 26 June 2014.