rdrandをtestU01とPractRandで検定する

結論

TestU01,BigCrush で検定すると、良い乱数でも時々異常と判定されることがある。
PractRand -te 0, -te 1 で検定すると、良い乱数でも時々 unusual と判定されるされることがある。
PractRand -te 10 で検定すると、良い乱数でも時々 suspicious や FAIL と判定されることがある。

rdrandの速さについて

rdrand は第３世代以降のＣＰＵで使用可能です。
rdrand は疑似乱数と比較して非常に遅いです。
第１０世代の Core i5-10500T＋Windows11 23H2 上で実行して、
このページの「testU01で検定する」の結果を得るのに３日半かかりました。
第１１世代のＣＰＵでは rdrand はさらに約４０倍時間がかかります。
つまり、結果を得るのに約４か月半かかることになります。
第１０世代以前のＣＰＵでも Windows11 24H2 以降の場合、C:\Windows\System32にある
mcupdate_GenuineIntel.dll をリネームする必要があります。リネームの方法は、
日記という名の不定記 https://griffonworks.net/~yamagin/diary/archive/6854.php
が参考になります。
ただし、リネームするとＯＳが起動しなくなる可能性があります。
Ubuntu の場合は、
ハードウエア乱数の生成速度 https://isaku-wada.my.coocan.jp/tips/hrng.html
で説明しました。
どちらの場合もＰＣが脆弱になるのでお勧めできません。
自己責任でお願いします。検定が終わったら元に戻しましょう。

rnc.h

gcc, clang, vc2012 .. vc2022 でコンパイル可能です。

define.h はこちら

/* rnc.h : coded by Isaku WADA(isaku@pb4.so-net.ne.jp) 2025
A hardware RNG.
gcc, clang, vc2012 .. vc2022
*/
#if !defined(RNC_H__)
 #define RNC_H__
 #include "define.h"

typedef char Rnc_t;

 #if defined(__clang__)
  #pragma clang attribute push(__attribute__((target("rdrnd"))), \
      apply_to=function)
 #elif defined(__GNUC__)
  #pragma GCC push_options
  #pragma GCC target("rdrnd")
 #endif

static __inline qword RncRand(Rnc_t*p) {
    qword r; (void)p;
 #if defined(__x86_64__) || defined(_M_X64)
    while (_rdrand64_step(&r)==0) _mm_pause();
 #else
    while (_rdrand32_step((unsigned*)&r)==0) _mm_pause();
    while (_rdrand32_step((unsigned*)&r+1)==0) _mm_pause();
 #endif
    return r;
}
 #if defined(__clang__)
  #pragma clang attribute pop
 #elif defined(__GNUC__)
  #pragma GCC pop_options
 #endif

static __inline int RncHaveRdrand(void) {
    int v[4]; CpuidEcx0(v,7);
    return (v[1]>>18)&1;
}

static __inline const char*RncInit(Rnc_t*p,qword s) {
    if (RncHaveRdrand()==0) return"RDRAND";
    (void)p; (void)s; return 0;
}

#endif

#if !defined(ExRand) && !defined(ROC_H__)
 #define ExRand   RncRand
 #define ExName  "RncRand"
 #define ExInit   RncInit
 #define Ex_t     Rnc_t
 #define ExAlign  1
#endif

roc.h

roc.h は rnc.h を OpenMP でマルチスレッド化したプログラムです。

gcc, clang, vc2012 .. vc2022 でコンパイル可能です。

直前の rnc.h と define.h が必要です。

define.h はこちら

/* roc.h : coded by Isaku WADA(isaku@pb4.so-net.ne.jp) 2025
A hardware RNG with OpenMP.
gcc, clang, vc2012 .. vc2022
msys2:clang -> clang -fopenmp foo.c -lwinmm -static -lomp
msys2:gcc   -> gcc   -fopenmp foo.c -lwinmm -static -lgomp
*/
#if !defined(ROC_H__)
 #define ROC_H__
 #if !defined(_OPENMP)
  #error OpenMP
 #endif
 #include <omp.h>
 #include "rnc.h"
 #define ROC_NUM 8192

