urng 0.4.6

Universal Random Number Generator
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178

use colored::Colorize;
use std::time::Instant;
use thousands::Separable;
#[rustfmt::skip]
#[allow(unused_imports)]
use urng::{
    cabi32::{
        mt19937_free, mt19937_new, mt19937_next_u32s, 
        pcg32_free, pcg32_new, pcg32_next_u32s,
        pcg32simd_free, pcg32simd_new, pcg32simd_next_u32s, 
        pcg32x8_free, pcg32x8_new, pcg32x8_next_u32s, 
        philox32_free, philox32_new, philox32_next_u32s, 
        philox32x4_free, philox32x4_new, philox32x4_next_u32s, 
        philox32x4x4_free, philox32x4x4_new, philox32x4x4_next_u32s, 
        sfmt607_free, sfmt607_new, sfmt607_next_u32s,
        sfmt1279_free, sfmt1279_new, sfmt1279_next_u32s,
        sfmt2281_free, sfmt2281_new, sfmt2281_next_u32s,
        sfmt4253_free, sfmt4253_new, sfmt4253_next_u32s,
        sfmt11213_free, sfmt11213_new, sfmt11213_next_u32s,
        sfmt19937_free, sfmt19937_new, sfmt19937_next_u32s,
        sfmt44497_free, sfmt44497_new, sfmt44497_next_u32s,
        sfmt86243_free, sfmt86243_new, sfmt86243_next_u32s,
        sfmt132049_free, sfmt132049_new, sfmt132049_next_u32s,
        sfmt216091_free, sfmt216091_new, sfmt216091_next_u32s,
        splitmix32_free, splitmix32_new, splitmix32_next_u32s, 
        splitmix32simd_free, splitmix32simd_new, splitmix32simd_next_u32s, 
        splitmix32x16_free, splitmix32x16_new, splitmix32x16_next_u32s, 
        squares32_free, squares32_new, squares32_next_u32s,
        squares32simd_free, squares32simd_new, squares32simd_next_u32s, 
        squares32x8_free, squares32x8_new, squares32x8_next_u32s, 
        threefry32x2_free, threefry32x2_new, threefry32x2_next_u32s, 
        threefry32x4_free, threefry32x4_new, threefry32x4_next_u32s,
        xoshiro128pp_free, xoshiro128pp_new, xoshiro128pp_next_u32s,
        xoshiro128ppx16_free, xoshiro128ppx16_new, xoshiro128ppx16_next_u32s,
        xoshiro128ss_free, xoshiro128ss_new, xoshiro128ss_next_u32s,
        xoshiro128ssx16_free, xoshiro128ssx16_new, xoshiro128ssx16_next_u32s,
        jsf32_free, jsf32_new, jsf32_next_u32s,
        jsf32x16_free, jsf32x16_new, jsf32x16_next_u32s,
    },
    cabi64::{
        mt1993764_free, mt1993764_new, mt1993764_next_u64s, 
        philox64_free, philox64_new, philox64_next_u64s, 
        sfmt1993764_free, sfmt1993764_new, sfmt1993764_next_u64s,
        splitmix64_free, splitmix64_new, splitmix64_next_u64s, 
        cet64_free, cet64_new, cet64_next_u64s,
        cet64x8_free, cet64x8_new, cet64x8_next_u64s,
        cet256_free, cet256_new, cet256_next_u64s,
        cet256x2_free, cet256x2_new, cet256x2_next_u64s,
        threefish256_free, threefish256_new, threefish256_next_u64s, 
        xoshiro256pp_free, xoshiro256pp_new, xoshiro256pp_next_u64s,
        xoshiro256ss_free, xoshiro256ss_new, xoshiro256ss_next_u64s, 
        xoshiro256ssx2_free, xoshiro256ssx2_new, xoshiro256ssx2_next_u64s,
        sfc64_free, sfc64_new, sfc64_next_u64s,
        sfc64x8_free, sfc64x8_new, sfc64x8_next_u64s,
        biski64_free, biski64_new, biski64_next_u64s,
        biski64x8_free, biski64x8_new, biski64x8_next_u64s,
    },
};

