seed_seq.h 5.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160
  1. //===----------------------------------------------------------------------===//
  2. //
  3. // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
  4. // See https://llvm.org/LICENSE.txt for license information.
  5. // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
  6. //
  7. //===----------------------------------------------------------------------===//
  8. #ifndef _LIBCPP___RANDOM_SEED_SEQ_H
  9. #define _LIBCPP___RANDOM_SEED_SEQ_H
  10. #include <__algorithm/copy.h>
  11. #include <__algorithm/fill.h>
  12. #include <__algorithm/max.h>
  13. #include <__config>
  14. #include <__iterator/iterator_traits.h>
  15. #include <cstdint>
  16. #include <initializer_list>
  17. #include <vector>
  18. #if !defined(_LIBCPP_HAS_NO_PRAGMA_SYSTEM_HEADER)
  19. # pragma GCC system_header
  20. #endif
  21. _LIBCPP_PUSH_MACROS
  22. #include <__undef_macros>
  23. _LIBCPP_BEGIN_NAMESPACE_STD
  24. class _LIBCPP_TEMPLATE_VIS seed_seq {
  25. public:
  26. // types
  27. typedef uint32_t result_type;
  28. // constructors
  29. _LIBCPP_HIDE_FROM_ABI seed_seq() _NOEXCEPT {}
  30. #ifndef _LIBCPP_CXX03_LANG
  31. template <class _Tp, __enable_if_t<is_integral<_Tp>::value, int> = 0>
  32. _LIBCPP_HIDE_FROM_ABI seed_seq(initializer_list<_Tp> __il) {
  33. __init(__il.begin(), __il.end());
  34. }
  35. #endif // _LIBCPP_CXX03_LANG
  36. template <class _InputIterator>
  37. _LIBCPP_HIDE_FROM_ABI seed_seq(_InputIterator __first, _InputIterator __last) {
  38. static_assert(is_integral<typename iterator_traits<_InputIterator>::value_type>::value,
  39. "Mandates: iterator_traits<InputIterator>::value_type is an integer type");
  40. __init(__first, __last);
  41. }
  42. // generating functions
  43. template <class _RandomAccessIterator>
  44. _LIBCPP_HIDE_FROM_ABI void generate(_RandomAccessIterator __first, _RandomAccessIterator __last);
  45. // property functions
  46. _LIBCPP_HIDE_FROM_ABI size_t size() const _NOEXCEPT { return __v_.size(); }
  47. template <class _OutputIterator>
  48. _LIBCPP_HIDE_FROM_ABI void param(_OutputIterator __dest) const {
  49. std::copy(__v_.begin(), __v_.end(), __dest);
  50. }
  51. seed_seq(const seed_seq&) = delete;
  52. void operator=(const seed_seq&) = delete;
  53. _LIBCPP_HIDE_FROM_ABI static result_type _Tp(result_type __x) { return __x ^ (__x >> 27); }
  54. private:
  55. template <class _InputIterator>
  56. _LIBCPP_HIDE_FROM_ABI void __init(_InputIterator __first, _InputIterator __last);
  57. vector<result_type> __v_;
  58. };
  59. template <class _InputIterator>
  60. void seed_seq::__init(_InputIterator __first, _InputIterator __last) {
  61. for (_InputIterator __s = __first; __s != __last; ++__s)
  62. __v_.push_back(*__s & 0xFFFFFFFF);
  63. }
  64. template <class _RandomAccessIterator>
  65. void seed_seq::generate(_RandomAccessIterator __first, _RandomAccessIterator __last) {
  66. if (__first != __last) {
  67. std::fill(__first, __last, 0x8b8b8b8b);
  68. const size_t __n = static_cast<size_t>(__last - __first);
  69. const size_t __s = __v_.size();
  70. const size_t __t = (__n >= 623) ? 11 : (__n >= 68) ? 7 : (__n >= 39) ? 5 : (__n >= 7) ? 3 : (__n - 1) / 2;
  71. const size_t __p = (__n - __t) / 2;
  72. const size_t __q = __p + __t;
  73. const size_t __m = std::max(__s + 1, __n);
  74. // __k = 0;
  75. {
  76. result_type __r = 1664525 * _Tp(__first[0] ^ __first[__p] ^ __first[__n - 1]);
  77. __first[__p] += __r;
  78. __r += __s;
  79. __first[__q] += __r;
  80. __first[0] = __r;
  81. }
  82. // Initialize indexing terms used with if statements as an optimization to
  83. // avoid calculating modulo n on every loop iteration for each term.
  84. size_t __kmodn = 0; // __k % __n
  85. size_t __k1modn = __n - 1; // (__k - 1) % __n
  86. size_t __kpmodn = __p % __n; // (__k + __p) % __n
  87. size_t __kqmodn = __q % __n; // (__k + __q) % __n
  88. for (size_t __k = 1; __k <= __s; ++__k) {
  89. if (++__kmodn == __n)
  90. __kmodn = 0;
  91. if (++__k1modn == __n)
  92. __k1modn = 0;
  93. if (++__kpmodn == __n)
  94. __kpmodn = 0;
  95. if (++__kqmodn == __n)
  96. __kqmodn = 0;
  97. result_type __r = 1664525 * _Tp(__first[__kmodn] ^ __first[__kpmodn] ^ __first[__k1modn]);
  98. __first[__kpmodn] += __r;
  99. __r += __kmodn + __v_[__k - 1];
  100. __first[__kqmodn] += __r;
  101. __first[__kmodn] = __r;
  102. }
  103. for (size_t __k = __s + 1; __k < __m; ++__k) {
  104. if (++__kmodn == __n)
  105. __kmodn = 0;
  106. if (++__k1modn == __n)
  107. __k1modn = 0;
  108. if (++__kpmodn == __n)
  109. __kpmodn = 0;
  110. if (++__kqmodn == __n)
  111. __kqmodn = 0;
  112. result_type __r = 1664525 * _Tp(__first[__kmodn] ^ __first[__kpmodn] ^ __first[__k1modn]);
  113. __first[__kpmodn] += __r;
  114. __r += __kmodn;
  115. __first[__kqmodn] += __r;
  116. __first[__kmodn] = __r;
  117. }
  118. for (size_t __k = __m; __k < __m + __n; ++__k) {
  119. if (++__kmodn == __n)
  120. __kmodn = 0;
  121. if (++__k1modn == __n)
  122. __k1modn = 0;
  123. if (++__kpmodn == __n)
  124. __kpmodn = 0;
  125. if (++__kqmodn == __n)
  126. __kqmodn = 0;
  127. result_type __r = 1566083941 * _Tp(__first[__kmodn] + __first[__kpmodn] + __first[__k1modn]);
  128. __first[__kpmodn] ^= __r;
  129. __r -= __kmodn;
  130. __first[__kqmodn] ^= __r;
  131. __first[__kmodn] = __r;
  132. }
  133. }
  134. }
  135. _LIBCPP_END_NAMESPACE_STD
  136. _LIBCPP_POP_MACROS
  137. #endif // _LIBCPP___RANDOM_SEED_SEQ_H