typedef Rnc_t Roc_t;

static __inline qword RocRand(Roc_t*p) {
    static qword buf[ROC_NUM]; static int idx=ROC_NUM; int i; (void)p;
    if (omp_in_parallel()) return RncRand(p);
    if (idx>=ROC_NUM) {
 #if defined(RO_NUMTH)
        #pragma omp parallel for num_threads(RO_NUMTH)
 #else
        #pragma omp parallel for
 #endif
        for (i=0;i<ROC_NUM;i++) buf[i]=RncRand(p);
        idx=0;
    }
    return buf[idx++];
}

 #define RocInit RncInit

#endif

#if !defined(ExRand)
 #define ExRand   RocRand
 #define ExName  "RocRand"
 #define ExInit   RocInit
 #define Ex_t     Roc_t
 #define ExAlign  1
#endif

rdrandをtestU01で検定する

以下にtestU01で検定した結果の要約を示します。

type                 , c, test                     , p-value 
RocRand rawhi Seed(1), B, ClosePairs mNP2,t=3      , 4.9e-4  
RocRand rawhi Seed(2), B, CollisionOver,t=3        , 3.4e-4  
RocRand rawhi Seed(3), B, Passed                   , -       
RocRand rawhi Seed(4), B, Passed                   , -       
RocRand rawhi Seed(5), B, Passed                   , -       
RocRand rawhi Seed(6), B, Passed                   , -       
RocRand rawlo Seed(1), B, Passed                   , -       
RocRand rawlo Seed(2), B, Passed                   , -       
RocRand rawlo Seed(3), B, Passed                   , -       
RocRand rawlo Seed(4), B, Passed                   , -       
RocRand rawlo Seed(5), B, Passed                   , -       
RocRand rawlo Seed(6), B, LongestHeadRun,r=27      , 0.9995  
RocRand revhi Seed(1), B, Passed                   , -       
RocRand revhi Seed(2), B, Passed                   , -       
RocRand revhi Seed(3), B, Passed                   , -       
RocRand revhi Seed(4), B, CollisionOver,t=3        , 0.9991  
RocRand revhi Seed(5), B, RandomWalk1 R(L=1000,r=0), 2.4e-4  
-                    , B, AutoCor,d=1,r=0          , 3.5e-4  
RocRand revhi Seed(6), B, Passed                   , -       
RocRand revlo Seed(1), B, Passed                   , -       
RocRand revlo Seed(2), B, CollisionOver,t=21       , 8.6e-4  
RocRand revlo Seed(3), B, CouponCollector,r=20     , 6.2e-4  
-                    , B, RandomWalk1 M(L=50,r=25) , 8.5e-4  
RocRand revlo Seed(4), B, Passed                   , -       
RocRand revlo Seed(5), B, Passed                   , -       
RocRand revlo Seed(6), B, LempelZiv,r=15           , 1-9.9e-5
-                    , B, Run of bits,r=27         , 1.0e-3  
RocRand swphi Seed(1), B, Passed                   , -       
RocRand swphi Seed(2), B, Passed                   , -       
RocRand swphi Seed(3), B, Passed                   , -       
RocRand swphi Seed(4), B, CollisionOver,t=14       , 1-6.0e-5
RocRand swphi Seed(5), B, CollisionOver,t=7        , 0.9997  
-                    , B, SampleProd,t=16          , 2.4e-4  
RocRand swphi Seed(6), B, Passed                   , -       
RocRand swplo Seed(1), B, PeriodsInStrings,r=0     , 9.5e-4  
RocRand swplo Seed(2), B, Passed                   , -       
RocRand swplo Seed(3), B, Passed                   , -       
RocRand swplo Seed(4), B, Passed                   , -       
RocRand swplo Seed(5), B, Passed                   , -       
RocRand swplo Seed(6), B, Passed                   , -

rawhi は64ビットの上位32ビットを検定の対象としたものです。
rawlo は64ビットの下位32ビットを検定の対象としたものです。
revhi, revlo はビット毎に順序を反転したものです。
swphi, swplo はバイト毎に順序を反転したものです。
B は BigCrush を意味します。
このように、何度も検定すると何回か異常と判断されます。

