randomNumbers_Mseries.h

Go to the documentation of this file.

#ifndef RANDOMNUMBERS_MSERIES_INCLUDED
#define RANDOMNUMBERS_MSERIES_INCLUDED

#include <assert.h>
#include <cmath>
#include <fstream>
#include <iostream>
#include <string>

#include "randomNumbers.h"
#include "communicator.h"

#include "bridgeIO.h"
using Bridge::vout;


class RandomNumbers_Mseries : public RandomNumbers
{
 protected:
  Bridge::VerboseLevel m_vl;

 private:
  static const int Np = 521, Nq = 32;
  double           Fnorm;
  int              w[Np];
  int              jr, kr;

  double sq2r;
  double pi, pi2;

 public:
  RandomNumbers_Mseries(const int ndelay)
  {
    Fnorm = 4.656612870908988e-10;
    initset(ndelay);
  }

  RandomNumbers_Mseries(const std::string filename)
  {
    Fnorm = 4.656612870908988e-10;
    set_readfile(filename);
  }

  double get()
  {
    w[jr] = w[jr] ^ w[kr];
    double rw = w[jr] * Fnorm;
    jr = jr + 1;
    if (jr >= Np) jr = jr - Np;
    kr = kr + 1;
    if (kr >= Np) kr = kr - Np;
    return rw;
  }

  void writefile(const std::string);

 private:
  void set_readfile(const std::string);

  void initset(const int ndelay);

  void delay3(const int ndelay);
};
#endif