stroud_rule, a Fortran90 code 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.
A few other rules have been collected as well, particularly for quadrature over the interior of a triangle, which is useful in finite element calculations.
Arthur Stroud published his vast collection of quadrature formulas for multidimensional regions in 1971. In a few cases, he printed sample Fortran77 programs to compute these integrals. Integration regions included:
We have added a few new terms for regions:
The information on this web page is distributed under the MIT license.
stroud_rule is available in a C version and a C++ version and a Fortran77 version and a Fortran90 version and a MATLAB version and an Octave version.
disk_rule, a Fortran90 code which computes quadrature rules for the unit disk in 2D, that is, the interior of the circle of radius 1 and center (0,0).
pyramid_rule, a Fortran90 code which computes a quadrature rule for a pyramid.
simplex_gm_rule, a Fortran90 code which defines Grundmann-Moeller quadrature rules over the interior of a triangle in 2D, a tetrahedron in 3D, or over the interior of the simplex in M dimensions.
sphere_design_rule, a Fortran90 code which returns point sets on the surface of the unit sphere, known as "designs", which can be useful for estimating integrals on the surface.
sphere_lebedev_rule, a Fortran90 code which computes Lebedev quadrature rules on the surface of the unit sphere in 3D.
tetrahedron_arbq_rule, a Fortran90 code 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_keast_rule, a Fortran90 code which defines ten quadrature rules, with exactness degrees 0 through 8, over the interior of a tetrahedron in 3D.
tetrahedron_ncc_rule, a Fortran90 code which defines Newton-Cotes Closed (NCC) quadrature rules over the interior of a tetrahedron in 3D.
tetrahedron_nco_rule, a Fortran90 code which defines Newton-Cotes Open (NCO) quadrature rules over the interior of a tetrahedron in 3D.
triangle_dunavant_rule, a Fortran90 code which defines Dunavant rules for quadrature over the interior of a triangle in 2D.
tetrahedron_keast_rule, a Fortran90 code which defines ten quadrature rules, with exactness degrees 0 through 8, over the interior of a tetrahedron in 3D.
triangle_dunavant_rule, a Fortran90 code which defines Dunavant rules for quadrature over the interior of a triangle in 2D.
triangle_fekete, a Fortran90 code which defines Fekete rules for quadrature or interpolation over the interior of a triangle in 2D.
triangle_lyness_rule, a Fortran90 code which returns Lyness-Jespersen quadrature rules over the interior of a triangle in 2D.
triangle_ncc_rule, a Fortran90 code which defines Newton-Cotes Closed (NCC) quadrature rules over the interior of a triangle in 2D.
triangle_nco_rule, a Fortran90 code which defines Newton-Cotes Open (NCO) quadrature rules over the interior of a triangle in 2D.
triangle_wandzura_rule, a Fortran90 code which returns quadrature rules of exactness 5, 10, 15, 20, 25 and 30 over the interior of the triangle in 2D.