randUtils.h

Go to the documentation of this file.

#ifndef RAND_UTILS
#define RAND_UTILS

#include "time.h"

/* FORTRAN functions */

extern "C"
void normalen_(int* n, double* numbers);

extern "C"
void zufalli_(int*);

extern "C"
void zufall_(int* n, double* numbers);

/* C++ wrappers */

inline void normalRand(int n, double* vals)
{
  int tm=time(NULL);
  zufalli_(&tm);
  normalen_(&n, vals);
}
inline void normalRand(double mean, double std, int n, double* vals)
{
  int tm=time(NULL);
  zufalli_(&tm);
  normalen_(&n, vals);

  for(int i=0; i<n; i++)
    vals[i] = mean+std*vals[i];
}
#endif