rdrandをPractRandで検定する

PractRandで rdrand を４０回 1GB, 2GB, 4GB の３つの長さで、異常のあったものを表示します。
この１２０回分を１つのテストとします。
-te 0 を te0、-te 1 を te1、-te 10 を te10、-tf 1 を tf1、-tf 2 を te2 と表します。
１つのテストを te0, te1, te10 と tf1, tf2 の６個の組み合わせをテストして合計７２０回調べました。

RocRand te0 tf1 Seed(4) 2GB  [Low16/64]BCFN(2+8,13-5,T) p=1-1.2e-5 unusual 
RocRand te0 tf1 Seed(4) 4GB  [Low4/64]mod3n(5):(2,9-6) p=5.3e-5 unusual 
RocRand te0 tf1 Seed(6) 1GB  [Low4/64]DC6-9x1Bytes-1 p=1-3.6e-3 unusual 
RocRand te0 tf1 Seed(6) 2GB  FPF-14+6/16:cross p=1-3.6e-4 unusual 
RocRand te0 tf1 Seed(8) 4GB  [Low4/64]Gap-16:A p=1-6.7e-4 unusual 
RocRand te0 tf1 Seed(10) 4GB  DC6-9x1Bytes-1 p=1-4.0e-3 unusual 
RocRand te0 tf1 Seed(15) 4GB  [Low16/64]DC6-9x1Bytes-1 p=2.1e-3 unusual 
RocRand te0 tf1 Seed(18) 4GB  [Low1/64]FPF-14+6/16:cross p=1-4.1e-4 unusual 
RocRand te0 tf1 Seed(23) 4GB  [Low4/64]mod3n(5):(1,9-6) p=1.0e-4 unusual 
RocRand te0 tf1 Seed(27) 1GB  [Low16/64]FPF-14+6/16:(14,14-8) p=4.6e-6 unusual 
RocRand te0 tf1 Seed(36) 2GB  [Low4/64]DC6-9x1Bytes-1 p=1.6e-3 unusual 
RocRand te0 tf1 Seed(37) 2GB  BCFN(2+0,13-0,T) p=1-9.6e-4 unusual

RocRand te0 tf2 Seed(6) 1GB  [Low4/16]BCFN(2+0,13-2,T) p=5.1e-5 unusual 
RocRand te0 tf2 Seed(9) 1GB  [Low1/16]BCFN(2+3,13-4,T) p=1-1.1e-5 unusual 
RocRand te0 tf2 Seed(9) 2GB  [Low1/32]Gap-16:A p=1-2.3e-4 unusual 
RocRand te0 tf2 Seed(18) 2GB  [Low8/32]BCFN(2+0,13-1,T) p=1-3.3e-4 unusual 
RocRand te0 tf2 Seed(22) 1GB  [Low4/64]DC6-9x1Bytes-1 p=1.1e-3 unusual 
RocRand te0 tf2 Seed(35) 4GB  [Low4/16]FPF-14+6/16:cross p=1.3e-4 unusual 
RocRand te0 tf2 Seed(36) 1GB  [Low1/32]BCFN(2+3,13-5,T) p=1.5e-5 unusual 
RocRand te0 tf2 Seed(36) 1GB  [Low1/32]BCFN(2+5,13-6,T) p=1.9e-5 unusual

RocRand te1 tf1 Seed(7) 4GB  FPF-14+6/32:all p=7.6e-5 unusual 
RocRand te1 tf1 Seed(24) 4GB  [Low16/64]mod3n(0):(3,9-6) p=1.3e-5 unusual 
RocRand te1 tf1 Seed(29) 1GB  [Low16/64]BCFN(0+0,13-1,T) p=1-1.5e-4 unusual 
RocRand te1 tf1 Seed(30) 4GB  [Low4/64]BCFN(0+3,13-2,T) p=1-1.5e-4 unusual 
RocRand te1 tf1 Seed(31) 1GB  [Low16/64]BCFN(0+4,13-3,T) p=1-8.0e-5 unusual 
RocRand te1 tf1 Seed(31) 4GB  [Low1/64]DC7-9x1Bytes-1:both p=1.1e-3 unusual 
RocRand te1 tf1 Seed(37) 1GB  BCFN(0+0,13-0,T) p=1-5.8e-5 unusual

