function sample( out, buffer, seed ) {
	var rand;
	var inds;
	var obuf;
	var buf;
	var sb1;
	var sb2;
	var so2;
	var oo;
	var s;
	var M;
	var N;
	var i;

	M = out.shape[ 0 ];
	N = out.shape[ 1 ];

	obuf = out.data;
	so2 = out.strides[ 1 ];
	oo = out.offset;

	buf = buffer.data;
	sb1 = buffer.strides[ 0 ];
	sb2 = buffer.strides[ 1 ];

	// Generate an array of data vector indices...
	inds = [];
	for ( i = 0; i < buffer.shape[ 0 ]; i++ ) {
		inds.push( i );
	}
	// Only randomly sample from the data buffer if the number of centroids is not equal to the number of data vectors...
	if ( M === inds.length ) {
		// Buffer already qualifies as a "random" sample:
		s = inds;
	} else {
		// Create a seeded random sampler (without replacement):
		rand = factory({
			'seed': seed,
			'size': M,
			'mutate': false,
			'replacement': false
		});

		// Generate a random sample:
		s = rand( inds );
	}
	// Update the centroids...
	for ( i = 0; i < M; i++ ) {
		// Note: the following is likely to be an "out-of-order" copy...
		dcopy( N, buf, sb2, sb1*s[i], obuf, so2, oo );
	}
	return out;
}
Initializes centroids by randomly sampling from a data buffer.

@private
@param {ndarray} out - output centroids `kxd` matrix
@param {ndarray} buffer - buffer from which to sample
@param {*} seed - PRNG seed
@returns {ndarray} centroids