qdp_random.h

Go to the documentation of this file.

// -*- C++ -*-
//
// QDP data parallel interface
//
// Random number support
 
#ifndef QDP_RANDOM_H
#define QDP_RANDOM_H
 
namespace QDP {
 


 
namespace RNG
{

  void initDefaultRNG(void);

  void initRNG(void);


  void setrn(const Seed& lseed);


  void savern(Seed& lseed);
 
  void finalizeRNG();

  float sranf(Seed& seed, Seed&, const Seed&);

  void sranf(float* d, int N, Seed& seed, ILatticeSeed&, const Seed&);
}

template<class T1, class T2>
inline void
fill_random(float& d, T1& seed, T2& skewed_seed, const T1& seed_mult)
{
  d = float(RNG::sranf(seed, skewed_seed, seed_mult));
}

template<class T1, class T2>
inline void
fill_random(double& d, T1& seed, T2& skewed_seed, const T1& seed_mult)
{
  d = double(RNG::sranf(seed, skewed_seed, seed_mult));
}
 

template<class T1, class T2, int N>
inline void
fill_random(float* d, T1& seed, T2& skewed_seed, const T1& seed_mult)
{
  RNG::sranf(d, N, seed, skewed_seed, seed_mult);
}
 

template<class T1, class T2, int N>
inline void
fill_random(double* d, T1& seed, T2& skewed_seed, const T1& seed_mult)
{
  float dd[N];
  RNG::sranf(dd, N, seed, skewed_seed, seed_mult);
  for(int i=0; i < N; ++i)
    d[i] = float(dd[i]);
}
 
} // namespace QDP
 
#endif