randomNumbers.h

Go to the documentation of this file.

#ifndef RANDOMNUMBERS_INCLUDED
#define RANDOMNUMBERS_INCLUDED

#include <cmath>
#include <string>

#include "field.h"

#include "bridgeIO.h"


class RandomNumbers
{
 public:
  static const std::string class_name;

 protected:
  Bridge::VerboseLevel m_vl;

 public:

  RandomNumbers()
    : m_vl(CommonParameters::Vlevel()) {}

  virtual ~RandomNumbers() {}

 private:
  // non-copyable
  RandomNumbers(const RandomNumbers&);
  RandomNumbers& operator=(const RandomNumbers&);

 public:

  void set_parameter_verboselevel(const Bridge::VerboseLevel vl) { m_vl = vl; }

  virtual double get() = 0;

  void gauss(double& rand1, double& rand2);


  virtual void gauss_lex_global(Field&);

  virtual void gauss_eo_global(Field&);

  virtual void uniform_lex_global(Field&);

 protected:

  class rand_gauss_even
  {
   public:
    rand_gauss_even(Field& f, RandomNumbers *rng)
      : m_rng(rng), m_ptr(f.ptr(0)), m_block(f.nin()) {}
    void operator()(const bool do_fill);
    size_t block_size() const;

   private:
    RandomNumbers *m_rng;
    double        *m_ptr;
    size_t        m_block;
  };

  class rand_gauss_odd
  {
   public:
    rand_gauss_odd(Field& f, RandomNumbers *rng)
      : m_rng(rng), m_ptr(f.ptr(0)), m_block(f.nin()) {}
    void operator()(const bool do_fill);
    size_t block_size() const;

   private:
    RandomNumbers *m_rng;
    double        *m_ptr;
    size_t        m_block;
  };

  class rand_uniform
  {
   public:
    rand_uniform(Field& f, RandomNumbers *rng)
      : m_rng(rng), m_ptr(f.ptr(0)), m_block(f.nin()) {}
    void operator()(const bool do_fill);
    size_t block_size() const;

   private:
    RandomNumbers *m_rng;
    double        *m_ptr;
    size_t        m_block;
  };

 private:
  template<typename InnerGenerator>
  void generate_global(Field& f);
};
#endif