1/*
2 * This file derives from SFMT 1.3.3
3 * (http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/SFMT/index.html), which was
4 * released under the terms of the following license:
5 *
6 *   Copyright (c) 2006,2007 Mutsuo Saito, Makoto Matsumoto and Hiroshima
7 *   University. All rights reserved.
8 *
9 *   Redistribution and use in source and binary forms, with or without
10 *   modification, are permitted provided that the following conditions are
11 *   met:
12 *
13 *       * Redistributions of source code must retain the above copyright
14 *         notice, this list of conditions and the following disclaimer.
15 *       * Redistributions in binary form must reproduce the above
16 *         copyright notice, this list of conditions and the following
17 *         disclaimer in the documentation and/or other materials provided
18 *         with the distribution.
19 *       * Neither the name of the Hiroshima University nor the names of
20 *         its contributors may be used to endorse or promote products
21 *         derived from this software without specific prior written
22 *         permission.
23 *
24 *   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
25 *   "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
26 *   LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
27 *   A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
28 *   OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
29 *   SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
30 *   LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
31 *   DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
32 *   THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
33 *   (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
34 *   OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
35 */
36#ifndef SFMT_PARAMS_H
37#define SFMT_PARAMS_H
38
39#if !defined(MEXP)
40#ifdef __GNUC__
41  #warning "MEXP is not defined. I assume MEXP is 19937."
42#endif
43  #define MEXP 19937
44#endif
45/*-----------------
46  BASIC DEFINITIONS
47  -----------------*/
48/** Mersenne Exponent. The period of the sequence
49 *  is a multiple of 2^MEXP-1.
50 * #define MEXP 19937 */
51/** SFMT generator has an internal state array of 128-bit integers,
52 * and N is its size. */
53#define N (MEXP / 128 + 1)
54/** N32 is the size of internal state array when regarded as an array
55 * of 32-bit integers.*/
56#define N32 (N * 4)
57/** N64 is the size of internal state array when regarded as an array
58 * of 64-bit integers.*/
59#define N64 (N * 2)
60
61/*----------------------
62  the parameters of SFMT
63  following definitions are in paramsXXXX.h file.
64  ----------------------*/
65/** the pick up position of the array.
66#define POS1 122
67*/
68
69/** the parameter of shift left as four 32-bit registers.
70#define SL1 18
71 */
72
73/** the parameter of shift left as one 128-bit register.
74 * The 128-bit integer is shifted by (SL2 * 8) bits.
75#define SL2 1
76*/
77
78/** the parameter of shift right as four 32-bit registers.
79#define SR1 11
80*/
81
82/** the parameter of shift right as one 128-bit register.
83 * The 128-bit integer is shifted by (SL2 * 8) bits.
84#define SR2 1
85*/
86
87/** A bitmask, used in the recursion.  These parameters are introduced
88 * to break symmetry of SIMD.
89#define MSK1 0xdfffffefU
90#define MSK2 0xddfecb7fU
91#define MSK3 0xbffaffffU
92#define MSK4 0xbffffff6U
93*/
94
95/** These definitions are part of a 128-bit period certification vector.
96#define PARITY1	0x00000001U
97#define PARITY2	0x00000000U
98#define PARITY3	0x00000000U
99#define PARITY4	0xc98e126aU
100*/
101
102#if MEXP == 607
103  #include "test/SFMT-params607.h"
104#elif MEXP == 1279
105  #include "test/SFMT-params1279.h"
106#elif MEXP == 2281
107  #include "test/SFMT-params2281.h"
108#elif MEXP == 4253
109  #include "test/SFMT-params4253.h"
110#elif MEXP == 11213
111  #include "test/SFMT-params11213.h"
112#elif MEXP == 19937
113  #include "test/SFMT-params19937.h"
114#elif MEXP == 44497
115  #include "test/SFMT-params44497.h"
116#elif MEXP == 86243
117  #include "test/SFMT-params86243.h"
118#elif MEXP == 132049
119  #include "test/SFMT-params132049.h"
120#elif MEXP == 216091
121  #include "test/SFMT-params216091.h"
122#else
123#ifdef __GNUC__
124  #error "MEXP is not valid."
125  #undef MEXP
126#else
127  #undef MEXP
128#endif
129
130#endif
131
132#endif /* SFMT_PARAMS_H */
133