function nearest = find_closest ( ndim, n, sample_num, s, r )

%% FIND_CLOSEST finds the nearest R point to each S point.
%
%  Discussion:
%
%    This routine finds the closest Voronoi cell generator by checking every
%    one.  For problems with many cells, this process can take the bulk
%    of the CPU time.  Other approaches, which group the cell generators into
%    bins, can run faster by a large factor.
%
%  Modified:
%
%    22 October 2004
%
%  Author:
%
%    John Burkardt
%
%  Parameters:
%
%    Input, integer NDIM, the spatial dimension.
%
%    Input, integer N, the number of cell generators.
%
%    Input, integer SAMPLE_NUM, the number of sample points.
%
%    Input, real S(NDIM,SAMPLE_NUM), the points to be checked.
%
%    Input, real R(NDIM,N), the cell generators.
%
%    Output, integer NEAREST(SAMPLE_NUM), the index of the nearest cell generators.
%
  for js = 1 : sample_num

    distance = Inf;
    nearest(js) = -1;

    for jr = 1 : n

      dist_sq = 0.0E+00;
      for i = 1 : ndim
        dist_sq = dist_sq + ( r(i,jr) - s(i,js) )^2;
      end

      if ( dist_sq < distance )
        distance = dist_sq;
        nearest(js) = jr;
      end

    end

  end