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