35 namespace std _GLIBCXX_VISIBILITY(default)
42 _GLIBCXX_BEGIN_NAMESPACE_VERSION
50 template<
typename _Tp, _Tp __m, _Tp __a, _Tp __c>
52 _Mod<_Tp, __m, __a, __c, false, true>::
59 static const _Tp __q = __m / __a;
60 static const _Tp __r = __m % __a;
62 _Tp __t1 = __a * (__x % __q);
63 _Tp __t2 = __r * (__x / __q);
67 __x = __m - __t2 + __t1;
72 const _Tp __d = __m - __x;
81 template<
typename _InputIterator,
typename _OutputIterator,
84 __normalize(_InputIterator __first, _InputIterator __last,
85 _OutputIterator __result,
const _Tp& __factor)
87 for (; __first != __last; ++__first, ++__result)
88 *__result = *__first / __factor;
92 _GLIBCXX_END_NAMESPACE_VERSION
95 _GLIBCXX_BEGIN_NAMESPACE_VERSION
97 template<
typename _UIntType, _UIntType __a, _UIntType __c, _UIntType __m>
99 linear_congruential_engine<_UIntType, __a, __c, __m>::multiplier;
101 template<
typename _UIntType, _UIntType __a, _UIntType __c, _UIntType __m>
103 linear_congruential_engine<_UIntType, __a, __c, __m>::increment;
105 template<
typename _UIntType, _UIntType __a, _UIntType __c, _UIntType __m>
107 linear_congruential_engine<_UIntType, __a, __c, __m>::modulus;
109 template<
typename _UIntType, _UIntType __a, _UIntType __c, _UIntType __m>
111 linear_congruential_engine<_UIntType, __a, __c, __m>::default_seed;
117 template<
typename _UIntType, _UIntType __a, _UIntType __c, _UIntType __m>
122 if ((__detail::__mod<_UIntType, __m>(__c) == 0)
123 && (__detail::__mod<_UIntType, __m>(__s) == 0))
126 _M_x = __detail::__mod<_UIntType, __m>(__s);
132 template<
typename _UIntType, _UIntType __a, _UIntType __c, _UIntType __m>
133 template<
typename _Sseq>
140 const _UIntType __k = (__k0 + 31) / 32;
141 uint_least32_t __arr[__k + 3];
142 __q.generate(__arr + 0, __arr + __k + 3);
143 _UIntType __factor = 1u;
144 _UIntType __sum = 0u;
145 for (
size_t __j = 0; __j < __k; ++__j)
147 __sum += __arr[__j + 3] * __factor;
148 __factor *= __detail::_Shift<_UIntType, 32>::__value;
153 template<
typename _UIntType, _UIntType __a, _UIntType __c, _UIntType __m,
154 typename _CharT,
typename _Traits>
156 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
158 __a, __c, __m>& __lcr)
161 typedef typename __ostream_type::ios_base __ios_base;
163 const typename __ios_base::fmtflags __flags = __os.
flags();
164 const _CharT __fill = __os.fill();
166 __os.fill(__os.widen(
' '));
175 template<
typename _UIntType, _UIntType __a, _UIntType __c, _UIntType __m,
176 typename _CharT,
typename _Traits>
179 linear_congruential_engine<_UIntType, __a, __c, __m>& __lcr)
182 typedef typename __istream_type::ios_base __ios_base;
184 const typename __ios_base::fmtflags __flags = __is.
flags();
194 template<
typename _UIntType,
195 size_t __w,
size_t __n,
size_t __m,
size_t __r,
196 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
197 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
200 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
201 __s, __b, __t, __c, __l, __f>::word_size;
203 template<
typename _UIntType,
204 size_t __w,
size_t __n,
size_t __m,
size_t __r,
205 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
206 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
209 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
210 __s, __b, __t, __c, __l, __f>::state_size;
212 template<
typename _UIntType,
213 size_t __w,
size_t __n,
size_t __m,
size_t __r,
214 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
215 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
218 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
219 __s, __b, __t, __c, __l, __f>::shift_size;
221 template<
typename _UIntType,
222 size_t __w,
size_t __n,
size_t __m,
size_t __r,
223 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
224 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
227 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
228 __s, __b, __t, __c, __l, __f>::mask_bits;
230 template<
typename _UIntType,
231 size_t __w,
size_t __n,
size_t __m,
size_t __r,
232 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
233 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
236 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
237 __s, __b, __t, __c, __l, __f>::xor_mask;
239 template<
typename _UIntType,
240 size_t __w,
size_t __n,
size_t __m,
size_t __r,
241 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
242 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
245 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
246 __s, __b, __t, __c, __l, __f>::tempering_u;
248 template<
typename _UIntType,
249 size_t __w,
size_t __n,
size_t __m,
size_t __r,
250 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
251 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
254 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
255 __s, __b, __t, __c, __l, __f>::tempering_d;
257 template<
typename _UIntType,
258 size_t __w,
size_t __n,
size_t __m,
size_t __r,
259 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
260 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
263 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
264 __s, __b, __t, __c, __l, __f>::tempering_s;
266 template<
typename _UIntType,
267 size_t __w,
size_t __n,
size_t __m,
size_t __r,
268 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
269 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
272 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
273 __s, __b, __t, __c, __l, __f>::tempering_b;
275 template<
typename _UIntType,
276 size_t __w,
size_t __n,
size_t __m,
size_t __r,
277 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
278 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
281 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
282 __s, __b, __t, __c, __l, __f>::tempering_t;
284 template<
typename _UIntType,
285 size_t __w,
size_t __n,
size_t __m,
size_t __r,
286 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
287 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
290 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
291 __s, __b, __t, __c, __l, __f>::tempering_c;
293 template<
typename _UIntType,
294 size_t __w,
size_t __n,
size_t __m,
size_t __r,
295 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
296 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
299 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
300 __s, __b, __t, __c, __l, __f>::tempering_l;
302 template<
typename _UIntType,
303 size_t __w,
size_t __n,
size_t __m,
size_t __r,
304 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
305 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
308 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
309 __s, __b, __t, __c, __l, __f>::
310 initialization_multiplier;
312 template<
typename _UIntType,
313 size_t __w,
size_t __n,
size_t __m,
size_t __r,
314 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
315 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
318 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
319 __s, __b, __t, __c, __l, __f>::default_seed;
321 template<
typename _UIntType,
322 size_t __w,
size_t __n,
size_t __m,
size_t __r,
323 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
324 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
327 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
328 __s, __b, __t, __c, __l, __f>::
329 seed(result_type __sd)
331 _M_x[0] = __detail::__mod<_UIntType,
332 __detail::_Shift<_UIntType, __w>::__value>(__sd);
334 for (
size_t __i = 1; __i < state_size; ++__i)
336 _UIntType __x = _M_x[__i - 1];
337 __x ^= __x >> (__w - 2);
339 __x += __detail::__mod<_UIntType, __n>(__i);
340 _M_x[__i] = __detail::__mod<_UIntType,
341 __detail::_Shift<_UIntType, __w>::__value>(__x);
346 template<
typename _UIntType,
347 size_t __w,
size_t __n,
size_t __m,
size_t __r,
348 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
349 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
351 template<
typename _Sseq>
353 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
354 __s, __b, __t, __c, __l, __f>::
357 const _UIntType __upper_mask = (~_UIntType()) << __r;
358 const size_t __k = (__w + 31) / 32;
359 uint_least32_t __arr[__n * __k];
360 __q.generate(__arr + 0, __arr + __n * __k);
363 for (
size_t __i = 0; __i < state_size; ++__i)
365 _UIntType __factor = 1u;
366 _UIntType __sum = 0u;
367 for (
size_t __j = 0; __j < __k; ++__j)
369 __sum += __arr[__k * __i + __j] * __factor;
370 __factor *= __detail::_Shift<_UIntType, 32>::__value;
372 _M_x[__i] = __detail::__mod<_UIntType,
373 __detail::_Shift<_UIntType, __w>::__value>(__sum);
379 if ((_M_x[0] & __upper_mask) != 0u)
382 else if (_M_x[__i] != 0u)
387 _M_x[0] = __detail::_Shift<_UIntType, __w - 1>::__value;
391 template<
typename _UIntType,
size_t __w,
392 size_t __n,
size_t __m,
size_t __r,
393 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
394 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
397 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
398 __s, __b, __t, __c, __l, __f>::
401 const _UIntType __upper_mask = (~_UIntType()) << __r;
402 const _UIntType __lower_mask = ~__upper_mask;
404 for (
size_t __k = 0; __k < (__n - __m); ++__k)
406 _UIntType __y = ((_M_x[__k] & __upper_mask)
407 | (_M_x[__k + 1] & __lower_mask));
408 _M_x[__k] = (_M_x[__k + __m] ^ (__y >> 1)
409 ^ ((__y & 0x01) ? __a : 0));
412 for (
size_t __k = (__n - __m); __k < (__n - 1); ++__k)
414 _UIntType __y = ((_M_x[__k] & __upper_mask)
415 | (_M_x[__k + 1] & __lower_mask));
416 _M_x[__k] = (_M_x[__k + (__m - __n)] ^ (__y >> 1)
417 ^ ((__y & 0x01) ? __a : 0));
420 _UIntType __y = ((_M_x[__n - 1] & __upper_mask)
421 | (_M_x[0] & __lower_mask));
422 _M_x[__n - 1] = (_M_x[__m - 1] ^ (__y >> 1)
423 ^ ((__y & 0x01) ? __a : 0));
427 template<
typename _UIntType,
size_t __w,
428 size_t __n,
size_t __m,
size_t __r,
429 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
430 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
433 mersenne_twister_engine<_UIntType, __w, __n, __m, __r, __a, __u, __d,
434 __s, __b, __t, __c, __l, __f>::
435 discard(
unsigned long long __z)
437 while (__z > state_size - _M_p)
439 __z -= state_size - _M_p;
445 template<
typename _UIntType,
size_t __w,
446 size_t __n,
size_t __m,
size_t __r,
447 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
448 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
452 __s, __b, __t, __c, __l, __f>::result_type
454 __s, __b, __t, __c, __l, __f>::
458 if (_M_p >= state_size)
462 result_type __z = _M_x[_M_p++];
463 __z ^= (__z >> __u) & __d;
464 __z ^= (__z << __s) & __b;
465 __z ^= (__z << __t) & __c;
471 template<
typename _UIntType,
size_t __w,
472 size_t __n,
size_t __m,
size_t __r,
473 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
474 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
475 _UIntType __f,
typename _CharT,
typename _Traits>
477 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
478 const mersenne_twister_engine<_UIntType, __w, __n, __m,
479 __r, __a, __u, __d, __s, __b, __t, __c, __l, __f>& __x)
482 typedef typename __ostream_type::ios_base __ios_base;
484 const typename __ios_base::fmtflags __flags = __os.
flags();
485 const _CharT __fill = __os.fill();
486 const _CharT __space = __os.widen(
' ');
490 for (
size_t __i = 0; __i < __n; ++__i)
491 __os << __x._M_x[__i] << __space;
499 template<
typename _UIntType,
size_t __w,
500 size_t __n,
size_t __m,
size_t __r,
501 _UIntType __a,
size_t __u, _UIntType __d,
size_t __s,
502 _UIntType __b,
size_t __t, _UIntType __c,
size_t __l,
503 _UIntType __f,
typename _CharT,
typename _Traits>
506 mersenne_twister_engine<_UIntType, __w, __n, __m,
507 __r, __a, __u, __d, __s, __b, __t, __c, __l, __f>& __x)
510 typedef typename __istream_type::ios_base __ios_base;
512 const typename __ios_base::fmtflags __flags = __is.
flags();
515 for (
size_t __i = 0; __i < __n; ++__i)
516 __is >> __x._M_x[__i];
524 template<
typename _UIntType,
size_t __w,
size_t __s,
size_t __r>
526 subtract_with_carry_engine<_UIntType, __w, __s, __r>::word_size;
528 template<
typename _UIntType,
size_t __w,
size_t __s,
size_t __r>
530 subtract_with_carry_engine<_UIntType, __w, __s, __r>::short_lag;
532 template<
typename _UIntType,
size_t __w,
size_t __s,
size_t __r>
534 subtract_with_carry_engine<_UIntType, __w, __s, __r>::long_lag;
536 template<
typename _UIntType,
size_t __w,
size_t __s,
size_t __r>
538 subtract_with_carry_engine<_UIntType, __w, __s, __r>::default_seed;
540 template<
typename _UIntType,
size_t __w,
size_t __s,
size_t __r>
542 subtract_with_carry_engine<_UIntType, __w, __s, __r>::
543 seed(result_type __value)
546 __lcg(__value == 0u ? default_seed : __value);
548 const size_t __n = (__w + 31) / 32;
550 for (
size_t __i = 0; __i < long_lag; ++__i)
552 _UIntType __sum = 0u;
553 _UIntType __factor = 1u;
554 for (
size_t __j = 0; __j < __n; ++__j)
556 __sum += __detail::__mod<uint_least32_t,
557 __detail::_Shift<uint_least32_t, 32>::__value>
558 (__lcg()) * __factor;
559 __factor *= __detail::_Shift<_UIntType, 32>::__value;
561 _M_x[__i] = __detail::__mod<_UIntType,
562 __detail::_Shift<_UIntType, __w>::__value>(__sum);
564 _M_carry = (_M_x[long_lag - 1] == 0) ? 1 : 0;
568 template<
typename _UIntType,
size_t __w,
size_t __s,
size_t __r>
569 template<
typename _Sseq>
571 subtract_with_carry_engine<_UIntType, __w, __s, __r>::
574 const size_t __k = (__w + 31) / 32;
575 uint_least32_t __arr[__r * __k];
576 __q.generate(__arr + 0, __arr + __r * __k);
578 for (
size_t __i = 0; __i < long_lag; ++__i)
580 _UIntType __sum = 0u;
581 _UIntType __factor = 1u;
582 for (
size_t __j = 0; __j < __k; ++__j)
584 __sum += __arr[__k * __i + __j] * __factor;
585 __factor *= __detail::_Shift<_UIntType, 32>::__value;
587 _M_x[__i] = __detail::__mod<_UIntType,
588 __detail::_Shift<_UIntType, __w>::__value>(__sum);
590 _M_carry = (_M_x[long_lag - 1] == 0) ? 1 : 0;
594 template<
typename _UIntType,
size_t __w,
size_t __s,
size_t __r>
595 typename subtract_with_carry_engine<_UIntType, __w, __s, __r>::
597 subtract_with_carry_engine<_UIntType, __w, __s, __r>::
601 long __ps = _M_p - short_lag;
609 if (_M_x[__ps] >= _M_x[_M_p] + _M_carry)
611 __xi = _M_x[__ps] - _M_x[_M_p] - _M_carry;
616 __xi = (__detail::_Shift<_UIntType, __w>::__value
617 - _M_x[_M_p] - _M_carry + _M_x[__ps]);
623 if (++_M_p >= long_lag)
629 template<
typename _UIntType,
size_t __w,
size_t __s,
size_t __r,
630 typename _CharT,
typename _Traits>
632 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
633 const subtract_with_carry_engine<_UIntType,
637 typedef typename __ostream_type::ios_base __ios_base;
639 const typename __ios_base::fmtflags __flags = __os.
flags();
640 const _CharT __fill = __os.fill();
641 const _CharT __space = __os.widen(
' ');
645 for (
size_t __i = 0; __i < __r; ++__i)
646 __os << __x._M_x[__i] << __space;
647 __os << __x._M_carry << __space << __x._M_p;
654 template<
typename _UIntType,
size_t __w,
size_t __s,
size_t __r,
655 typename _CharT,
typename _Traits>
658 subtract_with_carry_engine<_UIntType, __w, __s, __r>& __x)
661 typedef typename __istream_type::ios_base __ios_base;
663 const typename __ios_base::fmtflags __flags = __is.
flags();
666 for (
size_t __i = 0; __i < __r; ++__i)
667 __is >> __x._M_x[__i];
668 __is >> __x._M_carry;
676 template<
typename _RandomNumberEngine,
size_t __p,
size_t __r>
678 discard_block_engine<_RandomNumberEngine, __p, __r>::block_size;
680 template<
typename _RandomNumberEngine,
size_t __p,
size_t __r>
682 discard_block_engine<_RandomNumberEngine, __p, __r>::used_block;
684 template<
typename _RandomNumberEngine,
size_t __p,
size_t __r>
685 typename discard_block_engine<_RandomNumberEngine,
686 __p, __r>::result_type
690 if (_M_n >= used_block)
692 _M_b.discard(block_size - _M_n);
699 template<
typename _RandomNumberEngine,
size_t __p,
size_t __r,
700 typename _CharT,
typename _Traits>
702 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
707 typedef typename __ostream_type::ios_base __ios_base;
709 const typename __ios_base::fmtflags __flags = __os.
flags();
710 const _CharT __fill = __os.fill();
711 const _CharT __space = __os.widen(
' ');
715 __os << __x.base() << __space << __x._M_n;
722 template<
typename _RandomNumberEngine,
size_t __p,
size_t __r,
723 typename _CharT,
typename _Traits>
726 discard_block_engine<_RandomNumberEngine, __p, __r>& __x)
729 typedef typename __istream_type::ios_base __ios_base;
731 const typename __ios_base::fmtflags __flags = __is.
flags();
734 __is >> __x._M_b >> __x._M_n;
741 template<
typename _RandomNumberEngine,
size_t __w,
typename _UIntType>
742 typename independent_bits_engine<_RandomNumberEngine, __w, _UIntType>::
747 typedef typename _RandomNumberEngine::result_type _Eresult_type;
748 const _Eresult_type __r
750 ? _M_b.max() - _M_b.min() + 1 : 0);
752 const unsigned __m = __r ?
std::__lg(__r) : __edig;
759 __ctype __s0, __s1, __y0, __y1;
761 for (
size_t __i = 0; __i < 2; ++__i)
763 __n = (__w + __m - 1) / __m + __i;
764 __n0 = __n - __w % __n;
765 const unsigned __w0 = __w / __n;
771 __s0 = __ctype(1) << __w0;
779 __y0 = __s0 * (__r / __s0);
781 __y1 = __s1 * (__r / __s1);
783 if (__r - __y0 <= __y0 / __n)
791 for (
size_t __k = 0; __k < __n0; ++__k)
795 __u = _M_b() - _M_b.
min();
796 while (__y0 && __u >= __y0);
797 __sum = __s0 * __sum + (__s0 ? __u % __s0 : __u);
799 for (
size_t __k = __n0; __k < __n; ++__k)
803 __u = _M_b() - _M_b.min();
804 while (__y1 && __u >= __y1);
805 __sum = __s1 * __sum + (__s1 ? __u % __s1 : __u);
811 template<
typename _RandomNumberEngine,
size_t __k>
815 template<
typename _RandomNumberEngine,
size_t __k>
820 size_t __j = __k * ((_M_y - _M_b.min())
821 / (_M_b.max() - _M_b.min() + 1.0L));
828 template<
typename _RandomNumberEngine,
size_t __k,
829 typename _CharT,
typename _Traits>
831 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
835 typedef typename __ostream_type::ios_base __ios_base;
837 const typename __ios_base::fmtflags __flags = __os.
flags();
838 const _CharT __fill = __os.fill();
839 const _CharT __space = __os.widen(
' ');
844 for (
size_t __i = 0; __i < __k; ++__i)
845 __os << __space << __x._M_v[__i];
846 __os << __space << __x._M_y;
853 template<
typename _RandomNumberEngine,
size_t __k,
854 typename _CharT,
typename _Traits>
857 shuffle_order_engine<_RandomNumberEngine, __k>& __x)
860 typedef typename __istream_type::ios_base __ios_base;
862 const typename __ios_base::fmtflags __flags = __is.
flags();
866 for (
size_t __i = 0; __i < __k; ++__i)
867 __is >> __x._M_v[__i];
875 template<
typename _IntType>
876 template<
typename _UniformRandomNumberGenerator>
877 typename uniform_int_distribution<_IntType>::result_type
879 operator()(_UniformRandomNumberGenerator& __urng,
880 const param_type& __param)
882 typedef typename _UniformRandomNumberGenerator::result_type
884 typedef typename std::make_unsigned<result_type>::type __utype;
888 const __uctype __urngmin = __urng.min();
889 const __uctype __urngmax = __urng.max();
890 const __uctype __urngrange = __urngmax - __urngmin;
891 const __uctype __urange
892 = __uctype(__param.b()) - __uctype(__param.a());
896 if (__urngrange > __urange)
899 const __uctype __uerange = __urange + 1;
900 const __uctype __scaling = __urngrange / __uerange;
901 const __uctype __past = __uerange * __scaling;
903 __ret = __uctype(__urng()) - __urngmin;
904 while (__ret >= __past);
907 else if (__urngrange < __urange)
927 const __uctype __uerngrange = __urngrange + 1;
928 __tmp = (__uerngrange * operator()
929 (__urng, param_type(0, __urange / __uerngrange)));
930 __ret = __tmp + (__uctype(__urng()) - __urngmin);
932 while (__ret > __urange || __ret < __tmp);
935 __ret = __uctype(__urng()) - __urngmin;
937 return __ret + __param.a();
941 template<
typename _IntType>
942 template<
typename _ForwardIterator,
943 typename _UniformRandomNumberGenerator>
945 uniform_int_distribution<_IntType>::
946 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
947 _UniformRandomNumberGenerator& __urng,
948 const param_type& __param)
950 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
951 typedef typename _UniformRandomNumberGenerator::result_type
953 typedef typename std::make_unsigned<result_type>::type __utype;
954 typedef typename std::common_type<_Gresult_type, __utype>::type
957 const __uctype __urngmin = __urng.
min();
958 const __uctype __urngmax = __urng.
max();
959 const __uctype __urngrange = __urngmax - __urngmin;
960 const __uctype __urange
961 = __uctype(__param.b()) - __uctype(__param.a());
965 if (__urngrange > __urange)
967 if (__detail::_Power_of_2(__urngrange + 1)
968 && __detail::_Power_of_2(__urange + 1))
972 __ret = __uctype(__urng()) - __urngmin;
973 *__f++ = (__ret & __urange) + __param.a();
979 const __uctype __uerange = __urange + 1;
980 const __uctype __scaling = __urngrange / __uerange;
981 const __uctype __past = __uerange * __scaling;
985 __ret = __uctype(__urng()) - __urngmin;
986 while (__ret >= __past);
987 *__f++ = __ret / __scaling + __param.a();
991 else if (__urngrange < __urange)
1013 const __uctype __uerngrange = __urngrange + 1;
1014 __tmp = (__uerngrange * operator()
1015 (__urng, param_type(0, __urange / __uerngrange)));
1016 __ret = __tmp + (__uctype(__urng()) - __urngmin);
1018 while (__ret > __urange || __ret < __tmp);
1024 *__f++ = __uctype(__urng()) - __urngmin + __param.a();
1027 template<
typename _IntType,
typename _CharT,
typename _Traits>
1029 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
1033 typedef typename __ostream_type::ios_base __ios_base;
1035 const typename __ios_base::fmtflags __flags = __os.
flags();
1036 const _CharT __fill = __os.fill();
1037 const _CharT __space = __os.widen(
' ');
1041 __os << __x.a() << __space << __x.b();
1043 __os.flags(__flags);
1048 template<
typename _IntType,
typename _CharT,
typename _Traits>
1054 typedef typename __istream_type::ios_base __ios_base;
1056 const typename __ios_base::fmtflags __flags = __is.
flags();
1062 param_type(__a, __b));
1064 __is.
flags(__flags);
1069 template<
typename _RealType>
1070 template<
typename _ForwardIterator,
1071 typename _UniformRandomNumberGenerator>
1073 uniform_real_distribution<_RealType>::
1074 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
1075 _UniformRandomNumberGenerator& __urng,
1076 const param_type& __p)
1078 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
1079 __detail::_Adaptor<_UniformRandomNumberGenerator, result_type>
1081 auto __range = __p.b() - __p.a();
1083 *__f++ = __aurng() * __range + __p.a();
1086 template<typename _RealType, typename _CharT, typename _Traits>
1087 std::basic_ostream<_CharT, _Traits>&
1092 typedef typename __ostream_type::ios_base __ios_base;
1094 const typename __ios_base::fmtflags __flags = __os.
flags();
1095 const _CharT __fill = __os.fill();
1097 const _CharT __space = __os.widen(
' ');
1102 __os << __x.a() << __space << __x.b();
1104 __os.flags(__flags);
1106 __os.precision(__precision);
1110 template<
typename _RealType,
typename _CharT,
typename _Traits>
1116 typedef typename __istream_type::ios_base __ios_base;
1118 const typename __ios_base::fmtflags __flags = __is.
flags();
1124 param_type(__a, __b));
1126 __is.
flags(__flags);
1131 template<
typename _ForwardIterator,
1132 typename _UniformRandomNumberGenerator>
1134 std::bernoulli_distribution::
1135 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
1136 _UniformRandomNumberGenerator& __urng,
1137 const param_type& __p)
1139 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
1140 __detail::_Adaptor<_UniformRandomNumberGenerator,
double>
1142 auto __limit = __p.p() * (__aurng.max() - __aurng.min());
1145 *__f++ = (__aurng() - __aurng.min()) < __limit;
1148 template<typename _CharT, typename _Traits>
1149 std::basic_ostream<_CharT, _Traits>&
1154 typedef typename __ostream_type::ios_base __ios_base;
1156 const typename __ios_base::fmtflags __flags = __os.
flags();
1157 const _CharT __fill = __os.fill();
1160 __os.fill(__os.widen(
' '));
1165 __os.flags(__flags);
1167 __os.precision(__precision);
1172 template<
typename _IntType>
1173 template<
typename _UniformRandomNumberGenerator>
1174 typename geometric_distribution<_IntType>::result_type
1176 operator()(_UniformRandomNumberGenerator& __urng,
1177 const param_type& __param)
1181 const double __naf =
1184 const double __thr =
1186 __detail::_Adaptor<_UniformRandomNumberGenerator, double>
1191 __cand = std::floor(
std::log(1.0 - __aurng()) / __param._M_log_1_p);
1192 while (__cand >= __thr);
1194 return result_type(__cand + __naf);
1197 template<
typename _IntType>
1198 template<
typename _ForwardIterator,
1199 typename _UniformRandomNumberGenerator>
1201 geometric_distribution<_IntType>::
1202 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
1203 _UniformRandomNumberGenerator& __urng,
1204 const param_type& __param)
1206 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
1209 const
double __naf =
1210 (1 - std::numeric_limits<
double>::epsilon()) / 2;
1212 const
double __thr =
1213 std::numeric_limits<_IntType>::max() + __naf;
1214 __detail::_Adaptor<_UniformRandomNumberGenerator,
double>
1221 __cand = std::floor(
std::log(1.0 - __aurng())
1222 / __param._M_log_1_p);
1223 while (__cand >= __thr);
1225 *__f++ = __cand + __naf;
1229 template<
typename _IntType,
1230 typename _CharT,
typename _Traits>
1232 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
1236 typedef typename __ostream_type::ios_base __ios_base;
1238 const typename __ios_base::fmtflags __flags = __os.
flags();
1239 const _CharT __fill = __os.fill();
1242 __os.fill(__os.widen(
' '));
1247 __os.flags(__flags);
1249 __os.precision(__precision);
1253 template<
typename _IntType,
1254 typename _CharT,
typename _Traits>
1260 typedef typename __istream_type::ios_base __ios_base;
1262 const typename __ios_base::fmtflags __flags = __is.
flags();
1269 __is.
flags(__flags);
1274 template<
typename _IntType>
1275 template<
typename _UniformRandomNumberGenerator>
1276 typename negative_binomial_distribution<_IntType>::result_type
1278 operator()(_UniformRandomNumberGenerator& __urng)
1280 const double __y = _M_gd(__urng);
1284 return __poisson(__urng);
1287 template<
typename _IntType>
1288 template<
typename _UniformRandomNumberGenerator>
1289 typename negative_binomial_distribution<_IntType>::result_type
1291 operator()(_UniformRandomNumberGenerator& __urng,
1292 const param_type& __p)
1298 _M_gd(__urng, param_type(__p.k(), (1.0 - __p.p()) / __p.p()));
1301 return __poisson(__urng);
1304 template<
typename _IntType>
1305 template<
typename _ForwardIterator,
1306 typename _UniformRandomNumberGenerator>
1308 negative_binomial_distribution<_IntType>::
1309 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
1310 _UniformRandomNumberGenerator& __urng)
1312 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
1315 const double __y = _M_gd(__urng);
1319 *__f++ = __poisson(__urng);
1323 template<
typename _IntType>
1324 template<
typename _ForwardIterator,
1325 typename _UniformRandomNumberGenerator>
1327 negative_binomial_distribution<_IntType>::
1328 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
1329 _UniformRandomNumberGenerator& __urng,
1330 const param_type& __p)
1332 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
1333 typename std::gamma_distribution<result_type>::param_type
1334 __p2(__p.k(), (1.0 - __p.p()) / __p.p());
1338 const double __y = _M_gd(__urng, __p2);
1341 *__f++ = __poisson(__urng);
1345 template<
typename _IntType,
typename _CharT,
typename _Traits>
1347 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
1348 const negative_binomial_distribution<_IntType>& __x)
1351 typedef typename __ostream_type::ios_base __ios_base;
1353 const typename __ios_base::fmtflags __flags = __os.
flags();
1354 const _CharT __fill = __os.
fill();
1356 const _CharT __space = __os.
widen(
' ');
1361 __os << __x.k() << __space << __x.p()
1362 << __space << __x._M_gd;
1364 __os.
flags(__flags);
1370 template<
typename _IntType,
typename _CharT,
typename _Traits>
1373 negative_binomial_distribution<_IntType>& __x)
1376 typedef typename __istream_type::ios_base __ios_base;
1378 const typename __ios_base::fmtflags __flags = __is.
flags();
1383 __is >> __k >> __p >> __x._M_gd;
1384 __x.param(
typename negative_binomial_distribution<_IntType>::
1385 param_type(__k, __p));
1387 __is.
flags(__flags);
1392 template<
typename _IntType>
1394 poisson_distribution<_IntType>::param_type::
1397 #if _GLIBCXX_USE_C99_MATH_TR1
1400 const double __m = std::floor(_M_mean);
1402 _M_lfm = std::lgamma(__m + 1);
1405 const double __pi_4 = 0.7853981633974483096156608458198757L;
1409 const double __cx = 2 * __m + _M_d;
1414 _M_cb = 2 * __cx *
std::exp(-_M_d * _M_1cx * (1 + _M_d / 2))
1432 template<
typename _IntType>
1433 template<
typename _UniformRandomNumberGenerator>
1434 typename poisson_distribution<_IntType>::result_type
1437 const param_type& __param)
1439 __detail::_Adaptor<_UniformRandomNumberGenerator, double>
1441 #if _GLIBCXX_USE_C99_MATH_TR1
1442 if (__param.mean() >= 12)
1447 const double __naf =
1449 const double __thr =
1452 const double __m = std::floor(__param.mean());
1454 const double __spi_2 = 1.2533141373155002512078826424055226L;
1455 const double __c1 = __param._M_sm * __spi_2;
1456 const double __c2 = __param._M_c2b + __c1;
1457 const double __c3 = __c2 + 1;
1458 const double __c4 = __c3 + 1;
1460 const double __e178 = 1.0129030479320018583185514777512983L;
1461 const double __c5 = __c4 + __e178;
1462 const double __c = __param._M_cb + __c5;
1463 const double __2cx = 2 * (2 * __m + __param._M_d);
1465 bool __reject =
true;
1468 const double __u = __c * __aurng();
1469 const double __e = -
std::log(1.0 - __aurng());
1475 const double __n = _M_nd(__urng);
1476 const double __y = -
std::abs(__n) * __param._M_sm - 1;
1477 __x = std::floor(__y);
1478 __w = -__n * __n / 2;
1482 else if (__u <= __c2)
1484 const double __n = _M_nd(__urng);
1485 const double __y = 1 +
std::abs(__n) * __param._M_scx;
1486 __x = std::ceil(__y);
1487 __w = __y * (2 - __y) * __param._M_1cx;
1488 if (__x > __param._M_d)
1491 else if (__u <= __c3)
1495 else if (__u <= __c4)
1497 else if (__u <= __c5)
1501 const double __v = -
std::log(1.0 - __aurng());
1502 const double __y = __param._M_d
1503 + __v * __2cx / __param._M_d;
1504 __x = std::ceil(__y);
1505 __w = -__param._M_d * __param._M_1cx * (1 + __y / 2);
1508 __reject = (__w - __e - __x * __param._M_lm_thr
1509 > __param._M_lfm - std::lgamma(__x + __m + 1));
1511 __reject |= __x + __m >= __thr;
1521 double __prod = 1.0;
1525 __prod *= __aurng();
1528 while (__prod > __param._M_lm_thr);
1534 template<
typename _IntType>
1535 template<
typename _ForwardIterator,
1536 typename _UniformRandomNumberGenerator>
1540 _UniformRandomNumberGenerator& __urng,
1541 const param_type& __param)
1543 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
1546 *__f++ = this->operator()(__urng, __param);
1549 template<typename _IntType,
1550 typename _CharT, typename _Traits>
1556 typedef typename __ostream_type::ios_base __ios_base;
1558 const typename __ios_base::fmtflags __flags = __os.
flags();
1559 const _CharT __fill = __os.fill();
1561 const _CharT __space = __os.widen(
' ');
1566 __os << __x.mean() << __space << __x._M_nd;
1568 __os.flags(__flags);
1570 __os.precision(__precision);
1574 template<
typename _IntType,
1575 typename _CharT,
typename _Traits>
1578 poisson_distribution<_IntType>& __x)
1581 typedef typename __istream_type::ios_base __ios_base;
1583 const typename __ios_base::fmtflags __flags = __is.
flags();
1587 __is >> __mean >> __x._M_nd;
1588 __x.param(
typename poisson_distribution<_IntType>::param_type(__mean));
1590 __is.
flags(__flags);
1595 template<
typename _IntType>
1597 binomial_distribution<_IntType>::param_type::
1600 const double __p12 = _M_p <= 0.5 ? _M_p : 1.0 - _M_p;
1604 #if _GLIBCXX_USE_C99_MATH_TR1
1605 if (_M_t * __p12 >= 8)
1608 const double __np = std::floor(_M_t * __p12);
1609 const double __pa = __np / _M_t;
1610 const double __1p = 1 - __pa;
1612 const double __pi_4 = 0.7853981633974483096156608458198757L;
1613 const double __d1x =
1615 / (81 * __pi_4 * __1p)));
1616 _M_d1 = std::round(
std::max(1.0, __d1x));
1617 const double __d2x =
1619 / (__pi_4 * __pa)));
1620 _M_d2 = std::round(
std::max(1.0, __d2x));
1623 const double __spi_2 = 1.2533141373155002512078826424055226L;
1624 _M_s1 =
std::sqrt(__np * __1p) * (1 + _M_d1 / (4 * __np));
1625 _M_s2 =
std::sqrt(__np * __1p) * (1 + _M_d2 / (4 * _M_t * __1p));
1626 _M_c = 2 * _M_d1 / __np;
1627 _M_a1 =
std::exp(_M_c) * _M_s1 * __spi_2;
1628 const double __a12 = _M_a1 + _M_s2 * __spi_2;
1629 const double __s1s = _M_s1 * _M_s1;
1630 _M_a123 = __a12 + (
std::exp(_M_d1 / (_M_t * __1p))
1632 *
std::exp(-_M_d1 * _M_d1 / (2 * __s1s)));
1633 const double __s2s = _M_s2 * _M_s2;
1634 _M_s = (_M_a123 + 2 * __s2s / _M_d2
1635 *
std::exp(-_M_d2 * _M_d2 / (2 * __s2s)));
1636 _M_lf = (std::lgamma(__np + 1)
1637 + std::lgamma(_M_t - __np + 1));
1640 _M_q = -
std::log(1 - (__p12 - __pa) / __1p);
1647 template<
typename _IntType>
1648 template<
typename _UniformRandomNumberGenerator>
1649 typename binomial_distribution<_IntType>::result_type
1650 binomial_distribution<_IntType>::
1651 _M_waiting(_UniformRandomNumberGenerator& __urng, _IntType __t)
1655 __detail::_Adaptor<_UniformRandomNumberGenerator, double>
1662 const double __e = -
std::log(1.0 - __aurng());
1663 __sum += __e / (__t - __x);
1666 while (__sum <= _M_param._M_q);
1681 template<
typename _IntType>
1682 template<
typename _UniformRandomNumberGenerator>
1683 typename binomial_distribution<_IntType>::result_type
1686 const param_type& __param)
1689 const _IntType __t = __param.t();
1690 const double __p = __param.p();
1691 const double __p12 = __p <= 0.5 ? __p : 1.0 - __p;
1692 __detail::_Adaptor<_UniformRandomNumberGenerator, double>
1695 #if _GLIBCXX_USE_C99_MATH_TR1
1696 if (!__param._M_easy)
1701 const double __naf =
1703 const double __thr =
1706 const double __np = std::floor(__t * __p12);
1709 const double __spi_2 = 1.2533141373155002512078826424055226L;
1710 const double __a1 = __param._M_a1;
1711 const double __a12 = __a1 + __param._M_s2 * __spi_2;
1712 const double __a123 = __param._M_a123;
1713 const double __s1s = __param._M_s1 * __param._M_s1;
1714 const double __s2s = __param._M_s2 * __param._M_s2;
1719 const double __u = __param._M_s * __aurng();
1725 const double __n = _M_nd(__urng);
1726 const double __y = __param._M_s1 *
std::abs(__n);
1727 __reject = __y >= __param._M_d1;
1730 const double __e = -
std::log(1.0 - __aurng());
1731 __x = std::floor(__y);
1732 __v = -__e - __n * __n / 2 + __param._M_c;
1735 else if (__u <= __a12)
1737 const double __n = _M_nd(__urng);
1738 const double __y = __param._M_s2 *
std::abs(__n);
1739 __reject = __y >= __param._M_d2;
1742 const double __e = -
std::log(1.0 - __aurng());
1743 __x = std::floor(-__y);
1744 __v = -__e - __n * __n / 2;
1747 else if (__u <= __a123)
1749 const double __e1 = -
std::log(1.0 - __aurng());
1750 const double __e2 = -
std::log(1.0 - __aurng());
1752 const double __y = __param._M_d1
1753 + 2 * __s1s * __e1 / __param._M_d1;
1754 __x = std::floor(__y);
1755 __v = (-__e2 + __param._M_d1 * (1 / (__t - __np)
1756 -__y / (2 * __s1s)));
1761 const double __e1 = -
std::log(1.0 - __aurng());
1762 const double __e2 = -
std::log(1.0 - __aurng());
1764 const double __y = __param._M_d2
1765 + 2 * __s2s * __e1 / __param._M_d2;
1766 __x = std::floor(-__y);
1767 __v = -__e2 - __param._M_d2 * __y / (2 * __s2s);
1771 __reject = __reject || __x < -__np || __x > __t - __np;
1774 const double __lfx =
1775 std::lgamma(__np + __x + 1)
1776 + std::lgamma(__t - (__np + __x) + 1);
1777 __reject = __v > __param._M_lf - __lfx
1778 + __x * __param._M_lp1p;
1781 __reject |= __x + __np >= __thr;
1785 __x += __np + __naf;
1787 const _IntType __z = _M_waiting(__urng, __t - _IntType(__x));
1788 __ret = _IntType(__x) + __z;
1792 __ret = _M_waiting(__urng, __t);
1795 __ret = __t - __ret;
1799 template<
typename _IntType>
1800 template<
typename _ForwardIterator,
1801 typename _UniformRandomNumberGenerator>
1805 _UniformRandomNumberGenerator& __urng,
1806 const param_type& __param)
1808 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
1811 *__f++ = this->operator()(__urng, __param);
1814 template<typename _IntType,
1815 typename _CharT, typename _Traits>
1821 typedef typename __ostream_type::ios_base __ios_base;
1823 const typename __ios_base::fmtflags __flags = __os.
flags();
1824 const _CharT __fill = __os.fill();
1826 const _CharT __space = __os.widen(
' ');
1831 __os << __x.t() << __space << __x.p()
1832 << __space << __x._M_nd;
1834 __os.flags(__flags);
1836 __os.precision(__precision);
1840 template<
typename _IntType,
1841 typename _CharT,
typename _Traits>
1844 binomial_distribution<_IntType>& __x)
1847 typedef typename __istream_type::ios_base __ios_base;
1849 const typename __ios_base::fmtflags __flags = __is.
flags();
1854 __is >> __t >> __p >> __x._M_nd;
1855 __x.param(
typename binomial_distribution<_IntType>::
1856 param_type(__t, __p));
1858 __is.
flags(__flags);
1863 template<
typename _RealType>
1864 template<
typename _ForwardIterator,
1865 typename _UniformRandomNumberGenerator>
1869 _UniformRandomNumberGenerator& __urng,
1870 const param_type& __p)
1872 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
1873 __detail::_Adaptor<_UniformRandomNumberGenerator,
result_type>
1879 template<typename _RealType, typename _CharT, typename _Traits>
1880 std::basic_ostream<_CharT, _Traits>&
1885 typedef typename __ostream_type::ios_base __ios_base;
1887 const typename __ios_base::fmtflags __flags = __os.
flags();
1888 const _CharT __fill = __os.fill();
1891 __os.fill(__os.widen(
' '));
1894 __os << __x.lambda();
1896 __os.flags(__flags);
1898 __os.precision(__precision);
1902 template<
typename _RealType,
typename _CharT,
typename _Traits>
1908 typedef typename __istream_type::ios_base __ios_base;
1910 const typename __ios_base::fmtflags __flags = __is.
flags();
1916 param_type(__lambda));
1918 __is.
flags(__flags);
1929 template<
typename _RealType>
1930 template<
typename _UniformRandomNumberGenerator>
1931 typename normal_distribution<_RealType>::result_type
1934 const param_type& __param)
1937 __detail::_Adaptor<_UniformRandomNumberGenerator, result_type>
1940 if (_M_saved_available)
1942 _M_saved_available =
false;
1952 __r2 = __x * __x + __y * __y;
1954 while (__r2 > 1.0 || __r2 == 0.0);
1957 _M_saved = __x * __mult;
1958 _M_saved_available =
true;
1959 __ret = __y * __mult;
1962 __ret = __ret * __param.stddev() + __param.mean();
1966 template<
typename _RealType>
1967 template<
typename _ForwardIterator,
1968 typename _UniformRandomNumberGenerator>
1972 _UniformRandomNumberGenerator& __urng,
1973 const param_type& __param)
1975 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
1980 if (_M_saved_available)
1982 _M_saved_available =
false;
1983 *__f++ = _M_saved * __param.stddev() + __param.mean();
1989 __detail::_Adaptor<_UniformRandomNumberGenerator, result_type>
1992 while (__f + 1 < __t)
1999 __r2 = __x * __x + __y * __y;
2001 while (__r2 > 1.0 || __r2 == 0.0);
2004 *__f++ = __y * __mult * __param.stddev() + __param.mean();
2005 *__f++ = __x * __mult * __param.stddev() + __param.mean();
2015 __r2 = __x * __x + __y * __y;
2017 while (__r2 > 1.0 || __r2 == 0.0);
2020 _M_saved = __x * __mult;
2021 _M_saved_available =
true;
2022 *__f = __y * __mult * __param.stddev() + __param.mean();
2026 template<
typename _RealType>
2031 if (__d1._M_param == __d2._M_param
2032 && __d1._M_saved_available == __d2._M_saved_available)
2034 if (__d1._M_saved_available
2035 && __d1._M_saved == __d2._M_saved)
2037 else if(!__d1._M_saved_available)
2046 template<
typename _RealType,
typename _CharT,
typename _Traits>
2048 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
2049 const normal_distribution<_RealType>& __x)
2052 typedef typename __ostream_type::ios_base __ios_base;
2054 const typename __ios_base::fmtflags __flags = __os.
flags();
2055 const _CharT __fill = __os.
fill();
2057 const _CharT __space = __os.
widen(
' ');
2062 __os << __x.mean() << __space << __x.stddev()
2063 << __space << __x._M_saved_available;
2064 if (__x._M_saved_available)
2065 __os << __space << __x._M_saved;
2067 __os.
flags(__flags);
2073 template<
typename _RealType,
typename _CharT,
typename _Traits>
2076 normal_distribution<_RealType>& __x)
2079 typedef typename __istream_type::ios_base __ios_base;
2081 const typename __ios_base::fmtflags __flags = __is.
flags();
2084 double __mean, __stddev;
2085 __is >> __mean >> __stddev
2086 >> __x._M_saved_available;
2087 if (__x._M_saved_available)
2088 __is >> __x._M_saved;
2089 __x.param(
typename normal_distribution<_RealType>::
2090 param_type(__mean, __stddev));
2092 __is.
flags(__flags);
2097 template<
typename _RealType>
2098 template<
typename _ForwardIterator,
2099 typename _UniformRandomNumberGenerator>
2101 lognormal_distribution<_RealType>::
2102 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
2103 _UniformRandomNumberGenerator& __urng,
2104 const param_type& __p)
2106 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2108 *__f++ = std::
exp(__p.s() * _M_nd(__urng) + __p.m());
2111 template<typename _RealType, typename _CharT, typename _Traits>
2112 std::basic_ostream<_CharT, _Traits>&
2113 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
2114 const lognormal_distribution<_RealType>& __x)
2117 typedef typename __ostream_type::ios_base __ios_base;
2119 const typename __ios_base::fmtflags __flags = __os.
flags();
2120 const _CharT __fill = __os.fill();
2122 const _CharT __space = __os.widen(
' ');
2127 __os << __x.m() << __space << __x.s()
2128 << __space << __x._M_nd;
2130 __os.flags(__flags);
2132 __os.precision(__precision);
2136 template<
typename _RealType,
typename _CharT,
typename _Traits>
2139 lognormal_distribution<_RealType>& __x)
2142 typedef typename __istream_type::ios_base __ios_base;
2144 const typename __ios_base::fmtflags __flags = __is.
flags();
2148 __is >> __m >> __s >> __x._M_nd;
2149 __x.param(
typename lognormal_distribution<_RealType>::
2150 param_type(__m, __s));
2152 __is.
flags(__flags);
2156 template<
typename _RealType>
2157 template<
typename _ForwardIterator,
2158 typename _UniformRandomNumberGenerator>
2162 _UniformRandomNumberGenerator& __urng)
2164 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2166 *__f++ = 2 * _M_gd(__urng);
2169 template<typename _RealType>
2170 template<typename _ForwardIterator,
2171 typename _UniformRandomNumberGenerator>
2173 std::chi_squared_distribution<_RealType>::
2174 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
2175 _UniformRandomNumberGenerator& __urng,
2177 std::gamma_distribution<
result_type>::param_type& __p)
2179 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2181 *__f++ = 2 * _M_gd(__urng, __p);
2184 template<typename _RealType, typename _CharT, typename _Traits>
2185 std::basic_ostream<_CharT, _Traits>&
2186 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
2187 const chi_squared_distribution<_RealType>& __x)
2190 typedef typename __ostream_type::ios_base __ios_base;
2192 const typename __ios_base::fmtflags __flags = __os.
flags();
2193 const _CharT __fill = __os.fill();
2195 const _CharT __space = __os.widen(
' ');
2200 __os << __x.n() << __space << __x._M_gd;
2202 __os.flags(__flags);
2204 __os.precision(__precision);
2208 template<
typename _RealType,
typename _CharT,
typename _Traits>
2211 chi_squared_distribution<_RealType>& __x)
2214 typedef typename __istream_type::ios_base __ios_base;
2216 const typename __ios_base::fmtflags __flags = __is.
flags();
2220 __is >> __n >> __x._M_gd;
2221 __x.param(
typename chi_squared_distribution<_RealType>::
2224 __is.
flags(__flags);
2229 template<
typename _RealType>
2230 template<
typename _UniformRandomNumberGenerator>
2231 typename cauchy_distribution<_RealType>::result_type
2233 operator()(_UniformRandomNumberGenerator& __urng,
2234 const param_type& __p)
2236 __detail::_Adaptor<_UniformRandomNumberGenerator, result_type>
2243 const _RealType __pi = 3.1415926535897932384626433832795029L;
2244 return __p.a() + __p.b() *
std::tan(__pi * __u);
2247 template<
typename _RealType>
2248 template<
typename _ForwardIterator,
2249 typename _UniformRandomNumberGenerator>
2251 cauchy_distribution<_RealType>::
2252 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
2253 _UniformRandomNumberGenerator& __urng,
2254 const param_type& __p)
2256 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2257 const _RealType __pi = 3.1415926535897932384626433832795029L;
2258 __detail::_Adaptor<_UniformRandomNumberGenerator,
result_type>
2267 *__f++ = __p.a() + __p.b() *
std::tan(__pi * __u);
2271 template<
typename _RealType,
typename _CharT,
typename _Traits>
2273 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
2277 typedef typename __ostream_type::ios_base __ios_base;
2279 const typename __ios_base::fmtflags __flags = __os.
flags();
2280 const _CharT __fill = __os.fill();
2282 const _CharT __space = __os.widen(
' ');
2287 __os << __x.a() << __space << __x.b();
2289 __os.flags(__flags);
2291 __os.precision(__precision);
2295 template<
typename _RealType,
typename _CharT,
typename _Traits>
2301 typedef typename __istream_type::ios_base __ios_base;
2303 const typename __ios_base::fmtflags __flags = __is.
flags();
2309 param_type(__a, __b));
2311 __is.
flags(__flags);
2316 template<
typename _RealType>
2317 template<
typename _ForwardIterator,
2318 typename _UniformRandomNumberGenerator>
2322 _UniformRandomNumberGenerator& __urng)
2324 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2326 *__f++ = ((_M_gd_x(__urng) * n()) / (_M_gd_y(__urng) * m()));
2329 template<typename _RealType>
2330 template<typename _ForwardIterator,
2331 typename _UniformRandomNumberGenerator>
2333 std::fisher_f_distribution<_RealType>::
2334 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
2335 _UniformRandomNumberGenerator& __urng,
2336 const param_type& __p)
2338 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2339 typedef typename std::gamma_distribution<
result_type>::param_type
2341 param_type __p1(__p.m() / 2);
2342 param_type __p2(__p.n() / 2);
2344 *__f++ = ((_M_gd_x(__urng, __p1) * n())
2345 / (_M_gd_y(__urng, __p2) * m()));
2348 template<typename _RealType, typename _CharT, typename _Traits>
2349 std::basic_ostream<_CharT, _Traits>&
2350 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
2351 const fisher_f_distribution<_RealType>& __x)
2354 typedef typename __ostream_type::ios_base __ios_base;
2356 const typename __ios_base::fmtflags __flags = __os.
flags();
2357 const _CharT __fill = __os.fill();
2359 const _CharT __space = __os.widen(
' ');
2364 __os << __x.m() << __space << __x.n()
2365 << __space << __x._M_gd_x << __space << __x._M_gd_y;
2367 __os.flags(__flags);
2369 __os.precision(__precision);
2373 template<
typename _RealType,
typename _CharT,
typename _Traits>
2376 fisher_f_distribution<_RealType>& __x)
2379 typedef typename __istream_type::ios_base __ios_base;
2381 const typename __ios_base::fmtflags __flags = __is.
flags();
2385 __is >> __m >> __n >> __x._M_gd_x >> __x._M_gd_y;
2386 __x.param(
typename fisher_f_distribution<_RealType>::
2387 param_type(__m, __n));
2389 __is.
flags(__flags);
2394 template<
typename _RealType>
2395 template<
typename _ForwardIterator,
2396 typename _UniformRandomNumberGenerator>
2400 _UniformRandomNumberGenerator& __urng)
2402 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2404 *__f++ = _M_nd(__urng) * std::
sqrt(n() / _M_gd(__urng));
2407 template<typename _RealType>
2408 template<typename _ForwardIterator,
2409 typename _UniformRandomNumberGenerator>
2411 std::student_t_distribution<_RealType>::
2412 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
2413 _UniformRandomNumberGenerator& __urng,
2414 const param_type& __p)
2416 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2417 typename std::gamma_distribution<
result_type>::param_type
2418 __p2(__p.n() / 2, 2);
2420 *__f++ = _M_nd(__urng) * std::
sqrt(__p.n() / _M_gd(__urng, __p2));
2423 template<typename _RealType, typename _CharT, typename _Traits>
2424 std::basic_ostream<_CharT, _Traits>&
2425 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
2426 const student_t_distribution<_RealType>& __x)
2429 typedef typename __ostream_type::ios_base __ios_base;
2431 const typename __ios_base::fmtflags __flags = __os.
flags();
2432 const _CharT __fill = __os.fill();
2434 const _CharT __space = __os.widen(
' ');
2439 __os << __x.n() << __space << __x._M_nd << __space << __x._M_gd;
2441 __os.flags(__flags);
2443 __os.precision(__precision);
2447 template<
typename _RealType,
typename _CharT,
typename _Traits>
2450 student_t_distribution<_RealType>& __x)
2453 typedef typename __istream_type::ios_base __ios_base;
2455 const typename __ios_base::fmtflags __flags = __is.
flags();
2459 __is >> __n >> __x._M_nd >> __x._M_gd;
2460 __x.param(
typename student_t_distribution<_RealType>::param_type(__n));
2462 __is.
flags(__flags);
2467 template<
typename _RealType>
2469 gamma_distribution<_RealType>::param_type::
2472 _M_malpha = _M_alpha < 1.0 ? _M_alpha + _RealType(1.0) : _M_alpha;
2474 const _RealType __a1 = _M_malpha - _RealType(1.0) / _RealType(3.0);
2475 _M_a2 = _RealType(1.0) /
std::sqrt(_RealType(9.0) * __a1);
2483 template<
typename _RealType>
2484 template<
typename _UniformRandomNumberGenerator>
2485 typename gamma_distribution<_RealType>::result_type
2488 const param_type& __param)
2490 __detail::_Adaptor<_UniformRandomNumberGenerator, result_type>
2495 - _RealType(1.0) / _RealType(3.0));
2501 __n = _M_nd(__urng);
2506 __v = __v * __v * __v;
2509 while (__u >
result_type(1.0) - 0.331 * __n * __n * __n * __n
2510 && (
std::log(__u) > (0.5 * __n * __n + __a1
2513 if (__param.alpha() == __param._M_malpha)
2514 return __a1 * __v * __param.beta();
2522 * __a1 * __v * __param.beta());
2526 template<
typename _RealType>
2527 template<
typename _ForwardIterator,
2528 typename _UniformRandomNumberGenerator>
2532 _UniformRandomNumberGenerator& __urng,
2533 const param_type& __param)
2535 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2536 __detail::_Adaptor<_UniformRandomNumberGenerator,
result_type>
2541 - _RealType(1.0) / _RealType(3.0));
2543 if (__param.alpha() == __param._M_malpha)
2550 __n = _M_nd(__urng);
2555 __v = __v * __v * __v;
2558 while (__u >
result_type(1.0) - 0.331 * __n * __n * __n * __n
2559 && (
std::log(__u) > (0.5 * __n * __n + __a1
2562 *__f++ = __a1 * __v * __param.beta();
2571 __n = _M_nd(__urng);
2576 __v = __v * __v * __v;
2579 while (__u >
result_type(1.0) - 0.331 * __n * __n * __n * __n
2580 && (
std::log(__u) > (0.5 * __n * __n + __a1
2588 * __a1 * __v * __param.beta());
2592 template<
typename _RealType,
typename _CharT,
typename _Traits>
2594 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
2595 const gamma_distribution<_RealType>& __x)
2598 typedef typename __ostream_type::ios_base __ios_base;
2600 const typename __ios_base::fmtflags __flags = __os.
flags();
2601 const _CharT __fill = __os.
fill();
2603 const _CharT __space = __os.
widen(
' ');
2608 __os << __x.alpha() << __space << __x.beta()
2609 << __space << __x._M_nd;
2611 __os.
flags(__flags);
2617 template<
typename _RealType,
typename _CharT,
typename _Traits>
2620 gamma_distribution<_RealType>& __x)
2623 typedef typename __istream_type::ios_base __ios_base;
2625 const typename __ios_base::fmtflags __flags = __is.
flags();
2628 _RealType __alpha_val, __beta_val;
2629 __is >> __alpha_val >> __beta_val >> __x._M_nd;
2630 __x.param(
typename gamma_distribution<_RealType>::
2631 param_type(__alpha_val, __beta_val));
2633 __is.
flags(__flags);
2638 template<
typename _RealType>
2639 template<
typename _UniformRandomNumberGenerator>
2640 typename weibull_distribution<_RealType>::result_type
2642 operator()(_UniformRandomNumberGenerator& __urng,
2643 const param_type& __p)
2645 __detail::_Adaptor<_UniformRandomNumberGenerator, result_type>
2651 template<
typename _RealType>
2652 template<
typename _ForwardIterator,
2653 typename _UniformRandomNumberGenerator>
2655 weibull_distribution<_RealType>::
2656 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
2657 _UniformRandomNumberGenerator& __urng,
2658 const param_type& __p)
2660 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2661 __detail::_Adaptor<_UniformRandomNumberGenerator,
result_type>
2670 template<typename _RealType, typename _CharT, typename _Traits>
2671 std::basic_ostream<_CharT, _Traits>&
2676 typedef typename __ostream_type::ios_base __ios_base;
2678 const typename __ios_base::fmtflags __flags = __os.
flags();
2679 const _CharT __fill = __os.fill();
2681 const _CharT __space = __os.widen(
' ');
2686 __os << __x.a() << __space << __x.b();
2688 __os.flags(__flags);
2690 __os.precision(__precision);
2694 template<
typename _RealType,
typename _CharT,
typename _Traits>
2700 typedef typename __istream_type::ios_base __ios_base;
2702 const typename __ios_base::fmtflags __flags = __is.
flags();
2708 param_type(__a, __b));
2710 __is.
flags(__flags);
2715 template<
typename _RealType>
2716 template<
typename _UniformRandomNumberGenerator>
2717 typename extreme_value_distribution<_RealType>::result_type
2719 operator()(_UniformRandomNumberGenerator& __urng,
2720 const param_type& __p)
2722 __detail::_Adaptor<_UniformRandomNumberGenerator, result_type>
2728 template<
typename _RealType>
2729 template<
typename _ForwardIterator,
2730 typename _UniformRandomNumberGenerator>
2732 extreme_value_distribution<_RealType>::
2733 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
2734 _UniformRandomNumberGenerator& __urng,
2735 const param_type& __p)
2737 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2738 __detail::_Adaptor<_UniformRandomNumberGenerator,
result_type>
2746 template<typename _RealType, typename _CharT, typename _Traits>
2747 std::basic_ostream<_CharT, _Traits>&
2752 typedef typename __ostream_type::ios_base __ios_base;
2754 const typename __ios_base::fmtflags __flags = __os.
flags();
2755 const _CharT __fill = __os.fill();
2757 const _CharT __space = __os.widen(
' ');
2762 __os << __x.a() << __space << __x.b();
2764 __os.flags(__flags);
2766 __os.precision(__precision);
2770 template<
typename _RealType,
typename _CharT,
typename _Traits>
2776 typedef typename __istream_type::ios_base __ios_base;
2778 const typename __ios_base::fmtflags __flags = __is.
flags();
2784 param_type(__a, __b));
2786 __is.
flags(__flags);
2791 template<
typename _IntType>
2793 discrete_distribution<_IntType>::param_type::
2796 if (_M_prob.size() < 2)
2803 _M_prob.end(), 0.0);
2805 __detail::__normalize(_M_prob.begin(), _M_prob.end(), _M_prob.begin(),
2808 _M_cp.reserve(_M_prob.size());
2812 _M_cp[_M_cp.size() - 1] = 1.0;
2815 template<
typename _IntType>
2816 template<
typename _Func>
2817 discrete_distribution<_IntType>::param_type::
2818 param_type(
size_t __nw,
double __xmin,
double __xmax, _Func __fw)
2819 : _M_prob(), _M_cp()
2821 const size_t __n = __nw == 0 ? 1 : __nw;
2822 const double __delta = (__xmax - __xmin) / __n;
2824 _M_prob.reserve(__n);
2825 for (
size_t __k = 0; __k < __nw; ++__k)
2826 _M_prob.push_back(__fw(__xmin + __k * __delta + 0.5 * __delta));
2831 template<
typename _IntType>
2832 template<
typename _UniformRandomNumberGenerator>
2833 typename discrete_distribution<_IntType>::result_type
2834 discrete_distribution<_IntType>::
2835 operator()(_UniformRandomNumberGenerator& __urng,
2836 const param_type& __param)
2838 if (__param._M_cp.empty())
2839 return result_type(0);
2841 __detail::_Adaptor<_UniformRandomNumberGenerator, double>
2844 const double __p = __aurng();
2846 __param._M_cp.end(), __p);
2848 return __pos - __param._M_cp.begin();
2851 template<
typename _IntType>
2852 template<
typename _ForwardIterator,
2853 typename _UniformRandomNumberGenerator>
2855 discrete_distribution<_IntType>::
2856 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
2857 _UniformRandomNumberGenerator& __urng,
2858 const param_type& __param)
2860 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
2862 if (__param._M_cp.empty())
2865 *__f++ = result_type(0);
2869 __detail::_Adaptor<_UniformRandomNumberGenerator, double>
2874 const double __p = __aurng();
2876 __param._M_cp.end(), __p);
2878 *__f++ = __pos - __param._M_cp.begin();
2882 template<
typename _IntType,
typename _CharT,
typename _Traits>
2884 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
2885 const discrete_distribution<_IntType>& __x)
2888 typedef typename __ostream_type::ios_base __ios_base;
2890 const typename __ios_base::fmtflags __flags = __os.
flags();
2891 const _CharT __fill = __os.
fill();
2893 const _CharT __space = __os.
widen(
' ');
2899 __os << __prob.
size();
2900 for (
auto __dit = __prob.
begin(); __dit != __prob.
end(); ++__dit)
2901 __os << __space << *__dit;
2903 __os.
flags(__flags);
2909 template<
typename _IntType,
typename _CharT,
typename _Traits>
2912 discrete_distribution<_IntType>& __x)
2915 typedef typename __istream_type::ios_base __ios_base;
2917 const typename __ios_base::fmtflags __flags = __is.
flags();
2925 for (; __n != 0; --__n)
2932 __x.param(
typename discrete_distribution<_IntType>::
2933 param_type(__prob_vec.
begin(), __prob_vec.
end()));
2935 __is.
flags(__flags);
2940 template<
typename _RealType>
2942 piecewise_constant_distribution<_RealType>::param_type::
2945 if (_M_int.size() < 2
2946 || (_M_int.size() == 2
2947 && _M_int[0] == _RealType(0)
2948 && _M_int[1] == _RealType(1)))
2958 __detail::__normalize(_M_den.begin(), _M_den.end(), _M_den.begin(),
2961 _M_cp.reserve(_M_den.size());
2966 _M_cp[_M_cp.size() - 1] = 1.0;
2968 for (
size_t __k = 0; __k < _M_den.size(); ++__k)
2969 _M_den[__k] /= _M_int[__k + 1] - _M_int[__k];
2972 template<
typename _RealType>
2973 template<
typename _InputIteratorB,
typename _InputIteratorW>
2974 piecewise_constant_distribution<_RealType>::param_type::
2975 param_type(_InputIteratorB __bbegin,
2976 _InputIteratorB __bend,
2977 _InputIteratorW __wbegin)
2978 : _M_int(), _M_den(), _M_cp()
2980 if (__bbegin != __bend)
2984 _M_int.push_back(*__bbegin);
2986 if (__bbegin == __bend)
2989 _M_den.push_back(*__wbegin);
2997 template<
typename _RealType>
2998 template<
typename _Func>
2999 piecewise_constant_distribution<_RealType>::param_type::
3000 param_type(initializer_list<_RealType> __bl, _Func __fw)
3001 : _M_int(), _M_den(), _M_cp()
3003 _M_int.reserve(__bl.size());
3004 for (
auto __biter = __bl.begin(); __biter != __bl.end(); ++__biter)
3005 _M_int.push_back(*__biter);
3007 _M_den.reserve(_M_int.size() - 1);
3008 for (
size_t __k = 0; __k < _M_int.size() - 1; ++__k)
3009 _M_den.push_back(__fw(0.5 * (_M_int[__k + 1] + _M_int[__k])));
3014 template<
typename _RealType>
3015 template<
typename _Func>
3016 piecewise_constant_distribution<_RealType>::param_type::
3017 param_type(
size_t __nw, _RealType __xmin, _RealType __xmax, _Func __fw)
3018 : _M_int(), _M_den(), _M_cp()
3020 const size_t __n = __nw == 0 ? 1 : __nw;
3021 const _RealType __delta = (__xmax - __xmin) / __n;
3023 _M_int.reserve(__n + 1);
3024 for (
size_t __k = 0; __k <= __nw; ++__k)
3025 _M_int.push_back(__xmin + __k * __delta);
3027 _M_den.reserve(__n);
3028 for (
size_t __k = 0; __k < __nw; ++__k)
3029 _M_den.push_back(__fw(_M_int[__k] + 0.5 * __delta));
3034 template<
typename _RealType>
3035 template<
typename _UniformRandomNumberGenerator>
3036 typename piecewise_constant_distribution<_RealType>::result_type
3037 piecewise_constant_distribution<_RealType>::
3038 operator()(_UniformRandomNumberGenerator& __urng,
3039 const param_type& __param)
3041 __detail::_Adaptor<_UniformRandomNumberGenerator, double>
3044 const double __p = __aurng();
3045 if (__param._M_cp.empty())
3049 __param._M_cp.end(), __p);
3050 const size_t __i = __pos - __param._M_cp.begin();
3052 const double __pref = __i > 0 ? __param._M_cp[__i - 1] : 0.0;
3054 return __param._M_int[__i] + (__p - __pref) / __param._M_den[__i];
3057 template<
typename _RealType>
3058 template<
typename _ForwardIterator,
3059 typename _UniformRandomNumberGenerator>
3061 piecewise_constant_distribution<_RealType>::
3062 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
3063 _UniformRandomNumberGenerator& __urng,
3064 const param_type& __param)
3066 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
3067 __detail::_Adaptor<_UniformRandomNumberGenerator,
double>
3070 if (__param._M_cp.empty())
3079 const double __p = __aurng();
3082 __param._M_cp.end(), __p);
3083 const size_t __i = __pos - __param._M_cp.begin();
3085 const double __pref = __i > 0 ? __param._M_cp[__i - 1] : 0.0;
3087 *__f++ = (__param._M_int[__i]
3088 + (__p - __pref) / __param._M_den[__i]);
3092 template<
typename _RealType,
typename _CharT,
typename _Traits>
3094 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
3095 const piecewise_constant_distribution<_RealType>& __x)
3098 typedef typename __ostream_type::ios_base __ios_base;
3100 const typename __ios_base::fmtflags __flags = __os.
flags();
3101 const _CharT __fill = __os.
fill();
3103 const _CharT __space = __os.
widen(
' ');
3109 __os << __int.
size() - 1;
3111 for (
auto __xit = __int.
begin(); __xit != __int.
end(); ++__xit)
3112 __os << __space << *__xit;
3115 for (
auto __dit = __den.
begin(); __dit != __den.
end(); ++__dit)
3116 __os << __space << *__dit;
3118 __os.
flags(__flags);
3124 template<
typename _RealType,
typename _CharT,
typename _Traits>
3127 piecewise_constant_distribution<_RealType>& __x)
3130 typedef typename __istream_type::ios_base __ios_base;
3132 const typename __ios_base::fmtflags __flags = __is.
flags();
3140 for (
size_t __i = 0; __i <= __n; ++__i)
3149 for (
size_t __i = 0; __i < __n; ++__i)
3156 __x.param(
typename piecewise_constant_distribution<_RealType>::
3157 param_type(__int_vec.
begin(), __int_vec.
end(), __den_vec.
begin()));
3159 __is.
flags(__flags);
3164 template<
typename _RealType>
3166 piecewise_linear_distribution<_RealType>::param_type::
3169 if (_M_int.size() < 2
3170 || (_M_int.size() == 2
3171 && _M_int[0] == _RealType(0)
3172 && _M_int[1] == _RealType(1)
3173 && _M_den[0] == _M_den[1]))
3181 _M_cp.reserve(_M_int.size() - 1);
3182 _M_m.reserve(_M_int.size() - 1);
3183 for (
size_t __k = 0; __k < _M_int.size() - 1; ++__k)
3185 const _RealType __delta = _M_int[__k + 1] - _M_int[__k];
3186 __sum += 0.5 * (_M_den[__k + 1] + _M_den[__k]) * __delta;
3187 _M_cp.push_back(__sum);
3188 _M_m.push_back((_M_den[__k + 1] - _M_den[__k]) / __delta);
3192 __detail::__normalize(_M_den.begin(), _M_den.end(), _M_den.begin(),
3195 __detail::__normalize(_M_cp.begin(), _M_cp.end(), _M_cp.begin(), __sum);
3197 __detail::__normalize(_M_m.begin(), _M_m.end(), _M_m.begin(), __sum);
3200 _M_cp[_M_cp.size() - 1] = 1.0;
3203 template<
typename _RealType>
3204 template<
typename _InputIteratorB,
typename _InputIteratorW>
3205 piecewise_linear_distribution<_RealType>::param_type::
3206 param_type(_InputIteratorB __bbegin,
3207 _InputIteratorB __bend,
3208 _InputIteratorW __wbegin)
3209 : _M_int(), _M_den(), _M_cp(), _M_m()
3211 for (; __bbegin != __bend; ++__bbegin, ++__wbegin)
3213 _M_int.push_back(*__bbegin);
3214 _M_den.push_back(*__wbegin);
3220 template<
typename _RealType>
3221 template<
typename _Func>
3222 piecewise_linear_distribution<_RealType>::param_type::
3223 param_type(initializer_list<_RealType> __bl, _Func __fw)
3224 : _M_int(), _M_den(), _M_cp(), _M_m()
3226 _M_int.reserve(__bl.size());
3227 _M_den.reserve(__bl.size());
3228 for (
auto __biter = __bl.begin(); __biter != __bl.end(); ++__biter)
3230 _M_int.push_back(*__biter);
3231 _M_den.push_back(__fw(*__biter));
3237 template<
typename _RealType>
3238 template<
typename _Func>
3239 piecewise_linear_distribution<_RealType>::param_type::
3240 param_type(
size_t __nw, _RealType __xmin, _RealType __xmax, _Func __fw)
3241 : _M_int(), _M_den(), _M_cp(), _M_m()
3243 const size_t __n = __nw == 0 ? 1 : __nw;
3244 const _RealType __delta = (__xmax - __xmin) / __n;
3246 _M_int.reserve(__n + 1);
3247 _M_den.reserve(__n + 1);
3248 for (
size_t __k = 0; __k <= __nw; ++__k)
3250 _M_int.push_back(__xmin + __k * __delta);
3251 _M_den.push_back(__fw(_M_int[__k] + __delta));
3257 template<
typename _RealType>
3258 template<
typename _UniformRandomNumberGenerator>
3259 typename piecewise_linear_distribution<_RealType>::result_type
3260 piecewise_linear_distribution<_RealType>::
3261 operator()(_UniformRandomNumberGenerator& __urng,
3262 const param_type& __param)
3264 __detail::_Adaptor<_UniformRandomNumberGenerator, double>
3267 const double __p = __aurng();
3268 if (__param._M_cp.empty())
3272 __param._M_cp.end(), __p);
3273 const size_t __i = __pos - __param._M_cp.begin();
3275 const double __pref = __i > 0 ? __param._M_cp[__i - 1] : 0.0;
3277 const double __a = 0.5 * __param._M_m[__i];
3278 const double __b = __param._M_den[__i];
3279 const double __cm = __p - __pref;
3281 _RealType __x = __param._M_int[__i];
3286 const double __d = __b * __b + 4.0 * __a * __cm;
3287 __x += 0.5 * (
std::sqrt(__d) - __b) / __a;
3293 template<
typename _RealType>
3294 template<
typename _ForwardIterator,
3295 typename _UniformRandomNumberGenerator>
3297 piecewise_linear_distribution<_RealType>::
3298 __generate_impl(_ForwardIterator __f, _ForwardIterator __t,
3299 _UniformRandomNumberGenerator& __urng,
3300 const param_type& __param)
3302 __glibcxx_function_requires(_ForwardIteratorConcept<_ForwardIterator>)
3305 *__f++ = this->operator()(__urng, __param);
3308 template<typename _RealType, typename _CharT, typename _Traits>
3309 std::basic_ostream<_CharT, _Traits>&
3310 operator<<(std::basic_ostream<_CharT, _Traits>& __os,
3311 const piecewise_linear_distribution<_RealType>& __x)
3314 typedef typename __ostream_type::ios_base __ios_base;
3316 const typename __ios_base::fmtflags __flags = __os.
flags();
3317 const _CharT __fill = __os.fill();
3319 const _CharT __space = __os.widen(
' ');
3325 __os << __int.
size() - 1;
3327 for (
auto __xit = __int.
begin(); __xit != __int.
end(); ++__xit)
3328 __os << __space << *__xit;
3331 for (
auto __dit = __den.
begin(); __dit != __den.
end(); ++__dit)
3332 __os << __space << *__dit;
3334 __os.flags(__flags);
3336 __os.precision(__precision);
3340 template<
typename _RealType,
typename _CharT,
typename _Traits>
3343 piecewise_linear_distribution<_RealType>& __x)
3346 typedef typename __istream_type::ios_base __ios_base;
3348 const typename __ios_base::fmtflags __flags = __is.
flags();
3356 for (
size_t __i = 0; __i <= __n; ++__i)
3365 for (
size_t __i = 0; __i <= __n; ++__i)
3372 __x.param(
typename piecewise_linear_distribution<_RealType>::
3373 param_type(__int_vec.
begin(), __int_vec.
end(), __den_vec.
begin()));
3375 __is.
flags(__flags);
3380 template<
typename _IntType>
3383 for (
auto __iter = __il.begin(); __iter != __il.end(); ++__iter)
3384 _M_v.push_back(__detail::__mod<result_type,
3385 __detail::_Shift<result_type, 32>::__value>(*__iter));
3388 template<
typename _InputIterator>
3389 seed_seq::seed_seq(_InputIterator __begin, _InputIterator __end)
3391 for (_InputIterator __iter = __begin; __iter != __end; ++__iter)
3392 _M_v.push_back(__detail::__mod<result_type,
3393 __detail::_Shift<result_type, 32>::__value>(*__iter));
3396 template<
typename _RandomAccessIterator>
3398 seed_seq::generate(_RandomAccessIterator __begin,
3399 _RandomAccessIterator __end)
3401 typedef typename iterator_traits<_RandomAccessIterator>::value_type
3404 if (__begin == __end)
3407 std::fill(__begin, __end, _Type(0x8b8b8b8bu));
3409 const size_t __n = __end - __begin;
3410 const size_t __s = _M_v.size();
3411 const size_t __t = (__n >= 623) ? 11
3416 const size_t __p = (__n - __t) / 2;
3417 const size_t __q = __p + __t;
3418 const size_t __m =
std::max(
size_t(__s + 1), __n);
3420 for (
size_t __k = 0; __k < __m; ++__k)
3422 _Type __arg = (__begin[__k % __n]
3423 ^ __begin[(__k + __p) % __n]
3424 ^ __begin[(__k - 1) % __n]);
3425 _Type __r1 = __arg ^ (__arg >> 27);
3426 __r1 = __detail::__mod<_Type,
3427 __detail::_Shift<_Type, 32>::__value>(1664525u * __r1);
3431 else if (__k <= __s)
3432 __r2 += __k % __n + _M_v[__k - 1];
3435 __r2 = __detail::__mod<_Type,
3436 __detail::_Shift<_Type, 32>::__value>(__r2);
3437 __begin[(__k + __p) % __n] += __r1;
3438 __begin[(__k + __q) % __n] += __r2;
3439 __begin[__k % __n] = __r2;
3442 for (
size_t __k = __m; __k < __m + __n; ++__k)
3444 _Type __arg = (__begin[__k % __n]
3445 + __begin[(__k + __p) % __n]
3446 + __begin[(__k - 1) % __n]);
3447 _Type __r3 = __arg ^ (__arg >> 27);
3448 __r3 = __detail::__mod<_Type,
3449 __detail::_Shift<_Type, 32>::__value>(1566083941u * __r3);
3450 _Type __r4 = __r3 - __k % __n;
3451 __r4 = __detail::__mod<_Type,
3452 __detail::_Shift<_Type, 32>::__value>(__r4);
3453 __begin[(__k + __p) % __n] ^= __r3;
3454 __begin[(__k + __q) % __n] ^= __r4;
3455 __begin[__k % __n] = __r4;
3459 template<
typename _RealType,
size_t __bits,
3460 typename _UniformRandomNumberGenerator>
3467 const long double __r =
static_cast<long double>(__urng.max())
3468 - static_cast<long double>(__urng.min()) + 1.0L;
3470 size_t __k = std::max<size_t>(1UL, (__b + __log2r - 1UL) / __log2r);
3471 _RealType __sum = _RealType(0);
3472 _RealType __tmp = _RealType(1);
3473 for (; __k != 0; --__k)
3475 __sum += _RealType(__urng() - __urng.min()) * __tmp;
3478 return __sum / __tmp;
3481 _GLIBCXX_END_NAMESPACE_VERSION