function r = cvt_sample ( ndim, n, n_now, sample, initialize )

%*****************************************************************************80
%
%% cvt_sample() returns sample points.
%
%  Discussion:
%
%    N sample points are to be taken from the unit box of dimension NDIM.
%
%    It may not be practical to generate all the sample points in a 
%    single call.  For that reason, the routine allows the user to
%    request a total of N points, but to require that only N_NOW be
%    generated now (on this call).  
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    20 January 2021
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer NDIM, the spatial dimension.
%
%    integer N, the number of sample points to be generated.
%
%    integer N_NOW, the number of sample points to be generated
%    on this call.  N_NOW must be at least 1.
%
%    integer SAMPLE, ignored.
%
%    logical INITIALIZE: true if the pseudorandom process should be
%    reinitialized.
%
%  Output:
%
%    real R(NDIM,N_NOW), the initial value for the generators.
%
  if ( n_now < 1 )
    fprintf ( 1, '\n' );
    fprintf ( 1, 'cvt_sample(): Fatal error!\n' );
    fprintf ( 1, '  N_NOW < 1.\n' );
    error ( 'cvt_sample(): Fatal error!' );
  end

  r = rand ( ndim, n_now );

  return
end

