MT19937

class numpy::MT19937

Random number generation class.

Example

#include <numpy/np_ndarray.h>
using namespace numpy;

// Use MT19937
MT19937 obj;
// ... operations ...

Constructors

Signature	Location	Example
MT19937(uint64_tseed = 0)	NP_BITGEN.H:181

Random

Signature	Return Type	Location	Example
void seed(uint64_ts)	void	NP_BITGEN.H:194	View

Other Methods

Signature	Return Type	Location	Example
std::string name()	std::string	NP_BITGEN.H:259	View
double next_double()	double	NP_BITGEN.H:189	View
uint64_t next_uint64()	uint64_t	NP_BITGEN.H:185	View
void set_state(const std::vector<uint64_t>&state)	void	NP_BITGEN.H:225	View

Code Examples

The following examples are extracted from the test suite.

seed (np_test_1_all.cpp:14339)

    }
    // std::cout << "[OK] Generator integers method\n";

    std::cout << " -> tests passed\n";
  }

  void test_vsl_reproducibility() {
    std::cout << "========= test_vsl_reproducibility ====";

    // Test seed reproducibility
    seed(42);
    auto arr1 = normal<double>(0.0, 1.0, { {10} });

    seed(42);
    auto arr2 = normal<double>(0.0, 1.0, { {10} });



    for (size_t i = 0; i < 10; ++i) {
      double val1 = arr1.getElementAt({ i });
      double val2 = arr2.getElementAt({ i });

name (np_test_1_all.cpp:24865)

      uint64_t before_restore = rng.next_uint64();
      rng.set_state(state);
      uint64_t after_restore = rng.next_uint64();

      if (before_restore != after_restore) {
        std::cout << "  [FAIL] : in np_test_bitgen_mt19937() : state restore failed";
        throw std::runtime_error("np_test_bitgen_mt19937 failed: state");
      }

      // Test name
      if (rng.name() != "MT19937") {
        std::cout << "  [FAIL] : in np_test_bitgen_mt19937() : incorrect name";
        throw std::runtime_error("np_test_bitgen_mt19937 failed: name");
      }

      std::cout << " -> tests passed" << std::endl;
    }

    // ============================================================================
    // PCG64 TESTS
    // ============================================================================

next_double (np_test_1_all.cpp:24837)

    // ============================================================================

    void np_test_bitgen_mt19937() {
      std::cout << "========= BitGenerator: MT19937 =======================";

      // Test construction
      numpy::random::MT19937 rng(12345);

      // Test generation
      uint64_t val1 = rng.next_uint64();
      double dval = rng.next_double();

      if (dval < 0.0 || dval >= 1.0) {
        std::cout << "  [FAIL] : in np_test_bitgen_mt19937() : next_double() out of range [" << dval << "]";
        throw std::runtime_error("np_test_bitgen_mt19937 failed: double range");
      }

      // Test seeding reproducibility
      numpy::random::MT19937 rng2(12345);
      uint64_t val2 = rng2.next_uint64();

next_uint64 (np_test_1_all.cpp:24836)

    // MT19937 TESTS
    // ============================================================================

    void np_test_bitgen_mt19937() {
      std::cout << "========= BitGenerator: MT19937 =======================";

      // Test construction
      numpy::random::MT19937 rng(12345);

      // Test generation
      uint64_t val1 = rng.next_uint64();
      double dval = rng.next_double();

      if (dval < 0.0 || dval >= 1.0) {
        std::cout << "  [FAIL] : in np_test_bitgen_mt19937() : next_double() out of range [" << dval << "]";
        throw std::runtime_error("np_test_bitgen_mt19937 failed: double range");
      }

      // Test seeding reproducibility
      numpy::random::MT19937 rng2(12345);
      uint64_t val2 = rng2.next_uint64();

set_state (np_test_1_all.cpp:22363)

    }
    // std::cout << "[OK] Random bytes\n";

    // Test state functions
    std::string state = get_state();
    if (!(!state.empty())) {
        std::string description = std::string("testRandomUtilityFunctions():") + __FILE__ + ":" + std::to_string(__LINE__) + ": !(!state.empty())";
        std::cout << std::string("[FAIL] ") + description << std::endl;
        throw std::runtime_error(description);
    }
    set_state(12345);
    // std::cout << "[OK] State functions\n";

    std::cout << " -> tests passed\n";
  }

  void testGeneratorAPI() {
    std::cout << "========= testGeneratorAPI ====";

    // Test modern Generator API
    auto gen = create_generator(42);