sphere_cvt, a FORTRAN90 code which seeks to determine N well-separated sites on the unit sphere in 3D, using centroidal Voronoi tessellation (CVT) techniques.
The code assumes that good separation will follow automatically if the points are the centroids of their Voronoi regions. Thus, the code actually places N points at random on the sphere, and then applies probabilistic centroidal Voronoi tessellation techniques in an attempt to force the the CVT condition to be satisfied. The output of the program is an XYZ file containing the coordinates of the points.
According to Steven Fortune, it is possible to compute the Delaunay triangulation of points on a sphere by computing their convex hull. If the sphere is the unit sphere at the origin, the facet normals are the Voronoi vertices.
sphere_cvt is available in a FORTRAN90 version and a MATLAB version.
SPHERE_DELAUNAY, a FORTRAN90 code which computes the Delaunay triangulation of points on a sphere.
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, among other uses.
SPHERE_INTEGRALS, a FORTRAN90 code which defines test functions for integration over the surface of the unit sphere in 3D.
SPHERE_LEBEDEV_RULE, a dataset directory which contains sets of points on a sphere which can be used for quadrature rules of a known precision;
SPHERE_QUAD, a FORTRAN90 code which estimates the integral of a function defined on the sphere.
SPHERE_STEREOGRAPH, a FORTRAN90 code which computes the stereographic mapping between points on the unit sphere and points on the plane Z = 1; a generalized mapping is also available.
SPHERE_VORONOI, a FORTRAN90 code which computes and plots the Voronoi diagram of points on the unit sphere.
STRIPACK, a FORTRAN90 code which can determine the Voronoi diagram or Delaunay triangulation of a given set of points on the sphere.