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
|
/*
* Generate testing csv files
*
* cl mt19937-test-data-gen.c randomkit.c
* -IC:\Anaconda\Lib\site-packages\numpy\core\include -IC:\Anaconda\include
* Advapi32.lib Kernel32.lib C:\Anaconda\libs\python36.lib -DRK_NO_WINCRYPT=1
*
*/
#include "randomkit.h"
#include <inttypes.h>
#include <stdio.h>
#define N 1000
int main() {
uint64_t sum = 0;
uint32_t seed = 0xDEADBEAF;
int i;
rk_state state;
rk_seed(seed, &state);
uint64_t store[N];
for (i = 0; i < N; i++) {
store[i] = (uint64_t)rk_random(&state);
}
FILE *fp;
fp = fopen("mt19937-testset-1.csv", "w");
if (fp == NULL) {
printf("Couldn't open file\n");
return -1;
}
fprintf(fp, "seed, 0x%" PRIx32 "\n", seed);
for (i = 0; i < N; i++) {
fprintf(fp, "%d, 0x%" PRIx64 "\n", i, store[i]);
if (i == 999) {
printf("%d, 0x%" PRIx64 "\n", i, store[i]);
}
}
fclose(fp);
seed = 0;
rk_seed(seed, &state);
for (i = 0; i < N; i++) {
store[i] = (uint64_t)rk_random(&state);
}
fp = fopen("mt19937-testset-2.csv", "w");
if (fp == NULL) {
printf("Couldn't open file\n");
return -1;
}
fprintf(fp, "seed, 0x%" PRIx32 "\n", seed);
for (i = 0; i < N; i++) {
fprintf(fp, "%d, 0x%" PRIx64 "\n", i, store[i]);
if (i == 999) {
printf("%d, 0x%" PRIx64 "\n", i, store[i]);
}
}
fclose(fp);
}
|
