[857] | 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
|
---|