RocRand te1 tf2 Seed(9) 4GB  [Low4/32]DC7-5x4Bytes-1:indep p=1-1e-4 unusual 
RocRand te1 tf2 Seed(14) 1GB  [Low1/16]BCFN(0+4,13-4,T) p=1-6.0e-5 unusual 
RocRand te1 tf2 Seed(20) 4GB  [Low4/16]BCFN(0+3,13-1,T) p=9.4e-5 unusual 
RocRand te1 tf2 Seed(29) 1GB  [Low4/16]DC7-5x4Bytes-1:indep p=1-4e-5 unusual 
RocRand te1 tf2 Seed(34) 4GB  [Low1/64]DC7-9x1Bytes-1:odd p=1.9e-4 unusual

RocRand te10 tf1 Seed(4) 2GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf1 Seed(7) 4GB  BDayS128(26):all1 p=1e-3 unusual 
RocRand te10 tf1 Seed(7) 4GB  BDayS128(26):all2 p=6e-5 suspicious 
RocRand te10 tf1 Seed(12) 2GB  BDayS128(26):all2 p=3e-3 unusual 
RocRand te10 tf1 Seed(13) 4GB  [Low4/64]BDayS128(26):inc2 p=8e-6 suspicious

RocRand te10 tf2 Seed(1) 2GB  [Low4/32]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(1) 4GB  [Low1/16]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(4) 1GB  [Low1/8]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(5) 4GB  [Low1/8]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(7) 1GB  [Low4/16]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(7) 2GB  [Low1/32]BDayS128(26):inc1 p=2e-3 unusual 
RocRand te10 tf2 Seed(7) 2GB  [Low1/32]BDayS128(26):inc2 p=9e-12 FAIL 
RocRand te10 tf2 Seed(7) 2GB  [Low4/16]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(10) 1GB  [Low1/32]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(10) 1GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(10) 4GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(13) 1GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(13) 2GB  [Low1/16]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(15) 1GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(15) 2GB  [Low1/16]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(17) 1GB  BDayS128(26):all2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(17) 4GB  BDayS128(26):all2 p=2e-3 unusual 
RocRand te10 tf2 Seed(17) 4GB  [Low1/8]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(25) 2GB  [Low8/32]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(27) 1GB  [Low4/64]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(32) 4GB  [Low8/32]BDayS128(26):all2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(33) 2GB  [Low4/16]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(33) 2GB  [Low8/32]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(33) 4GB  [Low1/64]BDayS128(26):inc1 p=2e-3 unusual 
RocRand te10 tf2 Seed(33) 4GB  [Low1/64]BDayS128(26):inc2 p=9e-12 FAIL 
RocRand te10 tf2 Seed(34) 2GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(34) 4GB  [Low4/16]BDayS128(26):all2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(34) 4GB  [Low4/32]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(35) 2GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(35) 4GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(36) 4GB  [Low1/32]BDayS128(26):inc1 p=2e-3 unusual 
RocRand te10 tf2 Seed(36) 4GB  [Low1/32]BDayS128(26):inc2 p=9e-12 FAIL 
RocRand te10 tf2 Seed(38) 1GB  [Low1/8]BDayS128(26):inc1 p=2e-3 unusual 
RocRand te10 tf2 Seed(38) 1GB  [Low1/8]BDayS128(26):inc2 p=9e-12 FAIL 
RocRand te10 tf2 Seed(38) 1GB  [Low4/32]BDayS128(26):inc2 p=8e-6 suspicious 
RocRand te10 tf2 Seed(39) 8GB  [Low4/64]BDayS128(26):inc2 p=8e-6 suspicious

te10 を指定すると te0 や te1 に比べてかなりの頻度で非常に悪い結果が多く出現します。

rdrand に欠陥がある可能性があると考えて、PractRand でお勧めの sfc64 でも te10 でテストしてみました。
(gawk.exe を自分で用意してください。 example の「tee sfc64-te10-tf2.txt」は省略可能です。)

