Chris@16
|
1 ///////////////////////////////////////////////////////////////
|
Chris@16
|
2 // Copyright 2013 John Maddock. Distributed under the Boost
|
Chris@16
|
3 // Software License, Version 1.0. (See accompanying file
|
Chris@16
|
4 // LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_
|
Chris@16
|
5 //
|
Chris@16
|
6 // Comparison operators for cpp_int_backend:
|
Chris@16
|
7 //
|
Chris@16
|
8 #ifndef BOOST_MP_DETAIL_BITSCAN_HPP
|
Chris@16
|
9 #define BOOST_MP_DETAIL_BITSCAN_HPP
|
Chris@16
|
10
|
Chris@16
|
11 namespace boost{ namespace multiprecision{ namespace detail{
|
Chris@16
|
12
|
Chris@16
|
13 template <class Unsigned>
|
Chris@16
|
14 inline unsigned find_lsb(Unsigned mask, const mpl::int_<0>&)
|
Chris@16
|
15 {
|
Chris@16
|
16 unsigned result = 0;
|
Chris@16
|
17 while(!(mask & 1u))
|
Chris@16
|
18 {
|
Chris@16
|
19 mask >>= 1;
|
Chris@16
|
20 ++result;
|
Chris@16
|
21 }
|
Chris@16
|
22 return result;
|
Chris@16
|
23 }
|
Chris@16
|
24
|
Chris@16
|
25 template <class Unsigned>
|
Chris@16
|
26 inline unsigned find_msb(Unsigned mask, const mpl::int_<0>&)
|
Chris@16
|
27 {
|
Chris@16
|
28 unsigned index = 0;
|
Chris@16
|
29 while(mask)
|
Chris@16
|
30 {
|
Chris@16
|
31 ++index;
|
Chris@16
|
32 mask >>= 1;
|
Chris@16
|
33 }
|
Chris@16
|
34 return --index;
|
Chris@16
|
35 }
|
Chris@16
|
36
|
Chris@16
|
37 #if defined(BOOST_MSVC) && (defined(_M_IX86) || defined(_M_X64))
|
Chris@16
|
38 BOOST_FORCEINLINE unsigned find_lsb(unsigned long mask, const mpl::int_<1>&)
|
Chris@16
|
39 {
|
Chris@16
|
40 unsigned long result;
|
Chris@16
|
41 _BitScanForward(&result, mask);
|
Chris@16
|
42 return result;
|
Chris@16
|
43 }
|
Chris@16
|
44
|
Chris@16
|
45 BOOST_FORCEINLINE unsigned find_msb(unsigned long mask, const mpl::int_<1>&)
|
Chris@16
|
46 {
|
Chris@16
|
47 unsigned long result;
|
Chris@16
|
48 _BitScanReverse(&result, mask);
|
Chris@16
|
49 return result;
|
Chris@16
|
50 }
|
Chris@16
|
51 #ifdef _M_X64
|
Chris@16
|
52 BOOST_FORCEINLINE unsigned find_lsb(unsigned __int64 mask, const mpl::int_<2>&)
|
Chris@16
|
53 {
|
Chris@16
|
54 unsigned long result;
|
Chris@16
|
55 _BitScanForward64(&result, mask);
|
Chris@16
|
56 return result;
|
Chris@16
|
57 }
|
Chris@16
|
58 template <class Unsigned>
|
Chris@16
|
59 BOOST_FORCEINLINE unsigned find_msb(Unsigned mask, const mpl::int_<2>&)
|
Chris@16
|
60 {
|
Chris@16
|
61 unsigned long result;
|
Chris@16
|
62 _BitScanReverse64(&result, mask);
|
Chris@16
|
63 return result;
|
Chris@16
|
64 }
|
Chris@16
|
65 #endif
|
Chris@16
|
66
|
Chris@16
|
67 template <class Unsigned>
|
Chris@16
|
68 BOOST_FORCEINLINE unsigned find_lsb(Unsigned mask)
|
Chris@16
|
69 {
|
Chris@16
|
70 typedef typename make_unsigned<Unsigned>::type ui_type;
|
Chris@16
|
71 typedef typename mpl::if_c<
|
Chris@16
|
72 sizeof(Unsigned) <= sizeof(unsigned long),
|
Chris@16
|
73 mpl::int_<1>,
|
Chris@16
|
74 #ifdef _M_X64
|
Chris@16
|
75 typename mpl::if_c<
|
Chris@16
|
76 sizeof(Unsigned) <= sizeof(__int64),
|
Chris@16
|
77 mpl::int_<2>,
|
Chris@16
|
78 mpl::int_<0>
|
Chris@16
|
79 >::type
|
Chris@16
|
80 #else
|
Chris@16
|
81 mpl::int_<0>
|
Chris@16
|
82 #endif
|
Chris@16
|
83 >::type tag_type;
|
Chris@16
|
84 return find_lsb(static_cast<ui_type>(mask), tag_type());
|
Chris@16
|
85 }
|
Chris@16
|
86
|
Chris@16
|
87 template <class Unsigned>
|
Chris@16
|
88 BOOST_FORCEINLINE unsigned find_msb(Unsigned mask)
|
Chris@16
|
89 {
|
Chris@16
|
90 typedef typename make_unsigned<Unsigned>::type ui_type;
|
Chris@16
|
91 typedef typename mpl::if_c<
|
Chris@16
|
92 sizeof(Unsigned) <= sizeof(unsigned long),
|
Chris@16
|
93 mpl::int_<1>,
|
Chris@16
|
94 #ifdef _M_X64
|
Chris@16
|
95 typename mpl::if_c<
|
Chris@16
|
96 sizeof(Unsigned) <= sizeof(__int64),
|
Chris@16
|
97 mpl::int_<2>,
|
Chris@16
|
98 mpl::int_<0>
|
Chris@16
|
99 >::type
|
Chris@16
|
100 #else
|
Chris@16
|
101 mpl::int_<0>
|
Chris@16
|
102 #endif
|
Chris@16
|
103 >::type tag_type;
|
Chris@16
|
104 return find_msb(static_cast<ui_type>(mask), tag_type());
|
Chris@16
|
105 }
|
Chris@16
|
106
|
Chris@16
|
107 #elif defined(BOOST_GCC) || defined(__clang__) || (defined(BOOST_INTEL) && defined(__GNUC__))
|
Chris@16
|
108
|
Chris@16
|
109 BOOST_FORCEINLINE unsigned find_lsb(unsigned mask, mpl::int_<1> const&)
|
Chris@16
|
110 {
|
Chris@16
|
111 return __builtin_ctz(mask);
|
Chris@16
|
112 }
|
Chris@16
|
113 BOOST_FORCEINLINE unsigned find_lsb(unsigned long mask, mpl::int_<2> const&)
|
Chris@16
|
114 {
|
Chris@16
|
115 return __builtin_ctzl(mask);
|
Chris@16
|
116 }
|
Chris@16
|
117 BOOST_FORCEINLINE unsigned find_lsb(unsigned long long mask, mpl::int_<3> const&)
|
Chris@16
|
118 {
|
Chris@16
|
119 return __builtin_ctzll(mask);
|
Chris@16
|
120 }
|
Chris@16
|
121 BOOST_FORCEINLINE unsigned find_msb(unsigned mask, mpl::int_<1> const&)
|
Chris@16
|
122 {
|
Chris@16
|
123 return sizeof(unsigned) * CHAR_BIT - 1 - __builtin_clz(mask);
|
Chris@16
|
124 }
|
Chris@16
|
125 BOOST_FORCEINLINE unsigned find_msb(unsigned long mask, mpl::int_<2> const&)
|
Chris@16
|
126 {
|
Chris@16
|
127 return sizeof(unsigned long) * CHAR_BIT - 1 - __builtin_clzl(mask);
|
Chris@16
|
128 }
|
Chris@16
|
129 BOOST_FORCEINLINE unsigned find_msb(unsigned long long mask, mpl::int_<3> const&)
|
Chris@16
|
130 {
|
Chris@16
|
131 return sizeof(unsigned long long) * CHAR_BIT - 1 - __builtin_clzll(mask);
|
Chris@16
|
132 }
|
Chris@16
|
133
|
Chris@16
|
134 template <class Unsigned>
|
Chris@16
|
135 BOOST_FORCEINLINE unsigned find_lsb(Unsigned mask)
|
Chris@16
|
136 {
|
Chris@16
|
137 typedef typename make_unsigned<Unsigned>::type ui_type;
|
Chris@16
|
138 typedef typename mpl::if_c<
|
Chris@16
|
139 sizeof(Unsigned) <= sizeof(unsigned),
|
Chris@16
|
140 mpl::int_<1>,
|
Chris@16
|
141 typename mpl::if_c<
|
Chris@16
|
142 sizeof(Unsigned) <= sizeof(unsigned long),
|
Chris@16
|
143 mpl::int_<2>,
|
Chris@16
|
144 typename mpl::if_c<
|
Chris@16
|
145 sizeof(Unsigned) <= sizeof(unsigned long long),
|
Chris@16
|
146 mpl::int_<3>,
|
Chris@16
|
147 mpl::int_<0>
|
Chris@16
|
148 >::type
|
Chris@16
|
149 >::type
|
Chris@16
|
150 >::type tag_type;
|
Chris@16
|
151 return find_lsb(static_cast<ui_type>(mask), tag_type());
|
Chris@16
|
152 }
|
Chris@16
|
153 template <class Unsigned>
|
Chris@16
|
154 BOOST_FORCEINLINE unsigned find_msb(Unsigned mask)
|
Chris@16
|
155 {
|
Chris@16
|
156 typedef typename make_unsigned<Unsigned>::type ui_type;
|
Chris@16
|
157 typedef typename mpl::if_c<
|
Chris@16
|
158 sizeof(Unsigned) <= sizeof(unsigned),
|
Chris@16
|
159 mpl::int_<1>,
|
Chris@16
|
160 typename mpl::if_c<
|
Chris@16
|
161 sizeof(Unsigned) <= sizeof(unsigned long),
|
Chris@16
|
162 mpl::int_<2>,
|
Chris@16
|
163 typename mpl::if_c<
|
Chris@16
|
164 sizeof(Unsigned) <= sizeof(unsigned long long),
|
Chris@16
|
165 mpl::int_<3>,
|
Chris@16
|
166 mpl::int_<0>
|
Chris@16
|
167 >::type
|
Chris@16
|
168 >::type
|
Chris@16
|
169 >::type tag_type;
|
Chris@16
|
170 return find_msb(static_cast<ui_type>(mask), tag_type());
|
Chris@16
|
171 }
|
Chris@16
|
172 #elif defined(BOOST_INTEL)
|
Chris@16
|
173 BOOST_FORCEINLINE unsigned find_lsb(unsigned mask, mpl::int_<1> const&)
|
Chris@16
|
174 {
|
Chris@16
|
175 return _bit_scan_forward(mask);
|
Chris@16
|
176 }
|
Chris@16
|
177 BOOST_FORCEINLINE unsigned find_msb(unsigned mask, mpl::int_<1> const&)
|
Chris@16
|
178 {
|
Chris@16
|
179 return _bit_scan_reverse(mask);
|
Chris@16
|
180 }
|
Chris@16
|
181 template <class Unsigned>
|
Chris@16
|
182 BOOST_FORCEINLINE unsigned find_lsb(Unsigned mask)
|
Chris@16
|
183 {
|
Chris@16
|
184 typedef typename make_unsigned<Unsigned>::type ui_type;
|
Chris@16
|
185 typedef typename mpl::if_c<
|
Chris@16
|
186 sizeof(Unsigned) <= sizeof(unsigned),
|
Chris@16
|
187 mpl::int_<1>,
|
Chris@16
|
188 mpl::int_<0>
|
Chris@16
|
189 >::type tag_type;
|
Chris@16
|
190 return find_lsb(static_cast<ui_type>(mask), tag_type());
|
Chris@16
|
191 }
|
Chris@16
|
192 template <class Unsigned>
|
Chris@16
|
193 BOOST_FORCEINLINE unsigned find_msb(Unsigned mask)
|
Chris@16
|
194 {
|
Chris@16
|
195 typedef typename make_unsigned<Unsigned>::type ui_type;
|
Chris@16
|
196 typedef typename mpl::if_c<
|
Chris@16
|
197 sizeof(Unsigned) <= sizeof(unsigned),
|
Chris@16
|
198 mpl::int_<1>,
|
Chris@16
|
199 mpl::int_<0>
|
Chris@16
|
200 >::type tag_type;
|
Chris@16
|
201 return find_msb(static_cast<ui_type>(mask), tag_type());
|
Chris@16
|
202 }
|
Chris@16
|
203 #else
|
Chris@16
|
204 template <class Unsigned>
|
Chris@16
|
205 BOOST_FORCEINLINE unsigned find_lsb(Unsigned mask)
|
Chris@16
|
206 {
|
Chris@16
|
207 return find_lsb(mask, mpl::int_<0>());
|
Chris@16
|
208 }
|
Chris@16
|
209 template <class Unsigned>
|
Chris@16
|
210 BOOST_FORCEINLINE unsigned find_msb(Unsigned mask)
|
Chris@16
|
211 {
|
Chris@16
|
212 return find_msb(mask, mpl::int_<0>());
|
Chris@16
|
213 }
|
Chris@16
|
214 #endif
|
Chris@16
|
215
|
Chris@16
|
216 }}}
|
Chris@16
|
217
|
Chris@16
|
218 #endif
|
Chris@16
|
219
|