random.hpp

Go to the documentation of this file.

/*
  RANDOM.hpp  -  support for random number generation with controlled seed

  Copyright (C)         Lumiera.org
    2024,               Hermann Vosseler <Ichthyostega@web.de>

  This program is free software; you can redistribute it and/or
  modify it under the terms of the GNU General Public License as
  published by the Free Software Foundation; either version 2 of
  the License, or (at your option) any later version.

  This program is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with this program; if not, write to the Free Software
  Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.

*/


#ifndef LIB_RANDOM_H
#define LIB_RANDOM_H


#include "lib/integral.hpp"
#include "lib/nocopy.hpp"

#include <random>


namespace lib {
  
  class SeedNucleus
    {
    public:
      virtual ~SeedNucleus();     
      virtual uint64_t getSeed()  =0;
    };
  
  template<class GEN>
  class RandomSequencer
    : util::NonCopyable
    {
      std::uniform_int_distribution<int>      uniformI_;
      std::uniform_int_distribution<uint64_t> uniformU_;
      std::uniform_real_distribution<double>  uniformD_;
      
      GEN generator_;
      
    public:
      RandomSequencer(SeedNucleus&);
      
      int      i32() { return uniformI_(generator_); }
      uint64_t u64() { return uniformU_(generator_); }
      double   uni() { return uniformD_(generator_); }
      
      void randomise(SeedNucleus&);
    };
  
  using Random = RandomSequencer<std::mt19937_64>;
  
  
  extern Random defaultGen;
  
  extern Random entropyGen;
  
  
  /* ===== convenience accessors ===== */
  
  inline int      rani() { return defaultGen.i32(); }
  inline uint64_t ranu() { return defaultGen.u64(); }
  inline double   runi() { return defaultGen.uni(); }
  
  
  void randomiseRandomness();
  
  
  /* ===== Implementation details ===== */
  
  template<class GEN>
  inline
  RandomSequencer<GEN>::RandomSequencer (SeedNucleus& nucleus)
    : uniformI_{0}
    , uniformU_{0}
    , uniformD_{}
    , generator_{nucleus.getSeed()}
    { }
  
  
  template<class GEN>
  inline void
  RandomSequencer<GEN>::randomise (SeedNucleus& nucleus)
  {
    generator_.discard (nucleus.getSeed() % 55555);
  }
  
  
} // namespace lib
#endif /*LIB_RANDOM_H*/