// philox32x4x4    : 6.80 GS/s
// philox32x4      : 5.40 GS/s
// philox32        : 6.80 GS/s
// threefry32x4    : 5.28 GS/s
// threefry32x2    : 4.62 GS/s
// squares32       : 5.59 GS/s
// squares32x8     : 6.70 GS/s
// squares32simd   : 6.72 GS/s
// pcg32           : 1.06 GS/s
// pcg32x8         : 6.82 GS/s
// pcg32simd       : 6.84 GS/s
// splitmix32      : 2.33 GS/s
// splitmix32x16   : 6.80 GS/s
// splitmix32simd  : 6.68 GS/s
// mt19937         : 0.84 GS/s
// sfmt607         : 1.43 GS/s
// sfmt1279        : 1.57 GS/s
// sfmt2281        : 1.68 GS/s
// sfmt4253        : 1.32 GS/s
// sfmt11213       : 1.45 GS/s
// sfmt19937       : 1.39 GS/s
// sfmt44497       : 1.29 GS/s
// sfmt86243       : 1.17 GS/s
// sfmt132049      : 1.20 GS/s
// sfmt216091      : 1.19 GS/s
// xoshiro128pp    : 0.81 GS/s
// xoshiro128ppx16 : 6.84 GS/s
// xoshiro128ssx16 : 6.88 GS/s
// philox64        : 3.01 GS/s
// splitmix64      : 3.17 GS/s
// mt1993764       : 0.58 GS/s
// sfmt1993764     : 0.53 GS/s
// xoshiro256pp    : 3.13 GS/s
// xoshiro256ss    : 3.14 GS/s
// xoshiro256ssx2  : 3.40 GS/s

const N: usize = 100_000_000;
const G: f64 = 1_000_000_000f64;

#[allow(unused_macros)]
macro_rules! bench32 {
    ($name:ident) => {
        paste::paste! {
            let name = stringify!($name);
            let mut buffer = vec![0u32; N];
            let ptr = [<$name _new>](0);

            let start = Instant::now();
            [<$name _next_u32s>](ptr, buffer.as_mut_ptr(), N);
            let duration = start.elapsed();
            let thruput = N as f64 / duration.as_secs_f64() / G;
            println!(
                "{:<16}: {} {}",
                name.bright_green(),
                format!("{:.2}", thruput).bright_cyan().bold(),
                "GS/s".bright_black(),
            );
            [<$name _free>](ptr);
        }
    };
    ($($name:ident),*) => {
        $(bench32!($name);)*
    };
}

#[allow(unused_macros)]
macro_rules! bench64 {
    ($name:ident) => {
        paste::paste! {
            let name = stringify!($name);
            let mut buffer = vec![0u64; N];
            let ptr = [<$name _new>](0);

            let start = Instant::now();
            [<$name _next_u64s>](ptr, buffer.as_mut_ptr(), N);
            let duration = start.elapsed();
            let thruput = N as f64 / duration.as_secs_f64() / G;
            println!(
                "{:<16}: {} {}",
                name.bright_green(),
                format!("{:.2}", thruput).bright_cyan().bold(),
                "GS/s".bright_black(),
            );
            [<$name _free>](ptr);
        }
    };
    ($($name:ident),*) => {
        $(bench64!($name);)*
    };
}

fn main() {
    println!(
        "Benchmarking Random Number Generators (N = {})",
        N.separate_with_commas().bright_cyan().bold()
    );

    bench32!(philox32x4x4, philox32x4, philox32);
    bench32!(threefry32x4, threefry32x2);
    bench32!(squares32, squares32x8, squares32simd);
    bench32!(pcg32, pcg32x8, pcg32simd);
    bench32!(splitmix32, splitmix32x16, splitmix32simd);
    bench32!(mt19937, sfmt19937);
    bench32!(
        sfmt607, sfmt1279, sfmt2281, sfmt4253, sfmt11213, sfmt44497, sfmt86243, sfmt132049,
        sfmt216091
    );
    bench32!(xoshiro128pp, xoshiro128ppx16, xoshiro128ssx16);
    bench32!(jsf32, jsf32x16);

    bench64!(philox64);
    bench64!(splitmix64);
    bench64!(cet64, cet64x8, cet256, cet256x2);
    bench64!(mt1993764, sfmt1993764);
    bench64!(threefish256);
    bench64!(xoshiro256pp, xoshiro256ss, xoshiro256ssx2);
    bench64!(sfc64, sfc64x8);
    bench64!(biski64, biski64x8);
}