msci Code

#ifndef RANDOMIZE_H
#define RANDOMIZE_H 1

//provides a simple and uniform interface for random number generation...
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>

extern gsl_rng *libpetey_ran;

long seed_from_clock();

void ran_init(const gsl_rng_type *rt=gsl_rng_mt19937);

//uniform deviates:
inline double ranu() {
  return gsl_rng_uniform(libpetey_ran);
}

//gaussian deviates:
inline double rang() {
  return gsl_ran_ugaussian(libpetey_ran);
}

long *randomize(long n);

void ran_end();

#endif