: repeat2.bat
@echo off
if not '%3' == '' goto SkipUsage
:Usage
 echo/
 echo example:
 echo "    repeat2 sfc64 10 2 | tee sfc64-te10-tf2.txt"
 echo/
goto EndOfBat
:SkipUsage
for %%P in ( 0 1 10 ) do if '%%P' == '%2' goto TeOK
goto Usage
:TeOK
for %%P in ( 0 1 2 ) do if '%%P' == '%3' goto TfOK
goto Usage
:TfOK
set AWK="/seed =/{s=$6}/bytes/{t=$2}/R=/{gsub(/ +/,\" \");sub(/ *R=.*p[ ~]= */,\" p=\");print \"%1\",\"Seed \" s,\"te%2 tf%3\",t \"GB\",$0}"
set S2= 1  2  3  4  5  6  7  8  9  a  b  c  d  e  f 10 11 12 13 14
set S4=15 16 17 18 19 1a 1b 1c 1d 1e 1f 20 21 22 23 24 25 26 27 28
for %%P in ( %S2% %S4% ) do Rng_test %1 -seed %%P -te %2 -tf %3 -tlmin 1GB -tlmax 4GB -multithreaded | gawk %AWK%
:EndOfBat
if not '%0' == 'repeat2' pause

sfc64 Seed 0x2 te10 tf1 2GB  BDayS128(26):all2 p=3e-3 unusual 
sfc64 Seed 0x9 te10 tf1 2GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x11 te10 tf1 1GB  [Low16/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x16 te10 tf1 4GB  [Low16/64]BDayS128(26):all2 p=8e-6 suspicious 
sfc64 Seed 0x18 te10 tf1 2GB  [Low4/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x19 te10 tf1 2GB  BDayS128(26):all2 p=3e-3 unusual 
sfc64 Seed 0x19 te10 tf1 2GB  [Low16/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x1a te10 tf1 2GB  [Low4/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x1f te10 tf1 1GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x1f te10 tf1 2GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x26 te10 tf1 4GB  [Low4/64]BDayS128(26):inc2 p=8e-6 suspicious

sfc64 Seed 0x5 te10 tf2 4GB  [Low1/8]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x7 te10 tf2 2GB  [Low1/16]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x7 te10 tf2 4GB  [Low1/16]BDayS128(26):inc1 p=2e-3 unusual 
sfc64 Seed 0x7 te10 tf2 4GB  [Low1/16]BDayS128(26):inc2 p=9e-12 FAIL 
sfc64 Seed 0x9 te10 tf2 2GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0xa te10 tf2 2GB  [Low8/32]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x11 te10 tf2 2GB  [Low1/32]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x12 te10 tf2 4GB  [Low4/32]BDayS128(26):inc1 p=3e-4 unusual 
sfc64 Seed 0x12 te10 tf2 4GB  [Low4/32]BDayS128(26):inc2 p=6e-18 FAIL ! 
sfc64 Seed 0x15 te10 tf2 4GB  [Low1/8]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x16 te10 tf2 4GB  [Low8/32]BDayS128(26):all2 p=8e-6 suspicious 
sfc64 Seed 0x17 te10 tf2 1GB  [Low8/32]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x18 te10 tf2 1GB  [Low1/8]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x18 te10 tf2 2GB  [Low4/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x1a te10 tf2 2GB  [Low4/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x1b te10 tf2 4GB  [Low8/64]BDayS128(26):inc1 p=2e-3 unusual 
sfc64 Seed 0x1b te10 tf2 4GB  [Low8/64]BDayS128(26):inc2 p=9e-12 FAIL 
sfc64 Seed 0x1c te10 tf2 1GB  [Low4/32]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x1e te10 tf2 2GB  [Low1/16]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x1f te10 tf2 1GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x1f te10 tf2 2GB  [Low1/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x20 te10 tf2 4GB  [Low1/32]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x21 te10 tf2 1GB  [Low4/16]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x21 te10 tf2 2GB  [Low8/64]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x22 te10 tf2 2GB  [Low1/16]BDayS128(26):inc2 p=8e-6 suspicious 
sfc64 Seed 0x26 te10 tf2 1GB  [Low8/64]BDayS128(26):inc1 p=2e-3 unusual 
sfc64 Seed 0x26 te10 tf2 1GB  [Low8/64]BDayS128(26):inc2 p=9e-12 FAIL

te10 でテストした場合、suspicious や FAIL と判定されてもその乱数が悪いとは限らないことが分かりました。

rnb.h

多くのコンパイラでコンパイル可能です。

define.h, bincode.h はこちら

/* rnb.h : coded by Isaku WADA(isaku@pb4.so-net.ne.jp) 2025
A hardware RNG.
gcc, clang, bcc5.5.1, vc5.0 .. vc2022
*/
#if !defined(RNB_H__)
 #define RNB_H__
 #include "define.h"
 #include "bincode.h"

typedef char Rnb_t;

 #if defined(_M_IX86) || defined(__i386)
static BINCODE__(RnbRand__)={
    0x0f,0xc7,0xf0, /* rdrand eax */
    0x72,4,         /* jc +4 */
    0xf3,0x90,      /* pause */
    0xeb,255&-9,    /* jmp -9 */
    0x0f,0xc7,0xf2, /* rdrand edx */
    0x72,4,         /* jc +4 */
    0xf3,0x90,      /* pause */
    0xeb,255&-9,    /* jmp -9 */
    0xc3};          /* ret */
 #else
static BINCODE__(RnbRand__)={
    0x48,0x0f,0xc7,0xf0, /* rdrand rax */
    0x72,4,              /* jc +4 */
    0xf3,0x90,           /* pause */
    0xeb,255&-10,        /* jmp -10 */
    0xc3};               /* ret */
 #endif
 #define RnbRand ((qword(*)(Rnb_t*))(size_t)RnbRand__)

static __inline int RnbHaveRdrand(void) {
    int v[4]; CpuidEcx0(v,7);
    return (v[1]>>18)&1;
}

static __inline const char*RnbInit(Rnb_t*p,qword s) {
    if (RnbHaveRdrand()==0) return"RDRAND";
    (void)p; (void)s; return 0;
}

#endif

#if !defined(ExRand) && !defined(ROB_H__)
 #define ExRand   RnbRand
 #define ExName  "RnbRand"
 #define ExInit   RnbInit
 #define Ex_t     Rnb_t
 #define ExAlign  1
#endif

rob.h

OpenMP が必要なため、Visual C++ では 2005 からコンパイル可能です。
直前の rnb.h と define.h, bincode.h が必要です。

define.h, bincode.h はこちら

/* rob.h : coded by Isaku WADA(isaku@pb4.so-net.ne.jp) 2025
A hardware RNG with OpenMP.
gcc, clang, vc2005 .. vc2022
msys2:clang -> clang -fopenmp foo.c -lwinmm -static -lomp
msys2:gcc   -> gcc   -fopenmp foo.c -lwinmm -static -lgomp
*/
#if !defined(ROB_H__)
 #define ROB_H__
 #if !defined(_OPENMP)
  #error OpenMP
 #endif
 #include <omp.h>
 #include "rnb.h"
 #define ROB_NUM 8192

typedef Rnb_t Rob_t;

static __inline qword RobRand(Rob_t*p) {
    static qword buf[ROB_NUM]; static int idx=ROB_NUM; int i; (void)p;
    if (omp_in_parallel()) return RnbRand(p);
    if (idx>=ROB_NUM) {
 #if defined(RO_NUMTH)
        #pragma omp parallel for num_threads(RO_NUMTH)
 #else
        #pragma omp parallel for
 #endif
        for (i=0;i<ROB_NUM;i++) buf[i]=RnbRand(p);
        idx=0;
    }
    return buf[idx++];
}

 #define RobInit RnbInit

#endif

#if !defined(ExRand)
 #define ExRand   RobRand
 #define ExName  "RobRand"
 #define ExInit   RobInit
 #define Ex_t     Rob_t
 #define ExAlign  1
#endif

［ExRand目次］［次：SplitMix64について］