Pseudo-random numbers/Splitmix64: Difference between revisions

Line 851:

Seed = 987654321:

0: 20027, 1: 19892, 2: 20073, 3: 19978, 4: 20030</pre>

=={{header|Perl}}==

<syntaxhighlight lang="perl">use strict;

use warnings;

no warnings 'portable';

use feature 'say';

use Math::AnyNum qw(:overload);

package splitmix64 {

sub new {

my ($class, %opt) = @_;

bless {state => $opt{seed}}, $class;

}

sub next_int {

my ($self) = @_;

my $next = $self->{state} = ($self->{state} + 0x9e3779b97f4a7c15) & (2**64 - 1);

$next = ($next ^ ($next >> 30)) * 0xbf58476d1ce4e5b9 & (2**64 - 1);

$next = ($next ^ ($next >> 27)) * 0x94d049bb133111eb & (2**64 - 1);

($next ^ ($next >> 31)) & (2**64 - 1);

}

sub next_float {

my ($self) = @_;

$self->next_int / 2**64;

}

say 'Seed: 1234567, first 5 values:';

my $rng = splitmix64->new(seed => 1234567);

say $rng->next_int for 1 .. 5;

my %h;

say "\nSeed: 987654321, values histogram:";

$rng = splitmix64->new(seed => 987654321);

$h{int 5 * $rng->next_float}++ for 1 .. 100_000;

say "$_ $h{$_}" for sort keys %h;</syntaxhighlight>

<pre>Seed: 1234567, first 5 values:

6457827717110365317

3203168211198807973

9817491932198370423

4593380528125082431

16408922859458223821

Seed: 987654321, values histogram:

0 20027

1 19892

2 20073

3 19978

4 20030</pre>

=={{header|Phix}}==

As per [[Pseudo-random_numbers/PCG32#Phix]], resorting to mpfr/gmp

with javascript_semantics

include mpfr.e

mpz state = mpz_init(),

shift = mpz_init("0x9e3779b97f4a7c15"),

mult1 = mpz_init("0xbf58476d1ce4e5b9"),

mult2 = mpz_init("0x94d049bb133111eb"),

b64 = mpz_init("0x10000000000000000"), -- (truncate to 64 bits)

tmp = mpz_init(),

z = mpz_init()

procedure seed(integer num)

mpz_set_si(state,num)

end procedure

procedure next_int()

mpz_add(state, state, shift) -- state += shift

mpz_fdiv_r(state, state, b64) -- state := remainder(z,b64)

mpz_set(z, state) -- z := state

mpz_tdiv_q_2exp(tmp, z, 30) -- tmp := trunc(z/2^30)

mpz_xor(z, z, tmp) -- z := xor_bits(z,tmp)

mpz_mul(z, z, mult1) -- z *= mult1

mpz_fdiv_r(z, z, b64) -- z := remainder(z,b64)

mpz_tdiv_q_2exp(tmp, z, 27) -- tmp := trunc(z/2^27)

mpz_xor(z, z, tmp) -- z := xor_bits(z,tmp)

mpz_mul(z, z, mult2) -- z *= mult2

mpz_fdiv_r(z, z, b64) -- z := remainder(z,b64)

mpz_tdiv_q_2exp(tmp, z, 31) -- tmp := trunc(z/2^31)

mpz_xor(z, z, tmp) -- z := xor_bits(z,tmp)

-- (result left in z)

end procedure

function next_float()

next_int()

mpfr f = mpfr_init_set_z(z)

mpfr_div_z(f, f, b64)

return mpfr_get_d(f)

end function

seed(1234567)

for i=1 to 5 do

next_int()

printf(1,"%s\n",mpz_get_str(z))

end for

seed(987654321)

sequence r = repeat(0,5)

for i=1 to 100000 do

integer rdx = floor(next_float()*5)+1

r[rdx] += 1

end for

?r

<pre>

6457827717110365317

3203168211198807973

9817491932198370423

4593380528125082431

16408922859458223821

{20027,19892,20073,19978,20030}

</pre>

=={{header|Object Pascal}}==

Line 1,109:

Line 993:

16408922859458223821

0: 20027 1: 19892 2: 20073 3: 19978 4: 20030

</pre>

=={{header|Perl}}==

<syntaxhighlight lang="perl">use strict;

use warnings;

no warnings 'portable';

use feature 'say';

use Math::AnyNum qw(:overload);

package splitmix64 {

sub new {

my ($class, %opt) = @_;

bless {state => $opt{seed}}, $class;

}

sub next_int {

my ($self) = @_;

my $next = $self->{state} = ($self->{state} + 0x9e3779b97f4a7c15) & (2**64 - 1);

$next = ($next ^ ($next >> 30)) * 0xbf58476d1ce4e5b9 & (2**64 - 1);

$next = ($next ^ ($next >> 27)) * 0x94d049bb133111eb & (2**64 - 1);

($next ^ ($next >> 31)) & (2**64 - 1);

}

sub next_float {

my ($self) = @_;

$self->next_int / 2**64;

}

say 'Seed: 1234567, first 5 values:';

my $rng = splitmix64->new(seed => 1234567);

say $rng->next_int for 1 .. 5;

my %h;

say "\nSeed: 987654321, values histogram:";

$rng = splitmix64->new(seed => 987654321);

$h{int 5 * $rng->next_float}++ for 1 .. 100_000;

say "$_ $h{$_}" for sort keys %h;</syntaxhighlight>

<pre>Seed: 1234567, first 5 values:

6457827717110365317

3203168211198807973

9817491932198370423

4593380528125082431

16408922859458223821

Seed: 987654321, values histogram:

0 20027

1 19892

2 20073

3 19978

4 20030</pre>

=={{header|Phix}}==

As per [[Pseudo-random_numbers/PCG32#Phix]], resorting to mpfr/gmp

with javascript_semantics

include mpfr.e

mpz state = mpz_init(),

shift = mpz_init("0x9e3779b97f4a7c15"),

mult1 = mpz_init("0xbf58476d1ce4e5b9"),

mult2 = mpz_init("0x94d049bb133111eb"),

b64 = mpz_init("0x10000000000000000"), -- (truncate to 64 bits)

tmp = mpz_init(),

z = mpz_init()

procedure seed(integer num)

mpz_set_si(state,num)

end procedure

procedure next_int()

mpz_add(state, state, shift) -- state += shift

mpz_fdiv_r(state, state, b64) -- state := remainder(z,b64)

mpz_set(z, state) -- z := state

mpz_tdiv_q_2exp(tmp, z, 30) -- tmp := trunc(z/2^30)

mpz_xor(z, z, tmp) -- z := xor_bits(z,tmp)

mpz_mul(z, z, mult1) -- z *= mult1

mpz_fdiv_r(z, z, b64) -- z := remainder(z,b64)

mpz_tdiv_q_2exp(tmp, z, 27) -- tmp := trunc(z/2^27)

mpz_xor(z, z, tmp) -- z := xor_bits(z,tmp)

mpz_mul(z, z, mult2) -- z *= mult2

mpz_fdiv_r(z, z, b64) -- z := remainder(z,b64)

mpz_tdiv_q_2exp(tmp, z, 31) -- tmp := trunc(z/2^31)

mpz_xor(z, z, tmp) -- z := xor_bits(z,tmp)

-- (result left in z)

end procedure

function next_float()

next_int()

mpfr f = mpfr_init_set_z(z)

mpfr_div_z(f, f, b64)

return mpfr_get_d(f)

end function

seed(1234567)

for i=1 to 5 do

next_int()

printf(1,"%s\n",mpz_get_str(z))

end for

seed(987654321)

sequence r = repeat(0,5)

for i=1 to 100000 do

integer rdx = floor(next_float()*5)+1

r[rdx] += 1

end for

?r

<pre>

6457827717110365317

3203168211198807973

9817491932198370423

4593380528125082431

16408922859458223821

{20027,19892,20073,19978,20030}

</pre>