RNGManager.hpp

/*
 * (c) Copyright 2020 CORSIKA Project, corsika-project@lists.kit.edu
 *
 * This software is distributed under the terms of the GNU General Public
 * Licence version 3 (GPL Version 3). See file LICENSE for a full version of
 * the license.
 */

#pragma once

#include <map>
#include <cstdint>
#include <random>
#include <string>

#include <corsika/framework/utility/Singleton.hpp>
#include <corsika/framework/core/Logging.hpp>
#include <corsika/detail/framework/random/random_iterator/Stream.hpp>
namespace corsika {

  template <typename CBPRNG = random_iterator::philox>
  class RNGManager : public corsika::Singleton<RNGManager<CBPRNG>> {

    friend class corsika::Singleton<RNGManager<CBPRNG>>;

  public:
    typedef CBPRNG prng_type;
    typedef std::uint64_t seed_type;
    typedef std::string string_type;
    typedef std::map<std::string, prng_type> streams_type;

    RNGManager(RNGManager<CBPRNG> const&) = delete; // since it is a singleton

    RNGManager<CBPRNG>& operator=(RNGManager<CBPRNG> const&) = delete;

    inline void registerRandomStream(string_type const& streamName);

    inline prng_type& getRandomStream(string_type const& streamName);

    inline bool isRegistered(string_type const& streamName) const;

    inline std::stringstream dumpState() const;

    const streams_type& getRngs() const { return rngs_; }

    streams_type& Rngs() { return rngs_; }

    seed_type getSeed() const { return seed_; }

    void setSeed(seed_type seed) {
      seed_ = seed;
      // update the rgn states
      for (auto& [streamName, rng] : rngs_) rng.setSeed(seed_);
    }

  protected:
    RNGManager() = default;

  private:
    streams_type rngs_;
    seed_type seed_;
  };

  typedef typename RNGManager<>::prng_type default_prng_type;

} // namespace corsika

#include <corsika/detail/framework/random/RNGManager.inl>