TEST_TRI_INT
Triangle Integration Test Functions

TEST_TRI_INT is a library of C++ routines, using double precision arithmetic, which defines a number of test functions to be integrated over a triangle.

It is possible to invoke a particular function by number, or to try out all available functions, as demonstrated in the sample calling program.

For convenience, all the integrand functions have been scaled by a constant, so that the integral of the function over the specific domain is exactly 1.

The test functions include F(X,Y)=

1. 1 on the unit triangle.
2. x on the unit triangle.
3. y on the unit triangle.
4. x² on the unit triangle.
5. x*y on the unit triangle.
6. y² on the unit triangle.
7. x³ on the unit triangle.
8. x⁴ on the unit triangle.
9. x⁵ on the unit triangle.
10. x^(-0.2) remapped to (1,0), (5,0), (5,1).
11. (x+y)^(-0.2) on the unit triangle.
12. (1-x-y)^(-0.2) remapped to (-1,-3), (3,-2), (-1,2).
13. (x*y)^(-0.2) remapped to (0,0), (-7,0), (0,-3).
14. 1/sqrt(x) + 1/sqrt(y) + 1/sqrt(x+y) on the unit triangle.
15. 1/sqrt(1-x-y) on the unit triangle.
16. log(x*y) on the unit triangle.
17. 1/sqrt(|x-1/4|) + 1/sqrt(|y-1/2|) on the unit triangle.
18. log ( x + y ) on the unit triangle.
19. sin ( x ) cos ( 5 y ) on the unit triangle.
20. sin ( 11 x ) cos ( y ) on the unit triangle.
21. 1 / r on the unit triangle,
    r = sqrt ( x^2+y^2).
22. log ( r ) / r on the unit triangle,
    r = sqrt ( x^2+y^2).

The library includes a routines to define the integrand function, the triangle over which the integral is to be carried out, and a title for the problem. Thus, for each integrand function, four routines are supplied. For instance, for function #4, we have the routines:

• P04_FUN evaluates the integrand for problem 4.
• P04_TITLE returns the title for problem 4.
• P04_VERTICES returns the integration triangle for problem 4

Moreover, since the same interface is used for each function, if you wish to work with problem 16 instead, you simply change the "04" to "07" in your routine calls.

If you wish to call all of the functions, then you simply use the generic interface, which again has four routines, but which requires you to specify the problem number as an extra input argument:

• P00_FUN evaluates the integrand for problem PROB.
• P00_SINGULARITY warns if the integran for problem PROB has vertex, edge or interior singularities that might cause problems for some integration schemes.
• P00_TITLE returns the title for problem PROB.
• P00_VERTICES returns the integration triangle for problem PROB.

Finally, some demonstration routines are built in for simple quadrature methods. These routines include

• P00_WANDZURA05_SUB applies an order-5 accurate Wandzura rule, plus subdivision of the triangle.
• P00_MONTE_CARLO uses Monte Carlo sampling.
• P00_VERTEX_SUB averages values at vertices, and uses subdivision.

Related Data and Programs:

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

FEKETE is a C++ library of routines for defining a Fekete rule for quadrature or interpolation over a triangle.

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

INTLIB is a C++ library of routines for numerically estimating integrals in one dimension.

NCC_TRIANGLE is a C++ library defining Newton-Cotes closed quadrature rules on a triangle.

NCO_TRIANGLE is a C++ library defining Newton-Cotes open quadrature rules on a triangle.

NINT_EXACTNESS_TRI is an executable C++ program which investigates the polynomial exactness of a quadrature rule for the triangle.

NINTLIB is a C++ library of routines for numerically estimating integrals in multiple dimensions.

QUADPACK is a FORTRAN90 library of routines for estimating integrals in one dimension.

QUADRULE is a C++ library of routines for defining quadrature rules on a variety of intervals with different weight functions.

STRI_QUAD is a FORTRAN90 library of routines which can approximate the integral of a function over the surface of a sphere.

STROUD is a FORTRAN90 library of routines which includes some quadrature rules for triangles.

TEST_NINT is a C++ library of functions that can be used to test N-dimensional quadrature routines.

TEST_TRI_INT is also available in a FORTRAN90 version and a MATLAB version.

TOMS612 is a FORTRAN77 library of routines which can estimate the integral of a function over a triangle.

TOMS706 is a FORTRAN77 library which estimates the integral of a function over a triangulated region.

TRIANGULATION is a C++ library of routines for triangulations, which includes some quadrature rules on triangles.

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

Reference:

1. Elise deDoncker, Ian Robinson,
   An Algorithm for Automatic Integration Over a Triangle Using Nonlinear Extrapolation,
   ACM Transactions on Mathematical Software,
   Volume 10, Number 1, March 1984, pages 1-16.
2. Elise deDoncker, Ian Robinson,
   Algorithm 612: Integration over a Triangle Using Nonlinear Extrapolation,
   ACM Transactions on Mathematical Software,
   Volume 10, Number 1, March 1984, pages 17-22.
3. Stephen Wandzura, Hong Xiao,
   Symmetric Quadrature Rules on a Triangle,
   Computers and Mathematics with Applications,
   Volume 45, pages 1829-1840, 2003.

Source Code:

• test_tri_int.C, the source code.
• test_tri_int.H, the source code.
• test_tri_int.csh, commands to compile the source code.

