TEST_NEAREST is a FORTRAN90 program which tests several functions that carry out the nearest neighbor calculation.
In a nearest neighbor calculation, we are given:
Obviously, one method to determine the values in NEAREST is simply to compute every distance and take the index of the minimum. But even this simple idea can be implemented in many ways in MATLAB, and implementations will vary in their cost in memory and time.
Also, note that if the dimension M is small, and if the size of the R set is small relative to that of S, it may be much cheaper to compute the Delaunay triangulation of R (or its higher-dimensional generalization). Computing the triangulation is somewhat expensive, but makes the search procedure extremely quick.
Lloyd's version of the Centroidal Voronoi Tessellation (CVT) algorithm estimates the area or volume of the cells using sampling. This computation is a nearest neighbor calculation.
The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.
TEST_NEAREST is available in a C version and a C++ version and a FORTRAN77 version and a FORTRAN90 version and a MATLAB version.
CVT, a FORTRAN90 library which computes elements of a Centroidal Voronoi Tessellation (CVT).
NEAREST_INTERP_1D, a FORTRAN90 library which interpolates a set of data using a piecewise constant interpolant defined by the nearest neighbor criterion.
You can go up one level to the FORTRAN90 source codes.