Blame view

3rdparty/boost_1_81_0/boost/random/generate_canonical.hpp 2.85 KB
0b6a182c   Hu Chunming   添加无鉴权注册和注销
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
  /* boost random/generate_canonical.hpp header file
   *
   * Copyright Steven Watanabe 2011
   * Distributed under the Boost Software License, Version 1.0. (See
   * accompanying file LICENSE_1_0.txt or copy at
   * http://www.boost.org/LICENSE_1_0.txt)
   *
   * See http://www.boost.org for most recent version including documentation.
   *
   * $Id$
   *
   */
  
  #ifndef BOOST_RANDOM_GENERATE_CANONICAL_HPP
  #define BOOST_RANDOM_GENERATE_CANONICAL_HPP
  
  #include <algorithm>
  #include <boost/assert.hpp>
  #include <boost/config/no_tr1/cmath.hpp>
  #include <boost/limits.hpp>
  #include <boost/type_traits/is_integral.hpp>
  #include <boost/random/detail/signed_unsigned_tools.hpp>
  #include <boost/random/detail/generator_bits.hpp>
  
  namespace boost {
  namespace random {
  
  namespace detail {
  
  template<class RealType, std::size_t bits, class URNG>
  RealType generate_canonical_impl(URNG& g, boost::true_type /*is_integral*/)
  {
      using std::pow;
      typedef typename URNG::result_type base_result;
      std::size_t digits = std::numeric_limits<RealType>::digits;
      RealType R = RealType((g.max)()) - RealType((g.min)()) + 1;
      RealType mult = R;
      RealType limit =
          pow(RealType(2),
              RealType((std::min)(static_cast<std::size_t>(bits), digits)));
      RealType S = RealType(detail::subtract<base_result>()(g(), (g.min)()));
      while(mult < limit) {
          RealType inc = RealType(detail::subtract<base_result>()(g(), (g.min)()));
          S += inc * mult;
          mult *= R;
      }
      return S / mult;
  }
  
  template<class RealType, std::size_t bits, class URNG>
  RealType generate_canonical_impl(URNG& g, boost::false_type /*is_integral*/)
  {
      using std::pow;
      using std::floor;
      BOOST_ASSERT((g.min)() == 0);
      BOOST_ASSERT((g.max)() == 1);
      std::size_t digits = std::numeric_limits<RealType>::digits;
      std::size_t engine_bits = detail::generator_bits<URNG>::value();
      std::size_t b = (std::min)(bits, digits);
      RealType R = pow(RealType(2), RealType(engine_bits));
      RealType mult = R;
      RealType limit = pow(RealType(2), RealType(b));
      RealType S = RealType(g() - (g.min)());
      while(mult < limit) {
          RealType inc(floor((RealType(g()) - RealType((g.min)())) * R));
          S += inc * mult;
          mult *= R;
      }
      return S / mult;
  }
  
  }
  
  /**
   * Returns a value uniformly distributed in the range [0, 1)
   * with at least @c bits random bits.
   */
  template<class RealType, std::size_t bits, class URNG>
  RealType generate_canonical(URNG& g)
  {
      RealType result = detail::generate_canonical_impl<RealType, bits>(
          g, boost::random::traits::is_integral<typename URNG::result_type>());
      BOOST_ASSERT(result >= 0);
      BOOST_ASSERT(result <= 1);
      if(result == 1) {
          result -= std::numeric_limits<RealType>::epsilon() / 2;
          BOOST_ASSERT(result != 1);
      }
      return result;
  }
  
  } // namespace random
  } // namespace boost
  
  #endif // BOOST_RANDOM_GENERATE_CANONICAL_HPP