Estimate Integrals over Spherical Triangles

SPHERE_TRIANGLE_QUAD is a FORTRAN90 library which estimates the integral of a scalar function F(X,Y,Z) over a spherical triangle on the unit sphere.

Three methods of estimation are very crude and cannot be improved:

One method of estimation uses random sampling, the Monte Carlo method, whose accuracy can be improved by increasing the number of sample points.

Another method is based on the centroid rule, but allows the user to decompose the original spherical triangle into collection of smaller triangles. The accuracy of the estimate should improve as the size of these triangles decreases.


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


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

Related Data and Programs:

SPHERE_QUAD, a FORTRAN90 library which approximates an integral over the surface of the unit sphere by applying a triangulation to the surface;

SPHERE_TRIANGLE_MONTE_CARLO, a FORTRAN90 library which estimates the integral of a function over a spherical triangle using the Monte Carlo method.

STRIPACK, a FORTRAN90 library which computes the Voronoi diagram or Delaunay triangulation of pointsets on a sphere.

STROUD, a FORTRAN90 library which approximates the integral of a function on the surface or in the interior of a variety of geometric shapes.


Source Code:

Examples and Tests:

List of Routines:

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

Last revised on 23 April 2014.