]>
Commit | Line | Data |
---|---|---|
9bf714c2 UD |
1 | // Optimizations for random number extensions, x86 version -*- C++ -*- |
2 | ||
8d9254fc | 3 | // Copyright (C) 2012-2020 Free Software Foundation, Inc. |
9bf714c2 UD |
4 | // |
5 | // This file is part of the GNU ISO C++ Library. This library is free | |
6 | // software; you can redistribute it and/or modify it under the | |
7 | // terms of the GNU General Public License as published by the | |
8 | // Free Software Foundation; either version 3, or (at your option) | |
9 | // any later version. | |
10 | ||
11 | // This library is distributed in the hope that it will be useful, | |
12 | // but WITHOUT ANY WARRANTY; without even the implied warranty of | |
13 | // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
14 | // GNU General Public License for more details. | |
15 | ||
16 | // Under Section 7 of GPL version 3, you are granted additional | |
17 | // permissions described in the GCC Runtime Library Exception, version | |
18 | // 3.1, as published by the Free Software Foundation. | |
19 | ||
20 | // You should have received a copy of the GNU General Public License and | |
21 | // a copy of the GCC Runtime Library Exception along with this program; | |
22 | // see the files COPYING3 and COPYING.RUNTIME respectively. If not, see | |
23 | // <http://www.gnu.org/licenses/>. | |
24 | ||
25 | /** @file ext/random.tcc | |
26 | * This is an internal header file, included by other library headers. | |
27 | * Do not attempt to use it directly. @headername{ext/random} | |
28 | */ | |
29 | ||
30 | #ifndef _EXT_OPT_RANDOM_H | |
31 | #define _EXT_OPT_RANDOM_H 1 | |
32 | ||
33 | #pragma GCC system_header | |
34 | ||
eeeef8f4 PC |
35 | #if __BYTE_ORDER__ == __ORDER_LITTLE_ENDIAN__ |
36 | ||
ecfce0b8 | 37 | #ifdef __SSE2__ |
9bf714c2 UD |
38 | |
39 | namespace __gnu_cxx _GLIBCXX_VISIBILITY(default) | |
40 | { | |
41 | _GLIBCXX_BEGIN_NAMESPACE_VERSION | |
42 | ||
43 | namespace { | |
44 | ||
45 | template<size_t __sl1, size_t __sl2, size_t __sr1, size_t __sr2, | |
46 | uint32_t __msk1, uint32_t __msk2, uint32_t __msk3, uint32_t __msk4> | |
47 | inline __m128i __sse2_recursion(__m128i __a, __m128i __b, | |
48 | __m128i __c, __m128i __d) | |
49 | { | |
50 | __m128i __y = _mm_srli_epi32(__b, __sr1); | |
51 | __m128i __z = _mm_srli_si128(__c, __sr2); | |
52 | __m128i __v = _mm_slli_epi32(__d, __sl1); | |
53 | __z = _mm_xor_si128(__z, __a); | |
54 | __z = _mm_xor_si128(__z, __v); | |
55 | __m128i __x = _mm_slli_si128(__a, __sl2); | |
56 | __y = _mm_and_si128(__y, _mm_set_epi32(__msk4, __msk3, __msk2, __msk1)); | |
57 | __z = _mm_xor_si128(__z, __x); | |
58 | return _mm_xor_si128(__z, __y); | |
59 | } | |
60 | ||
61 | } | |
62 | ||
63 | ||
64 | #define _GLIBCXX_OPT_HAVE_RANDOM_SFMT_GEN_READ 1 | |
65 | template<typename _UIntType, size_t __m, | |
66 | size_t __pos1, size_t __sl1, size_t __sl2, | |
67 | size_t __sr1, size_t __sr2, | |
68 | uint32_t __msk1, uint32_t __msk2, | |
69 | uint32_t __msk3, uint32_t __msk4, | |
70 | uint32_t __parity1, uint32_t __parity2, | |
71 | uint32_t __parity3, uint32_t __parity4> | |
72 | void simd_fast_mersenne_twister_engine<_UIntType, __m, | |
73 | __pos1, __sl1, __sl2, __sr1, __sr2, | |
74 | __msk1, __msk2, __msk3, __msk4, | |
75 | __parity1, __parity2, __parity3, | |
76 | __parity4>:: | |
77 | _M_gen_rand(void) | |
78 | { | |
79 | __m128i __r1 = _mm_load_si128(&_M_state[_M_nstate - 2]); | |
80 | __m128i __r2 = _mm_load_si128(&_M_state[_M_nstate - 1]); | |
81 | ||
82 | size_t __i; | |
83 | for (__i = 0; __i < _M_nstate - __pos1; ++__i) | |
84 | { | |
85 | __m128i __r = __sse2_recursion<__sl1, __sl2, __sr1, __sr2, | |
86 | __msk1, __msk2, __msk3, __msk4> | |
87 | (_M_state[__i], _M_state[__i + __pos1], __r1, __r2); | |
88 | _mm_store_si128(&_M_state[__i], __r); | |
89 | __r1 = __r2; | |
90 | __r2 = __r; | |
91 | } | |
92 | for (; __i < _M_nstate; ++__i) | |
93 | { | |
94 | __m128i __r = __sse2_recursion<__sl1, __sl2, __sr1, __sr2, | |
95 | __msk1, __msk2, __msk3, __msk4> | |
96 | (_M_state[__i], _M_state[__i + __pos1 - _M_nstate], __r1, __r2); | |
97 | _mm_store_si128(&_M_state[__i], __r); | |
98 | __r1 = __r2; | |
99 | __r2 = __r; | |
100 | } | |
101 | ||
102 | _M_pos = 0; | |
103 | } | |
104 | ||
105 | ||
106 | #define _GLIBCXX_OPT_HAVE_RANDOM_SFMT_OPERATOREQUAL 1 | |
107 | template<typename _UIntType, size_t __m, | |
108 | size_t __pos1, size_t __sl1, size_t __sl2, | |
109 | size_t __sr1, size_t __sr2, | |
110 | uint32_t __msk1, uint32_t __msk2, | |
111 | uint32_t __msk3, uint32_t __msk4, | |
112 | uint32_t __parity1, uint32_t __parity2, | |
113 | uint32_t __parity3, uint32_t __parity4> | |
114 | bool | |
115 | operator==(const __gnu_cxx::simd_fast_mersenne_twister_engine<_UIntType, | |
116 | __m, __pos1, __sl1, __sl2, __sr1, __sr2, | |
117 | __msk1, __msk2, __msk3, __msk4, | |
118 | __parity1, __parity2, __parity3, __parity4>& __lhs, | |
119 | const __gnu_cxx::simd_fast_mersenne_twister_engine<_UIntType, | |
120 | __m, __pos1, __sl1, __sl2, __sr1, __sr2, | |
121 | __msk1, __msk2, __msk3, __msk4, | |
122 | __parity1, __parity2, __parity3, __parity4>& __rhs) | |
123 | { | |
124 | __m128i __res = _mm_cmpeq_epi8(__lhs._M_state[0], __rhs._M_state[0]); | |
125 | for (size_t __i = 1; __i < __lhs._M_nstate; ++__i) | |
126 | __res = _mm_and_si128(__res, _mm_cmpeq_epi8(__lhs._M_state[__i], | |
127 | __rhs._M_state[__i])); | |
128 | return (_mm_movemask_epi8(__res) == 0xffff | |
129 | && __lhs._M_pos == __rhs._M_pos); | |
130 | } | |
131 | ||
132 | ||
133 | _GLIBCXX_END_NAMESPACE_VERSION | |
134 | } // namespace | |
135 | ||
ecfce0b8 | 136 | #endif // __SSE2__ |
9bf714c2 | 137 | |
eeeef8f4 PC |
138 | #endif // __BYTE_ORDER__ == __ORDER_LITTLE_ENDIAN__ |
139 | ||
9bf714c2 | 140 | #endif // _EXT_OPT_RANDOM_H |