sidef/scripts/RosettaCode/MRG32k3a.sf

#!/usr/bin/ruby

# https://rosettacode.org/wiki/Pseudo-random_numbers/Combined_recursive_generator_MRG32k3a

class MRG32k3a(seed) {

    define(
        m1 = (2**32 - 209)
        m2 = (2**32 - 22853)
    )

    define(
        a1 = %n<     0 1403580  -810728>
        a2 = %n<527612       0 -1370589>
    )

    has x1 = [seed, 0, 0]
    has x2 = x1.clone

    method next_int {
        x1.unshift(a1.map_kv {|k,v| v * x1[k] }.sum % m1); x1.pop
        x2.unshift(a2.map_kv {|k,v| v * x2[k] }.sum % m2); x2.pop
        (x1[0] - x2[0]) % (m1 + 1)
    }

    method next_float { self.next_int / (m1 + 1) -> float }
}

say "Seed: 1234567, first 5 values:"
var rng = MRG32k3a(seed: 1234567)
var arr = 5.of{ rng.next_int }
assert_eq(arr,[1459213977,2827710106,4245671317,3877608661,2595287583])
say arr