ytRNG.h

/*
  math/ytRNG.{h,c} : Random Number Generator
  Copyright (C) 2018, Yoshinori Tamada <tamada  A T  ytlab.jp>
  All rights reserved.
 
  See LICENSE.txt for details of the licensing agreement.
*/
 
#ifndef __YTLIB_RNG_H
#define __YTLIB_RNG_H
 
#if defined(USE_FUJITSU_FCC_SSL)
 
#elif !defined(USE_DSFMT) && defined(USE_INTEL_MKL)
#include <mkl_vsl.h>
#else
#include "ytRNG_dSFMT.h"
#endif
 
typedef struct {
#if defined(USE_FUJITSU_FCC_SSL)
  int ssl_dvrau4_seed;
  double * ssl_dvrau4_work;
  int ssl_dvrau4_work_size;
 
  int ssl_dvran3_seed;
  double * ssl_dvran3_work;
  int ssl_dvran3_work_size;
#elif !defined(USE_DSFMT) && defined(USE_INTEL_MKL)
  VSLStreamStatePtr mkl_vsl_stream;
#else
  ytRNG_dSFMT * dsfmt;
  double work[DSFMT_N64];
#endif
} ytRNG;
 
ytRNG *   ytRNG_new();
void      ytRNG_delete(ytRNG * rng);
void      ytRNG_setSeed(ytRNG * rng, int seed);
void      ytRNG_normal(ytRNG * rng, double * ar, int n, double mean, double sigma);
void      ytRNG_uniform(ytRNG * rng, double * ar, int n, double min, double max);
void      ytRNG_integer(ytRNG * rng, int * ar, int n, int min, int max);
void      ytRNG_perm(ytRNG * rng, int * ar, int n);
void      ytRNG_shuffle(ytRNG * rng, int * ar, int n);
int       ytRNG_debug(int argc, char * argv[]);
 
#endif /* __YTLIB_RNG_H */