Examples and Tests:

• test_tri_int_prb.C, a sample calling program.
• test_tri_int_prb.csh, commands to compile and run the sample program.
• test_tri_int_prb.out, the output from a run of the sample program.

List of Routines:

• GET_PROBLEM_NUM returns the number of test integration problems.
• GET_SEED returns a random seed for the random number generator.
• I4_HUGE returns a "huge" I4.
• I4_POWER returns the value of I^J.
• P00_FUN evaluates the integrand for any problem.
• P00_MONTE_CARLO applies the Monte Carlo rule to integrate a function.
• P00_SINGULARITY warns of common singularities for any problem.
• P00_TITLE returns the title for any problem.
• 00_VERTEX_SUB approximates an integral in a triangle by subdivision.
• P00_VERTICES returns the vertices for any problem.
• P00_WANDZURA05_SUB uses subdivision and a Wandzura rule.
• P01_FUN evaluates the integrand for problem 1.
• P01_TITLE returns the name of problem 1.
• P01_VERTICES returns the vertices for problem 1.
• P02_FUN evaluates the integrand for problem 2.
• P02_TITLE returns the name of problem 2.
• P02_VERTICES returns the vertices for problem 2.
• P03_FUN evaluates the integrand for problem 3.
• P03_TITLE returns the name of problem 3.
• P03_VERTICES returns the vertices for problem 3.
• P04_FUN evaluates the integrand for problem 4.
• P04_TITLE returns the name of problem 4.
• P04_VERTICES returns the vertices for problem 4.
• P05_FUN evaluates the integrand for problem 5.
• P05_TITLE returns the name of problem 5.
• P05_VERTICES returns the vertices for problem 5.
• P06_FUN evaluates the integrand for problem 6.
• P06_TITLE returns the name of problem 6.
• P06_VERTICES returns the vertices for problem 6.
• P07_FUN evaluates the integrand for problem 7.
• P07_TITLE returns the name of problem 7.
• P07_VERTICES returns the vertices for problem 7.
• P08_FUN evaluates the integrand for problem 8.
• P08_TITLE returns the name of problem 8.
• P08_VERTICES returns the vertices for problem 8.
• P09_FUN evaluates the integrand for problem 9.
• P09_TITLE returns the name of problem 9.
• P09_VERTICES returns the vertices for problem 9.
• P10_FUN evaluates the integrand for problem 10.
• P10_TITLE returns the name of problem 10.
• P10_VERTICES returns the vertices for problem 10.
• P11_FUN evaluates the integrand for problem 11.
• P11_TITLE returns the name of problem 11.
• P11_VERTICES returns the vertices for problem 11.
• P12_FUN evaluates the integrand for problem 12.
• P12_TITLE returns the name of problem 12.
• P12_VERTICES returns the vertices for problem 12.
• P13_FUN evaluates the integrand for problem 13.
• P13_TITLE returns the name of problem 13.
• P13_VERTICES returns the vertices for problem 13.
• P14_FUN evaluates the integrand for problem 14.
• P14_TITLE returns the name of problem 14.
• P14_VERTICES returns the vertices for problem 14.
• P15_FUN evaluates the integrand for problem 15.
• P15_TITLE returns the name of problem 15.
• P15_VERTICES returns the vertices for problem 15.
• P16_FUN evaluates the integrand for problem 16.
• P16_TITLE returns the name of problem 16.
• P16_VERTICES returns the vertices for problem 16.
• P17_FUN evaluates the integrand for problem 17.
• P17_TITLE returns the name of problem 17.
• P17_VERTICES returns the vertices for problem 17.
• P18_FUN evaluates the integrand for problem 18.
• P18_TITLE returns the name of problem 18.
• P18_VERTICES returns the vertices for problem 18.
• P19_FUN evaluates the integrand for problem 19.
• P19_TITLE returns the name of problem 19.
• P19_VERTICES returns the vertices for problem 19.
• P20_FUN evaluates the integrand for problem 20.
• P20_TITLE returns the name of problem 20.
• P20_VERTICES returns the vertices for problem 20.
• P21_FUN evaluates the integrand for problem 21.
• P21_TITLE returns the name of problem 21.
• P21_VERTICES returns the vertices for problem 21.
• P22_FUN evaluates the integrand for problem 22.
• P22_TITLE returns the name of problem 22.
• P22_VERTICES returns the vertices for problem 22.
• R8_ABS returns the absolute value of an R8.
• R8_UNIFORM_01 is a unit pseudorandom R8.
• R8MAT_MXM multiplies two matrices.
• R8MAT_MXV multiplies a matrix times a vector.
• R8VEC_COPY copies an R8VEC.
• R8VEC_DOT computes the dot product of a pair of R8VEC's.
• R8VEC_SUM returns the sum of an R8VEC.
• R8VEC_UNIFORM_01 returns a unit pseudorandom R8VEC.
• RANDOM_INITIALIZE initializes the RANDOM random number generator.
• SUBTRIANGLE_NEXT computes the next subtriangle of a triangle.
• TIMESTAMP prints the current YMDHMS date as a time stamp.
• TRIANGLE_AREA computes the area of a triangle.
• TRIANGLE_SAMPLE returns random points in a triangle.

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