|
1 /* boost random/xor_combine.hpp header file |
|
2 * |
|
3 * Copyright Jens Maurer 2002 |
|
4 * Distributed under the Boost Software License, Version 1.0. (See |
|
5 * accompanying file LICENSE_1_0.txt or copy at |
|
6 * http://www.boost.org/LICENSE_1_0.txt) |
|
7 * |
|
8 * See http://www.boost.org for most recent version including documentation. |
|
9 * |
|
10 * $Id: xor_combine.hpp,v 1.13 2005/05/21 15:57:00 dgregor Exp $ |
|
11 * |
|
12 */ |
|
13 |
|
14 #ifndef BOOST_RANDOM_XOR_COMBINE_HPP |
|
15 #define BOOST_RANDOM_XOR_COMBINE_HPP |
|
16 |
|
17 #include <iostream> |
|
18 #include <cassert> |
|
19 #include <algorithm> // for std::min and std::max |
|
20 #include <boost/config.hpp> |
|
21 #include <boost/limits.hpp> |
|
22 #include <boost/static_assert.hpp> |
|
23 #include <boost/cstdint.hpp> // uint32_t |
|
24 |
|
25 |
|
26 namespace boost { |
|
27 namespace random { |
|
28 |
|
29 template<class URNG1, int s1, class URNG2, int s2, |
|
30 #ifndef BOOST_NO_DEPENDENT_TYPES_IN_TEMPLATE_VALUE_PARAMETERS |
|
31 typename URNG1::result_type |
|
32 #else |
|
33 uint32_t |
|
34 #endif |
|
35 val = 0> |
|
36 class xor_combine |
|
37 { |
|
38 public: |
|
39 typedef URNG1 base1_type; |
|
40 typedef URNG2 base2_type; |
|
41 typedef typename base1_type::result_type result_type; |
|
42 |
|
43 BOOST_STATIC_CONSTANT(bool, has_fixed_range = false); |
|
44 BOOST_STATIC_CONSTANT(int, shift1 = s1); |
|
45 BOOST_STATIC_CONSTANT(int, shfit2 = s2); |
|
46 |
|
47 xor_combine() : _rng1(), _rng2() |
|
48 { } |
|
49 xor_combine(const base1_type & rng1, const base2_type & rng2) |
|
50 : _rng1(rng1), _rng2(rng2) { } |
|
51 template<class It> xor_combine(It& first, It last) |
|
52 : _rng1(first, last), _rng2( /* advanced by other call */ first, last) { } |
|
53 void seed() { _rng1.seed(); _rng2.seed(); } |
|
54 template<class It> void seed(It& first, It last) |
|
55 { |
|
56 _rng1.seed(first, last); |
|
57 _rng2.seed(first, last); |
|
58 } |
|
59 |
|
60 const base1_type& base1() { return _rng1; } |
|
61 const base2_type& base2() { return _rng2; } |
|
62 |
|
63 result_type operator()() |
|
64 { |
|
65 // MSVC fails BOOST_STATIC_ASSERT with std::numeric_limits at class scope |
|
66 #if !defined(BOOST_NO_LIMITS_COMPILE_TIME_CONSTANTS) && !(defined(BOOST_MSVC) && BOOST_MSVC <= 1300) |
|
67 BOOST_STATIC_ASSERT(std::numeric_limits<typename base1_type::result_type>::is_integer); |
|
68 BOOST_STATIC_ASSERT(std::numeric_limits<typename base2_type::result_type>::is_integer); |
|
69 BOOST_STATIC_ASSERT(std::numeric_limits<typename base1_type::result_type>::digits >= std::numeric_limits<typename base2_type::result_type>::digits); |
|
70 #endif |
|
71 return (_rng1() << s1) ^ (_rng2() << s2); |
|
72 } |
|
73 |
|
74 result_type min BOOST_PREVENT_MACRO_SUBSTITUTION () const { return std::min BOOST_PREVENT_MACRO_SUBSTITUTION((_rng1.min)(), (_rng2.min)()); } |
|
75 result_type max BOOST_PREVENT_MACRO_SUBSTITUTION () const { return std::max BOOST_PREVENT_MACRO_SUBSTITUTION((_rng1.min)(), (_rng2.max)()); } |
|
76 static bool validation(result_type x) { return val == x; } |
|
77 |
|
78 #ifndef BOOST_NO_OPERATORS_IN_NAMESPACE |
|
79 |
|
80 #ifndef BOOST_NO_MEMBER_TEMPLATE_FRIENDS |
|
81 template<class CharT, class Traits> |
|
82 friend std::basic_ostream<CharT,Traits>& |
|
83 operator<<(std::basic_ostream<CharT,Traits>& os, const xor_combine& s) |
|
84 { |
|
85 os << s._rng1 << " " << s._rng2 << " "; |
|
86 return os; |
|
87 } |
|
88 |
|
89 template<class CharT, class Traits> |
|
90 friend std::basic_istream<CharT,Traits>& |
|
91 operator>>(std::basic_istream<CharT,Traits>& is, xor_combine& s) |
|
92 { |
|
93 is >> s._rng1 >> std::ws >> s._rng2 >> std::ws; |
|
94 return is; |
|
95 } |
|
96 #endif |
|
97 |
|
98 friend bool operator==(const xor_combine& x, const xor_combine& y) |
|
99 { return x._rng1 == y._rng1 && x._rng2 == y._rng2; } |
|
100 friend bool operator!=(const xor_combine& x, const xor_combine& y) |
|
101 { return !(x == y); } |
|
102 #else |
|
103 // Use a member function; Streamable concept not supported. |
|
104 bool operator==(const xor_combine& rhs) const |
|
105 { return _rng1 == rhs._rng1 && _rng2 == rhs._rng2; } |
|
106 bool operator!=(const xor_combine& rhs) const |
|
107 { return !(*this == rhs); } |
|
108 #endif |
|
109 |
|
110 private: |
|
111 base1_type _rng1; |
|
112 base2_type _rng2; |
|
113 }; |
|
114 |
|
115 #ifndef BOOST_NO_INCLASS_MEMBER_INITIALIZATION |
|
116 // A definition is required even for integral static constants |
|
117 template<class URNG1, int s1, class URNG2, int s2, |
|
118 #ifndef BOOST_NO_DEPENDENT_TYPES_IN_TEMPLATE_VALUE_PARAMETERS |
|
119 typename URNG1::result_type |
|
120 #else |
|
121 uint32_t |
|
122 #endif |
|
123 val> |
|
124 const bool xor_combine<URNG1, s1, URNG2, s2, val>::has_fixed_range; |
|
125 #endif |
|
126 |
|
127 } // namespace random |
|
128 } // namespace boost |
|
129 |
|
130 #endif // BOOST_RANDOM_XOR_COMBINE_HPP |