Branch data Line data Source code
1 : : /* Unit tests for grand
2 : : * Copyright (C) 1995-1997 Peter Mattis, Spencer Kimball and Josh MacDonald
3 : : *
4 : : * SPDX-License-Identifier: LicenseRef-old-glib-tests
5 : : *
6 : : * This work is provided "as is"; redistribution and modification
7 : : * in whole or in part, in any medium, physical or electronic is
8 : : * permitted without restriction.
9 : : *
10 : : * This work is distributed in the hope that it will be useful,
11 : : * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 : : * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
13 : : *
14 : : * In no event shall the authors or contributors be liable for any
15 : : * direct, indirect, incidental, special, exemplary, or consequential
16 : : * damages (including, but not limited to, procurement of substitute
17 : : * goods or services; loss of use, data, or profits; or business
18 : : * interruption) however caused and on any theory of liability, whether
19 : : * in contract, strict liability, or tort (including negligence or
20 : : * otherwise) arising in any way out of the use of this software, even
21 : : * if advised of the possibility of such damage.
22 : : */
23 : :
24 : : #include "glib.h"
25 : :
26 : : /* Outputs tested against the reference implementation mt19937ar.c from
27 : : * http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/MT2002/emt19937ar.html
28 : : */
29 : :
30 : : /* Tests for a simple seed, first number is the seed */
31 : : const guint32 first_numbers[] =
32 : : {
33 : : 0x7a7a7a7a,
34 : : 0xfdcc2d54,
35 : : 0x3a279ceb,
36 : : 0xc4d39c33,
37 : : 0xf31895cd,
38 : : 0x46ca0afc,
39 : : 0x3f5484ff,
40 : : 0x54bc9557,
41 : : 0xed2c24b1,
42 : : 0x84062503,
43 : : 0x8f6404b3,
44 : : 0x599a94b3,
45 : : 0xe46d03d5,
46 : : 0x310beb78,
47 : : 0x7bee5d08,
48 : : 0x760d09be,
49 : : 0x59b6e163,
50 : : 0xbf6d16ec,
51 : : 0xcca5fb54,
52 : : 0x5de7259b,
53 : : 0x1696330c,
54 : : };
55 : :
56 : : /* array seed */
57 : : const guint32 seed_array[] =
58 : : {
59 : : 0x6553375f,
60 : : 0xd6b8d43b,
61 : : 0xa1e7667f,
62 : : 0x2b10117c
63 : : };
64 : :
65 : : /* tests for the array seed */
66 : : const guint32 array_outputs[] =
67 : : {
68 : : 0xc22b7dc3,
69 : : 0xfdecb8ae,
70 : : 0xb4af0738,
71 : : 0x516bc6e1,
72 : : 0x7e372e91,
73 : : 0x2d38ff80,
74 : : 0x6096494a,
75 : : 0xd162d5a8,
76 : : 0x3c0aaa0d,
77 : : 0x10e736ae
78 : : };
79 : :
80 : : static void
81 : 1 : test_rand (void)
82 : : {
83 : : guint n;
84 : : guint ones;
85 : : double proportion;
86 : : GRand *rand;
87 : : GRand *copy;
88 : :
89 : 1 : rand = g_rand_new_with_seed (first_numbers[0]);
90 : :
91 [ + + ]: 21 : for (n = 1; n < G_N_ELEMENTS (first_numbers); n++)
92 : 20 : g_assert_cmpuint (first_numbers[n], ==, g_rand_int (rand));
93 : :
94 : 1 : g_rand_set_seed (rand, 2);
95 : 1 : g_rand_set_seed_array (rand, seed_array, G_N_ELEMENTS (seed_array));
96 : :
97 [ + + ]: 11 : for (n = 0; n < G_N_ELEMENTS (array_outputs); n++)
98 : 10 : g_assert_cmpuint (array_outputs[n], ==, g_rand_int (rand));
99 : :
100 : 1 : copy = g_rand_copy (rand);
101 [ + + ]: 101 : for (n = 0; n < 100; n++)
102 : 100 : g_assert_cmpuint (g_rand_int (copy), ==, g_rand_int (rand));
103 : :
104 [ + + ]: 100000 : for (n = 1; n < 100000; n++)
105 : : {
106 : : gint32 i;
107 : : gdouble d;
108 : : gboolean b;
109 : :
110 : 99999 : i = g_rand_int_range (rand, 8,16);
111 : 99999 : g_assert_cmpint (i, >=, 8);
112 : 99999 : g_assert_cmpint (i, <, 16);
113 : :
114 : 99999 : i = g_random_int_range (8,16);
115 : 99999 : g_assert_cmpint (i, >=, 8);
116 : 99999 : g_assert_cmpint (i, <, 16);
117 : :
118 : 99999 : d = g_rand_double (rand);
119 : 99999 : g_assert_cmpfloat (d, >=, 0.0);
120 : 99999 : g_assert_cmpfloat (d, <, 1.0);
121 : :
122 : 99999 : d = g_random_double ();
123 : 99999 : g_assert_cmpfloat (d, >=, 0.0);
124 : 99999 : g_assert_cmpfloat (d, <, 1.0);
125 : :
126 : 99999 : d = g_rand_double_range (rand, -8, 32);
127 : 99999 : g_assert_cmpfloat (d, >=, -8.0);
128 : 99999 : g_assert_cmpfloat (d, <, 32.0);
129 : :
130 : 99999 : d = g_random_double_range (-8, 32);
131 : 99999 : g_assert_cmpfloat (d, >=, -8.0);
132 : 99999 : g_assert_cmpfloat (d, <, 32.0);
133 : :
134 : 99999 : b = g_random_boolean ();
135 : 99999 : g_assert_true (b == TRUE || b == FALSE);
136 : :
137 : 99999 : b = g_rand_boolean (rand);
138 : 99999 : g_assert_true (b == TRUE || b == FALSE);
139 : : }
140 : :
141 : : /* Statistical sanity check, count the number of ones
142 : : * when getting random numbers in range [0,3) and see
143 : : * that it must be semi-close to 0.25 with a VERY large
144 : : * probability */
145 : 1 : ones = 0;
146 [ + + ]: 100000 : for (n = 1; n < 100000; n++)
147 : : {
148 [ + + ]: 99999 : if (g_random_int_range (0, 4) == 1)
149 : 24847 : ones ++;
150 : : }
151 : :
152 : 1 : proportion = (double)ones / (double)100000;
153 : : /* 0.025 is overkill, but should suffice to test for some unreasonability */
154 : 1 : g_assert_cmpfloat (ABS (proportion - 0.25), <, 0.025);
155 : :
156 : 1 : g_rand_free (rand);
157 : 1 : g_rand_free (copy);
158 : 1 : }
159 : :
160 : : static void
161 : 1 : test_double_range (void)
162 : : {
163 : : gdouble d;
164 : :
165 : 1 : g_test_bug ("https://bugzilla.gnome.org/show_bug.cgi?id=502560");
166 : :
167 : 1 : d = g_random_double_range (-G_MAXDOUBLE, G_MAXDOUBLE);
168 : :
169 : 1 : g_assert_cmpfloat (-G_MAXDOUBLE, <=, d);
170 : 1 : g_assert_cmpfloat (d, <, G_MAXDOUBLE);
171 : 1 : }
172 : :
173 : : int
174 : 1 : main (int argc,
175 : : char *argv[])
176 : : {
177 : 1 : g_test_init (&argc, &argv, NULL);
178 : :
179 : 1 : g_test_add_func ("/rand/test-rand", test_rand);
180 : 1 : g_test_add_func ("/rand/double-range", test_double_range);
181 : :
182 : 1 : return g_test_run();
183 : : }
|