Program Listing for File RandomDraw.hpp

Return to documentation for file (include/sgpl/operations/unary/RandomDraw.hpp)

#pragma once
#ifndef SGPL_OPERATIONS_UNARY_RANDOMDRAW_HPP_INCLUDE
#define SGPL_OPERATIONS_UNARY_RANDOMDRAW_HPP_INCLUDE

#include <cassert>
#include <map>
#include <set>
#include <string>

#include "../../../../third-party/conduit/include/uit_emp/tools/string_utils.hpp"

#include "../../hardware/Core.hpp"
#include "../../program/Instruction.hpp"
#include "../../program/Program.hpp"
#include "../../utility/ThreadLocalRandom.hpp"

namespace sgpl {

class RandomDraw {

  template<typename Spec>
  static double map_between_plusminus_one(
    const typename Spec::tag_t& tag
  ) noexcept {

    constexpr double max = 1.0;
    constexpr double min = -1.0;
    constexpr double max_double = Spec::tag_t::MaxDouble();

    return  (tag.GetDouble() / max_double) * (max - min) + min;

  }

  template<typename Spec>
  static bool is_odd(const typename Spec::tag_t& tag) noexcept {
    return tag.Get(0);
  }

  static double map_up( const double plusminus_unit_val ) {
    assert( plusminus_unit_val != 0 );

    return 1.0 / (
      plusminus_unit_val * plusminus_unit_val * plusminus_unit_val
    );

  }

  template<typename Spec>
  static float map_draw() noexcept {
    const typename Spec::tag_t tag( sgpl::tlrand.Get() );
    return is_odd<Spec>( tag )
      ? map_up( map_between_plusminus_one<Spec>(tag) )
      : map_between_plusminus_one<Spec>( tag )
    ;
  }

public:

  template<typename Spec>
  static void run(
    sgpl::Core<Spec>& core,
    const sgpl::Instruction<Spec>& inst,
    const sgpl::Program<Spec>&,
    typename Spec::peripheral_t&
  ) noexcept {
    core.registers[ inst.args[0] ] = map_draw<Spec>();
  }

  static std::string name() { return "Draw Random Value"; }

  static size_t prevalence() { return 1; }

  template<typename Spec>
  static auto descriptors( const sgpl::Instruction<Spec>& inst ) {

    return std::map<std::string, std::string>{
      { "argument a", uit_emp::to_string( static_cast<int>( inst.args[0] ) ) },
      { "summary", "a = weighted draw from random number generator" },
    };
  }

  template<typename Spec>
  static std::set<std::string> categories(const sgpl::Instruction<Spec>&) {
    return {
      "calculation",
      "intrinsic",
      "op",
      "stochastic",
    };
  }

};

} // namespace sgpl

#endif // #ifndef SGPL_OPERATIONS_UNARY_RANDOMDRAW_HPP_INCLUDE