KEAST
Quadrature Rules for a Tetrahedron

KEAST is a MATLAB library which defines a set of quadrature rules for the tetrahedron.

Related Data and Programs:

DUNAVANT is a MATLAB library of routines for defining Dunavant rules for quadrature on a triangle.

FEKETE is a MATLAB library of routines for defining Fekete rules for quadrature or interpolation on a triangle.

FELIPPA is a MATHEMATICA library which defines quadrature rules for lines, triangles, quadrilaterals, pyramids, wedges, tetrahedrons and hexahedrons.

GM_RULES is a MATLAB library of routines for defining a Grundmann-Moeller rule for quadrature over a triangle, tetrahedron, or general M-dimensional simplex.

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

NCC_TETRAHEDRON is a MATLAB library defining Newton-Cotes closed quadrature rules on a tetrahedron.

NCO_TETRAHEDRON is a MATLAB library defining Newton-Cotes open quadrature rules on a tetrahedron.

NINTLIB is a MATLAB library containing a variety of routines for numerical estimation of integrals in multiple dimensions.

QUADRATURE_RULES_TET is a dataset directory of triples of files defining various quadrature rules on tetrahedrons.

QUADRULE is a MATLAB library which includes a library of routines for defining quadrature rules on a variety of intervals with different weight functions.

STROUD is a MATLAB library containing quadrature rules for a variety of unusual areas, surfaces and volumes in 2D, 3D and N-dimensions.

WANDZURA is a MATLAB library of routines for defining Wandzura rules for quadrature on a triangle.

Reference:

1. Patrick Keast,
   Moderate Degree Tetrahedral Quadrature Formulas,
   Computer Methods in Applied Mechanics and Engineering,
   Volume 55, Number 3, May 1986, pages 339-348.

Tar File:

A GZIP'ed TAR file of the contents of this directory is available. This is only done as a convenience for users who want ALL the files, and don't want to download them individually. This is not a convenience for me, so don't be surprised if the tar file is somewhat out of date.

• keast_m_src.tar.csh the commands used to create the GZIP'ed TAR file.
• keast_m_src.tar.gz the GZIP'ed TAR file.

Source Code:

• comp_next.m returns the next composition of an integer into a given number of parts.
• i4_modp.m returns the nonnegative remainder of integer division.
• i4_wrap.m forces an integer to lie between given limits by wrapping.
• keast_degree.m returns the degree of a Keast rule for the tetrahedron.
• keast_order_num.m returns the order of a Keast rule for the tetrahedron.
• keast_rule.m returns the points and weights of a Keast rule.
• keast_rule_num.m returns the number of Keast rules for the tetrahedron.
• keast_suborder.m returns the suborders for a Keast rule.
• keast_suborder_num.m returns the number of suborders for a Keast rule.
• keast_subrule.m returns a compressed Keast rule.
• monomial_value.m evaluates a monomial function.
• r8mat_det_4d.m computes the determinant of a 4 by 4 R8MAT.
• tetrahedron_reference_to_physical.m maps T4 reference points to physical points.
• tetrahedron_volume.m computes the volume of a tetrahedron in 3D.
• timestamp.m prints the current YMDHMS date as a time stamp.

Examples and Tests:

• keast_test.m, calls all the tests.
• keast_test01.m, tests KEAST_RULE_NUM, KEAST_DEGREE, KEAST_ORDER_NUM.
• keast_test02.m, tests KEAST_RULE.
• keast_test03.m, tests KEAST_RULE.
• keast_test04.m, tests TETRAHEDRON_REFERENCE_TO_PHYSICAL.
• keast_test05.m, applies all the rules to a set of monomial integrands.
• keast_test07.m, prints out a rule.
• keast_test.out, the output from a run of the sample program.

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