msci Code

//provides a simple and uniform interface for random number generation...
#include <sys/timeb.h>

#include "peteys_tmpl_lib.h"
#include "randomize.h"

namespace libpetey {

  long libpetey_ran_seed;
  gsl_rng *libpetey_ran;

  long seed_from_clock() {
    long seed;
    timeb now;
    ftime(&now);
    seed= (long) now.time + ((long) now.millitm) << 7;
    //printf("seed = %d\n", seed);
    return seed;
  }

  void ran_init(const gsl_rng_type *rt) {
    libpetey_ran=gsl_rng_alloc(rt);
    gsl_rng_set(libpetey_ran, seed_from_clock());
  }

  long *randomize(long n) {
    double *flt;
    long *ind;

    flt=new double[n];
    for (long i=0; i<n; i++) {
      flt[i]=ranu();
      printf("%g\n", flt[i]);
    }
    ind=heapsort(flt, n);

    delete [] flt;
    return ind;
  }

  //not really necessary, since one should only call "ran_init"
  //once in a given program...
  void ran_end() {
    gsl_rng_free(libpetey_ran);
  }

} //end namespace libpetey