1*6467f958SSadaf Ebrahimi //
2*6467f958SSadaf Ebrahimi // Copyright (c) 2017 The Khronos Group Inc.
3*6467f958SSadaf Ebrahimi //
4*6467f958SSadaf Ebrahimi // Licensed under the Apache License, Version 2.0 (the "License");
5*6467f958SSadaf Ebrahimi // you may not use this file except in compliance with the License.
6*6467f958SSadaf Ebrahimi // You may obtain a copy of the License at
7*6467f958SSadaf Ebrahimi //
8*6467f958SSadaf Ebrahimi // http://www.apache.org/licenses/LICENSE-2.0
9*6467f958SSadaf Ebrahimi //
10*6467f958SSadaf Ebrahimi // Unless required by applicable law or agreed to in writing, software
11*6467f958SSadaf Ebrahimi // distributed under the License is distributed on an "AS IS" BASIS,
12*6467f958SSadaf Ebrahimi // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13*6467f958SSadaf Ebrahimi // See the License for the specific language governing permissions and
14*6467f958SSadaf Ebrahimi // limitations under the License.
15*6467f958SSadaf Ebrahimi //
16*6467f958SSadaf Ebrahimi #include "mt19937.h"
17*6467f958SSadaf Ebrahimi #include <stdio.h>
18*6467f958SSadaf Ebrahimi
main(void)19*6467f958SSadaf Ebrahimi int main(void)
20*6467f958SSadaf Ebrahimi {
21*6467f958SSadaf Ebrahimi MTdata d = init_genrand(42);
22*6467f958SSadaf Ebrahimi int i;
23*6467f958SSadaf Ebrahimi const cl_uint reference[16] = {
24*6467f958SSadaf Ebrahimi 0x5fe1dc66, 0x8b255210, 0x0380b0c8, 0xc87d2ce4, 0x55c31f24, 0x8bcd21ab,
25*6467f958SSadaf Ebrahimi 0x14d5fef5, 0x9416d2b6, 0xdf875de9, 0x00517d76, 0xd861c944, 0xa7676404,
26*6467f958SSadaf Ebrahimi 0x5491aff4, 0x67616209, 0xc368b3fb, 0x929dfc92
27*6467f958SSadaf Ebrahimi };
28*6467f958SSadaf Ebrahimi int errcount = 0;
29*6467f958SSadaf Ebrahimi
30*6467f958SSadaf Ebrahimi for (i = 0; i < 65536; i++)
31*6467f958SSadaf Ebrahimi {
32*6467f958SSadaf Ebrahimi cl_uint u = genrand_int32(d);
33*6467f958SSadaf Ebrahimi if (0 == (i & 4095))
34*6467f958SSadaf Ebrahimi {
35*6467f958SSadaf Ebrahimi if (u != reference[i >> 12])
36*6467f958SSadaf Ebrahimi {
37*6467f958SSadaf Ebrahimi printf("ERROR: expected *0x%8.8x at %d. Got 0x%8.8x\n",
38*6467f958SSadaf Ebrahimi reference[i >> 12], i, u);
39*6467f958SSadaf Ebrahimi errcount++;
40*6467f958SSadaf Ebrahimi }
41*6467f958SSadaf Ebrahimi }
42*6467f958SSadaf Ebrahimi }
43*6467f958SSadaf Ebrahimi
44*6467f958SSadaf Ebrahimi free_mtdata(d);
45*6467f958SSadaf Ebrahimi
46*6467f958SSadaf Ebrahimi if (errcount)
47*6467f958SSadaf Ebrahimi printf("mt19937 test failed.\n");
48*6467f958SSadaf Ebrahimi else
49*6467f958SSadaf Ebrahimi printf("mt19937 test passed.\n");
50*6467f958SSadaf Ebrahimi
51*6467f958SSadaf Ebrahimi
52*6467f958SSadaf Ebrahimi return 0;
53*6467f958SSadaf Ebrahimi }