* include/ext/random: New file.
* include/ext/random.tcc: New file.
* include/Makefile.am (ext_headers): Add random and random.tcc.
* testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
operators/inequal.cc: New file.
* testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
operators/equal.cc: New file.
* testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
operators/serialize.cc: New file.
* testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
cons/seed2.cc: New file.
* testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
cons/default.cc: New file.
* testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
cons/seed_seq.cc: New file.
* testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
cons/copy.cc: New file.
* testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
cons/seed1.cc: New file.
From-SVN: r190783
+2012-08-29 Ulrich Drepper <drepper@gmail.com>
+
+ Add simd_fast_mersenne_twister_engine class.
+ * include/ext/random: New file.
+ * include/ext/random.tcc: New file.
+ * include/Makefile.am (ext_headers): Add random and random.tcc.
+ * testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
+ operators/inequal.cc: New file.
+ * testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
+ operators/equal.cc: New file.
+ * testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
+ operators/serialize.cc: New file.
+ * testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
+ cons/seed2.cc: New file.
+ * testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
+ cons/default.cc: New file.
+ * testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
+ cons/seed_seq.cc: New file.
+ * testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
+ cons/copy.cc: New file.
+ * testsuite/26_numerics/random/simd_fast_mersenne_twister_engine/
+ cons/seed1.cc: New file.
+
2012-08-28 Benjamin Kosnik <bkoz@redhat.com>
PR libstdc++/54102, part 2
${ext_srcdir}/pointer.h \
${ext_srcdir}/pool_allocator.h \
${ext_srcdir}/rb_tree \
+ ${ext_srcdir}/random \
+ ${ext_srcdir}/random.tcc \
${ext_srcdir}/rope \
${ext_srcdir}/ropeimpl.h \
${ext_srcdir}/slist \
visibility=`cat stamp-visibility` ;\
externtemplate=`cat stamp-extern-template` ;\
ldbl_compat='s,g,g,' ;\
- grep "^[ ]*#[ ]*define[ ][ ]*_GLIBCXX_LONG_DOUBLE_COMPAT[ ][ ]*1[ ]*$$" \
+ grep "^[ ]*#[ ]*define[ ][ ]*_GLIBCXX_LONG_DOUBLE_COMPAT[ ][ ]*1[ ]*$$" \
${CONFIG_HEADER} > /dev/null 2>&1 \
&& ldbl_compat='s,^#undef _GLIBCXX_LONG_DOUBLE_COMPAT$$,#define _GLIBCXX_LONG_DOUBLE_COMPAT 1,' ;\
sed -e "s,define __GLIBCXX__,define __GLIBCXX__ $$date," \
-e 's/VERSION/_GLIBCXX_VERSION/g' \
-e 's/WORDS_/_GLIBCXX_WORDS_/g' \
-e 's/ICONV_CONST/_GLIBCXX_ICONV_CONST/g' \
- -e '/[ ]_GLIBCXX_LONG_DOUBLE_COMPAT[ ]/d' \
+ -e '/[ ]_GLIBCXX_LONG_DOUBLE_COMPAT[ ]/d' \
< ${CONFIG_HEADER} >> $@ ;\
echo "" >> $@ ;\
echo "#endif // _GLIBCXX_CXX_CONFIG_H" >> $@
-# Makefile.in generated by automake 1.11.1 from Makefile.am.
+# Makefile.in generated by automake 1.11.3 from Makefile.am.
# @configure_input@
# Copyright (C) 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002,
-# 2003, 2004, 2005, 2006, 2007, 2008, 2009 Free Software Foundation,
-# Inc.
+# 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011 Free Software
+# Foundation, Inc.
# This Makefile.in is free software; the Free Software Foundation
# gives unlimited permission to copy and/or distribute it,
# with or without modifications, as long as this notice is preserved.
${ext_srcdir}/pointer.h \
${ext_srcdir}/pool_allocator.h \
${ext_srcdir}/rb_tree \
+ ${ext_srcdir}/random \
+ ${ext_srcdir}/random.tcc \
${ext_srcdir}/rope \
${ext_srcdir}/ropeimpl.h \
${ext_srcdir}/slist \
echo ' cd $(top_builddir) && $(SHELL) ./config.status $(subdir)/$@ $(am__depfiles_maybe)'; \
cd $(top_builddir) && $(SHELL) ./config.status $(subdir)/$@ $(am__depfiles_maybe);; \
esac;
+$(top_srcdir)/fragment.am:
$(top_builddir)/config.status: $(top_srcdir)/configure $(CONFIG_STATUS_DEPENDENCIES)
cd $(top_builddir) && $(MAKE) $(AM_MAKEFLAGS) am--refresh
installcheck: installcheck-am
install-strip:
- $(MAKE) $(AM_MAKEFLAGS) INSTALL_PROGRAM="$(INSTALL_STRIP_PROGRAM)" \
- install_sh_PROGRAM="$(INSTALL_STRIP_PROGRAM)" INSTALL_STRIP_FLAG=-s \
- `test -z '$(STRIP)' || \
- echo "INSTALL_PROGRAM_ENV=STRIPPROG='$(STRIP)'"` install
+ if test -z '$(STRIP)'; then \
+ $(MAKE) $(AM_MAKEFLAGS) INSTALL_PROGRAM="$(INSTALL_STRIP_PROGRAM)" \
+ install_sh_PROGRAM="$(INSTALL_STRIP_PROGRAM)" INSTALL_STRIP_FLAG=-s \
+ install; \
+ else \
+ $(MAKE) $(AM_MAKEFLAGS) INSTALL_PROGRAM="$(INSTALL_STRIP_PROGRAM)" \
+ install_sh_PROGRAM="$(INSTALL_STRIP_PROGRAM)" INSTALL_STRIP_FLAG=-s \
+ "INSTALL_PROGRAM_ENV=STRIPPROG='$(STRIP)'" install; \
+ fi
mostlyclean-generic:
clean-generic:
visibility=`cat stamp-visibility` ;\
externtemplate=`cat stamp-extern-template` ;\
ldbl_compat='s,g,g,' ;\
- grep "^[ ]*#[ ]*define[ ][ ]*_GLIBCXX_LONG_DOUBLE_COMPAT[ ][ ]*1[ ]*$$" \
+ grep "^[ ]*#[ ]*define[ ][ ]*_GLIBCXX_LONG_DOUBLE_COMPAT[ ][ ]*1[ ]*$$" \
${CONFIG_HEADER} > /dev/null 2>&1 \
&& ldbl_compat='s,^#undef _GLIBCXX_LONG_DOUBLE_COMPAT$$,#define _GLIBCXX_LONG_DOUBLE_COMPAT 1,' ;\
sed -e "s,define __GLIBCXX__,define __GLIBCXX__ $$date," \
-e 's/VERSION/_GLIBCXX_VERSION/g' \
-e 's/WORDS_/_GLIBCXX_WORDS_/g' \
-e 's/ICONV_CONST/_GLIBCXX_ICONV_CONST/g' \
- -e '/[ ]_GLIBCXX_LONG_DOUBLE_COMPAT[ ]/d' \
+ -e '/[ ]_GLIBCXX_LONG_DOUBLE_COMPAT[ ]/d' \
< ${CONFIG_HEADER} >> $@ ;\
echo "" >> $@ ;\
echo "#endif // _GLIBCXX_CXX_CONFIG_H" >> $@
--- /dev/null
+// Random number extensions -*- C++ -*-
+
+// Copyright (C) 2012 Free Software Foundation, Inc.
+//
+// This file is part of the GNU ISO C++ Library. This library is free
+// software; you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the
+// Free Software Foundation; either version 3, or (at your option)
+// any later version.
+
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+
+// Under Section 7 of GPL version 3, you are granted additional
+// permissions described in the GCC Runtime Library Exception, version
+// 3.1, as published by the Free Software Foundation.
+
+// You should have received a copy of the GNU General Public License and
+// a copy of the GCC Runtime Library Exception along with this program;
+// see the files COPYING3 and COPYING.RUNTIME respectively. If not, see
+// <http://www.gnu.org/licenses/>.
+
+/** @file ext/random
+ * This file is a GNU extension to the Standard C++ Library.
+ */
+
+#ifndef _EXT_RANDOM
+#define _EXT_RANDOM 1
+
+#pragma GCC system_header
+
+#include <random>
+#ifdef __SSE2__
+# include <x86intrin.h>
+#endif
+
+
+namespace __gnu_cxx _GLIBCXX_VISIBILITY(default)
+{
+_GLIBCXX_BEGIN_NAMESPACE_VERSION
+
+ /* Mersenne twister implementation optimized for vector operations.
+ *
+ * Reference: http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/SFMT/
+ */
+ template<typename _UIntType, size_t __m,
+ size_t __pos1, size_t __sl1, size_t __sl2,
+ size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2,
+ uint32_t __msk3, uint32_t __msk4,
+ uint32_t __parity1, uint32_t __parity2,
+ uint32_t __parity3, uint32_t __parity4>
+ class simd_fast_mersenne_twister_engine
+ {
+ static_assert(std::is_unsigned<_UIntType>::value, "template argument "
+ "substituting _UIntType not an unsigned integral type");
+ static_assert(__sr1 < 32, "first right shift too large");
+ static_assert(__sr2 < 16, "second right shift too large");
+ static_assert(__sl1 < 32, "first left shift too large");
+ static_assert(__sl2 < 16, "second left shift too large");
+
+ public:
+ typedef _UIntType result_type;
+
+ private:
+ static constexpr size_t m_w = sizeof(result_type) * 8;
+ static constexpr size_t _M_nstate = __m / 128 + 1;
+ static constexpr size_t _M_nstate32 = _M_nstate * 4;
+
+ static_assert(std::is_unsigned<_UIntType>::value, "template argument "
+ "substituting _UIntType not an unsigned integral type");
+ static_assert(__pos1 < _M_nstate, "POS1 not smaller than state size");
+ static_assert(16 % sizeof(_UIntType) == 0,
+ "UIntType size must divide 16");
+
+ public:
+ static constexpr size_t state_size = _M_nstate * (16
+ / sizeof(result_type));
+ static constexpr result_type default_seed = 5489u;
+
+ // constructors and member function
+ explicit
+ simd_fast_mersenne_twister_engine(result_type __sd = default_seed)
+ { seed(__sd); }
+
+ template<typename _Sseq, typename = typename
+ std::enable_if<!std::is_same<_Sseq, simd_fast_mersenne_twister_engine>::value>
+ ::type>
+ explicit
+ simd_fast_mersenne_twister_engine(_Sseq& __q)
+ { seed(__q); }
+
+ void
+ seed(result_type __sd = default_seed);
+
+ template<typename _Sseq>
+ typename std::enable_if<std::is_class<_Sseq>::value>::type
+ seed(_Sseq& __q);
+
+ static constexpr result_type
+ min()
+ { return 0; };
+
+ static constexpr result_type
+ max()
+ { return std::numeric_limits<result_type>::max(); }
+
+ void
+ discard(unsigned long long __z);
+
+ result_type
+ operator()()
+ {
+ if (__builtin_expect(_M_pos >= state_size, 0))
+ _M_gen_rand();
+
+ return _M_stateT[_M_pos++];
+ }
+
+#ifdef __SSE2__
+ friend bool
+ operator==(const simd_fast_mersenne_twister_engine& __lhs,
+ const simd_fast_mersenne_twister_engine& __rhs)
+ { __m128i __res = _mm_cmpeq_epi8(__lhs._M_state[0], __rhs._M_state[0]);
+ for (size_t __i = 1; __i < __lhs._M_nstate; ++__i)
+ __res = _mm_and_si128(__res, _mm_cmpeq_epi8(__lhs._M_state[__i],
+ __rhs._M_state[__i]));
+ return (_mm_movemask_epi8(__res) == 0xffff
+ && __lhs._M_pos == __rhs._M_pos); }
+#else
+ friend bool
+ operator==(const simd_fast_mersenne_twister_engine& __lhs,
+ const simd_fast_mersenne_twister_engine& __rhs)
+ { return (std::equal(__lhs._M_stateT, __lhs._M_stateT + state_size,
+ __rhs._M_stateT)
+ && __lhs._M_pos == __rhs._M_pos); }
+#endif
+
+ template<typename _UIntType_2, size_t __m_2,
+ size_t __pos1_2, size_t __sl1_2, size_t __sl2_2,
+ size_t __sr1_2, size_t __sr2_2,
+ uint32_t __msk1_2, uint32_t __msk2_2,
+ uint32_t __msk3_2, uint32_t __msk4_2,
+ uint32_t __parity1_2, uint32_t __parity2_2,
+ uint32_t __parity3_2, uint32_t __parity4_2,
+ typename _CharT, typename _Traits>
+ friend std::basic_ostream<_CharT, _Traits>&
+ operator<<(std::basic_ostream<_CharT, _Traits>& __os,
+ const __gnu_cxx::simd_fast_mersenne_twister_engine<_UIntType_2,
+ __m_2, __pos1_2, __sl1_2, __sl2_2, __sr1_2, __sr2_2,
+ __msk1_2, __msk2_2, __msk3_2, __msk4_2,
+ __parity1_2, __parity2_2, __parity3_2, __parity4_2>& __x);
+
+ template<typename _UIntType_2, size_t __m_2,
+ size_t __pos1_2, size_t __sl1_2, size_t __sl2_2,
+ size_t __sr1_2, size_t __sr2_2,
+ uint32_t __msk1_2, uint32_t __msk2_2,
+ uint32_t __msk3_2, uint32_t __msk4_2,
+ uint32_t __parity1_2, uint32_t __parity2_2,
+ uint32_t __parity3_2, uint32_t __parity4_2,
+ typename _CharT, typename _Traits>
+ friend std::basic_istream<_CharT, _Traits>&
+ operator>>(std::basic_istream<_CharT, _Traits>& __is,
+ __gnu_cxx::simd_fast_mersenne_twister_engine<_UIntType_2,
+ __m_2, __pos1_2, __sl1_2, __sl2_2, __sr1_2, __sr2_2,
+ __msk1_2, __msk2_2, __msk3_2, __msk4_2,
+ __parity1_2, __parity2_2, __parity3_2, __parity4_2>& __x);
+
+ private:
+ union
+ {
+#ifdef __SSE2__
+ __m128i _M_state[_M_nstate];
+#endif
+ uint32_t _M_state32[_M_nstate32];
+ result_type _M_stateT[state_size];
+ } __attribute__ ((__aligned__ (16)));
+ size_t _M_pos;
+
+ void _M_gen_rand(void);
+ void _M_period_certification();
+ };
+
+
+ template<typename _UIntType, size_t __m,
+ size_t __pos1, size_t __sl1, size_t __sl2,
+ size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2,
+ uint32_t __msk3, uint32_t __msk4,
+ uint32_t __parity1, uint32_t __parity2,
+ uint32_t __parity3, uint32_t __parity4>
+ inline bool
+ operator!=(const __gnu_cxx::simd_fast_mersenne_twister_engine<_UIntType,
+ __m, __pos1, __sl1, __sl2, __sr1, __sr2, __msk1, __msk2, __msk3,
+ __msk4, __parity1, __parity2, __parity3, __parity4>& __lhs,
+ const __gnu_cxx::simd_fast_mersenne_twister_engine<_UIntType,
+ __m, __pos1, __sl1, __sl2, __sr1, __sr2, __msk1, __msk2, __msk3,
+ __msk4, __parity1, __parity2, __parity3, __parity4>& __rhs)
+ { return !(__lhs == __rhs); }
+
+
+ /* Definitions for the SIMD-oriented Fast Mersenne Twister as defined
+ * in the C implementation by Daito and Matsumoto, as both a 32-bit
+ * and 64-bit version.
+ */
+ typedef simd_fast_mersenne_twister_engine<uint32_t, 607, 2,
+ 15, 3, 13, 3,
+ 0xfdff37ffU, 0xef7f3f7dU,
+ 0xff777b7dU, 0x7ff7fb2fU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x5986f054U>
+ sfmt607;
+
+ typedef simd_fast_mersenne_twister_engine<uint64_t, 607, 2,
+ 15, 3, 13, 3,
+ 0xfdff37ffU, 0xef7f3f7dU,
+ 0xff777b7dU, 0x7ff7fb2fU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x5986f054U>
+ sfmt607_64;
+
+
+ typedef simd_fast_mersenne_twister_engine<uint32_t, 1279, 7,
+ 14, 3, 5, 1,
+ 0xf7fefffdU, 0x7fefcfffU,
+ 0xaff3ef3fU, 0xb5ffff7fU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x20000000U>
+ sfmt1279;
+
+ typedef simd_fast_mersenne_twister_engine<uint64_t, 1279, 7,
+ 14, 3, 5, 1,
+ 0xf7fefffdU, 0x7fefcfffU,
+ 0xaff3ef3fU, 0xb5ffff7fU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x20000000U>
+ sfmt1279_64;
+
+
+ typedef simd_fast_mersenne_twister_engine<uint32_t, 2281, 12,
+ 19, 1, 5, 1,
+ 0xbff7ffbfU, 0xfdfffffeU,
+ 0xf7ffef7fU, 0xf2f7cbbfU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x41dfa600U>
+ sfmt2281;
+
+ typedef simd_fast_mersenne_twister_engine<uint64_t, 2281, 12,
+ 19, 1, 5, 1,
+ 0xbff7ffbfU, 0xfdfffffeU,
+ 0xf7ffef7fU, 0xf2f7cbbfU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x41dfa600U>
+ sfmt2281_64;
+
+
+ typedef simd_fast_mersenne_twister_engine<uint32_t, 4253, 17,
+ 20, 1, 7, 1,
+ 0x9f7bffffU, 0x9fffff5fU,
+ 0x3efffffbU, 0xfffff7bbU,
+ 0xa8000001U, 0xaf5390a3U,
+ 0xb740b3f8U, 0x6c11486dU>
+ sfmt4253;
+
+ typedef simd_fast_mersenne_twister_engine<uint64_t, 4253, 17,
+ 20, 1, 7, 1,
+ 0x9f7bffffU, 0x9fffff5fU,
+ 0x3efffffbU, 0xfffff7bbU,
+ 0xa8000001U, 0xaf5390a3U,
+ 0xb740b3f8U, 0x6c11486dU>
+ sfmt4253_64;
+
+
+ typedef simd_fast_mersenne_twister_engine<uint32_t, 11213, 68,
+ 14, 3, 7, 3,
+ 0xeffff7fbU, 0xffffffefU,
+ 0xdfdfbfffU, 0x7fffdbfdU,
+ 0x00000001U, 0x00000000U,
+ 0xe8148000U, 0xd0c7afa3U>
+ sfmt11213;
+
+ typedef simd_fast_mersenne_twister_engine<uint64_t, 11213, 68,
+ 14, 3, 7, 3,
+ 0xeffff7fbU, 0xffffffefU,
+ 0xdfdfbfffU, 0x7fffdbfdU,
+ 0x00000001U, 0x00000000U,
+ 0xe8148000U, 0xd0c7afa3U>
+ sfmt11213_64;
+
+
+ typedef simd_fast_mersenne_twister_engine<uint32_t, 19937, 122,
+ 18, 1, 11, 1,
+ 0xdfffffefU, 0xddfecb7fU,
+ 0xbffaffffU, 0xbffffff6U,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x13c9e684U>
+ sfmt19937;
+
+ typedef simd_fast_mersenne_twister_engine<uint64_t, 19937, 122,
+ 18, 1, 11, 1,
+ 0xdfffffefU, 0xddfecb7fU,
+ 0xbffaffffU, 0xbffffff6U,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x13c9e684U>
+ sfmt19937_64;
+
+
+ typedef simd_fast_mersenne_twister_engine<uint32_t, 44497, 330,
+ 5, 3, 9, 3,
+ 0xeffffffbU, 0xdfbebfffU,
+ 0xbfbf7befU, 0x9ffd7bffU,
+ 0x00000001U, 0x00000000U,
+ 0xa3ac4000U, 0xecc1327aU>
+ sfmt44497;
+
+ typedef simd_fast_mersenne_twister_engine<uint64_t, 44497, 330,
+ 5, 3, 9, 3,
+ 0xeffffffbU, 0xdfbebfffU,
+ 0xbfbf7befU, 0x9ffd7bffU,
+ 0x00000001U, 0x00000000U,
+ 0xa3ac4000U, 0xecc1327aU>
+ sfmt44497_64;
+
+
+ typedef simd_fast_mersenne_twister_engine<uint32_t, 86243, 366,
+ 6, 7, 19, 1,
+ 0xfdbffbffU, 0xbff7ff3fU,
+ 0xfd77efffU, 0xbf9ff3ffU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0xe9528d85U>
+ sfmt86243;
+
+ typedef simd_fast_mersenne_twister_engine<uint64_t, 86243, 366,
+ 6, 7, 19, 1,
+ 0xfdbffbffU, 0xbff7ff3fU,
+ 0xfd77efffU, 0xbf9ff3ffU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0xe9528d85U>
+ sfmt86243_64;
+
+
+ typedef simd_fast_mersenne_twister_engine<uint32_t, 132049, 110,
+ 19, 1, 21, 1,
+ 0xffffbb5fU, 0xfb6ebf95U,
+ 0xfffefffaU, 0xcff77fffU,
+ 0x00000001U, 0x00000000U,
+ 0xcb520000U, 0xc7e91c7dU>
+ sfmt132049;
+
+ typedef simd_fast_mersenne_twister_engine<uint64_t, 132049, 110,
+ 19, 1, 21, 1,
+ 0xffffbb5fU, 0xfb6ebf95U,
+ 0xfffefffaU, 0xcff77fffU,
+ 0x00000001U, 0x00000000U,
+ 0xcb520000U, 0xc7e91c7dU>
+ sfmt132049_64;
+
+
+ typedef simd_fast_mersenne_twister_engine<uint32_t, 216091, 627,
+ 11, 3, 10, 1,
+ 0xbff7bff7U, 0xbfffffffU,
+ 0xbffffa7fU, 0xffddfbfbU,
+ 0xf8000001U, 0x89e80709U,
+ 0x3bd2b64bU, 0x0c64b1e4U>
+ sfmt216091;
+
+ typedef simd_fast_mersenne_twister_engine<uint64_t, 216091, 627,
+ 11, 3, 10, 1,
+ 0xbff7bff7U, 0xbfffffffU,
+ 0xbffffa7fU, 0xffddfbfbU,
+ 0xf8000001U, 0x89e80709U,
+ 0x3bd2b64bU, 0x0c64b1e4U>
+ sfmt216091_64;
+
+_GLIBCXX_END_NAMESPACE_VERSION
+} // namespace std
+
+#include "random.tcc"
+
+#endif /* _EXT_RANDOM */
--- /dev/null
+// Random number extensions -*- C++ -*-
+
+// Copyright (C) 2012 Free Software Foundation, Inc.
+//
+// This file is part of the GNU ISO C++ Library. This library is free
+// software; you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the
+// Free Software Foundation; either version 3, or (at your option)
+// any later version.
+
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+
+// Under Section 7 of GPL version 3, you are granted additional
+// permissions described in the GCC Runtime Library Exception, version
+// 3.1, as published by the Free Software Foundation.
+
+// You should have received a copy of the GNU General Public License and
+// a copy of the GCC Runtime Library Exception along with this program;
+// see the files COPYING3 and COPYING.RUNTIME respectively. If not, see
+// <http://www.gnu.org/licenses/>.
+
+/** @file ext/random.tcc
+ * This is an internal header file, included by other library headers.
+ * Do not attempt to use it directly. @headername{ext/random}
+ */
+
+#ifndef _EXT_RANDOM_TCC
+#define _EXT_RANDOM_TCC 1
+
+#pragma GCC system_header
+
+
+namespace __gnu_cxx _GLIBCXX_VISIBILITY(default)
+{
+_GLIBCXX_BEGIN_NAMESPACE_VERSION
+
+
+ template<typename _UIntType, size_t __m,
+ size_t __pos1, size_t __sl1, size_t __sl2,
+ size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2,
+ uint32_t __msk3, uint32_t __msk4,
+ uint32_t __parity1, uint32_t __parity2,
+ uint32_t __parity3, uint32_t __parity4>
+ void simd_fast_mersenne_twister_engine<_UIntType, __m,
+ __pos1, __sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4,
+ __parity1, __parity2, __parity3,
+ __parity4>::
+ seed(_UIntType __seed)
+ {
+ _M_state32[0] = static_cast<uint32_t>(__seed);
+ for (size_t __i = 1; __i < _M_nstate32; ++__i)
+ _M_state32[__i] = (1812433253UL
+ * (_M_state32[__i - 1] ^ (_M_state32[__i - 1] >> 30))
+ + __i);
+ _M_pos = state_size;
+ _M_period_certification();
+ }
+
+
+ namespace {
+
+ inline uint32_t _Func1(uint32_t __x)
+ {
+ return (__x ^ (__x >> 27)) * UINT32_C(1664525);
+ }
+
+ inline uint32_t _Func2(uint32_t __x)
+ {
+ return (__x ^ (__x >> 27)) * UINT32_C(1566083941);
+ }
+
+ }
+
+
+ template<typename _UIntType, size_t __m,
+ size_t __pos1, size_t __sl1, size_t __sl2,
+ size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2,
+ uint32_t __msk3, uint32_t __msk4,
+ uint32_t __parity1, uint32_t __parity2,
+ uint32_t __parity3, uint32_t __parity4>
+ template<typename _Sseq>
+ typename std::enable_if<std::is_class<_Sseq>::value>::type
+ simd_fast_mersenne_twister_engine<_UIntType, __m,
+ __pos1, __sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4,
+ __parity1, __parity2, __parity3,
+ __parity4>::
+ seed(_Sseq& __q)
+ {
+ size_t __lag;
+
+ if (_M_nstate32 >= 623)
+ __lag = 11;
+ else if (_M_nstate32 >= 68)
+ __lag = 7;
+ else if (_M_nstate32 >= 39)
+ __lag = 5;
+ else
+ __lag = 3;
+ const size_t __mid = (_M_nstate32 - __lag) / 2;
+
+ std::fill(_M_state32, _M_state32 + _M_nstate32, UINT32_C(0x8b8b8b8b));
+ uint32_t __arr[_M_nstate32];
+ __q.generate(__arr + 0, __arr + _M_nstate32);
+
+ uint32_t __r = _Func1(_M_state32[0] ^ _M_state32[__mid]
+ ^ _M_state32[_M_nstate32 - 1]);
+ _M_state32[__mid] += __r;
+ __r += _M_nstate32;
+ _M_state32[__mid + __lag] += __r;
+ _M_state32[0] = __r;
+
+ for (size_t __i = 1, __j = 0; __j < _M_nstate32; ++__j)
+ {
+ __r = _Func1(_M_state32[__i]
+ ^ _M_state32[(__i + __mid) % _M_nstate32]
+ ^ _M_state32[(__i + _M_nstate32 - 1) % _M_nstate32]);
+ _M_state32[(__i + __mid) % _M_nstate32] += __r;
+ __r += __arr[__j] + __i;
+ _M_state32[(__i + __mid + __lag) % _M_nstate32] += __r;
+ _M_state32[__i] = __r;
+ __i = (__i + 1) % _M_nstate32;
+ }
+ for (size_t __j = 0; __j < _M_nstate32; ++__j)
+ {
+ const size_t __i = (__j + 1) % _M_nstate32;
+ __r = _Func2(_M_state32[__i]
+ + _M_state32[(__i + __mid) % _M_nstate32]
+ + _M_state32[(__i + _M_nstate32 - 1) % _M_nstate32]);
+ _M_state32[(__i + __mid) % _M_nstate32] ^= __r;
+ __r -= __i;
+ _M_state32[(__i + __mid + __lag) % _M_nstate32] ^= __r;
+ _M_state32[__i] = __r;
+ }
+
+ _M_pos = state_size;
+ _M_period_certification();
+ }
+
+
+ template<typename _UIntType, size_t __m,
+ size_t __pos1, size_t __sl1, size_t __sl2,
+ size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2,
+ uint32_t __msk3, uint32_t __msk4,
+ uint32_t __parity1, uint32_t __parity2,
+ uint32_t __parity3, uint32_t __parity4>
+ void simd_fast_mersenne_twister_engine<_UIntType, __m,
+ __pos1, __sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4,
+ __parity1, __parity2, __parity3,
+ __parity4>::
+ _M_period_certification(void)
+ {
+ static const uint32_t __parity[4] = { __parity1, __parity2,
+ __parity3, __parity4 };
+ uint32_t __inner = 0;
+ for (size_t __i = 0; __i < 4; ++__i)
+ if (__parity[__i] != 0)
+ __inner ^= _M_state32[__i] & __parity[__i];
+
+ if (__builtin_parity(__inner) & 1)
+ return;
+ for (size_t __i = 0; __i < 4; ++__i)
+ if (__parity[__i] != 0)
+ {
+ _M_state32[__i] ^= 1 << (__builtin_ffs(__parity[__i]) - 1);
+ return;
+ }
+ __builtin_unreachable();
+ }
+
+
+ template<typename _UIntType, size_t __m,
+ size_t __pos1, size_t __sl1, size_t __sl2,
+ size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2,
+ uint32_t __msk3, uint32_t __msk4,
+ uint32_t __parity1, uint32_t __parity2,
+ uint32_t __parity3, uint32_t __parity4>
+ void simd_fast_mersenne_twister_engine<_UIntType, __m,
+ __pos1, __sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4,
+ __parity1, __parity2, __parity3,
+ __parity4>::
+ discard(unsigned long long __z)
+ {
+ while (__z > state_size - _M_pos)
+ {
+ __z -= state_size - _M_pos;
+
+ _M_gen_rand();
+ }
+
+ _M_pos += __z;
+ }
+
+
+#ifdef __SSE2__
+
+ namespace {
+
+ template<size_t __sl1, size_t __sl2, size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2, uint32_t __msk3, uint32_t __msk4>
+ inline __m128i __sse2_recursion(__m128i __a, __m128i __b,
+ __m128i __c, __m128i __d)
+ {
+ __m128i __y = _mm_srli_epi32(__b, __sr1);
+ __m128i __z = _mm_srli_si128(__c, __sr2);
+ __m128i __v = _mm_slli_epi32(__d, __sl1);
+ __z = _mm_xor_si128(__z, __a);
+ __z = _mm_xor_si128(__z, __v);
+ __m128i __x = _mm_slli_si128(__a, __sl2);
+ __y = _mm_and_si128(__y, _mm_set_epi32(__msk4, __msk3, __msk2, __msk1));
+ __z = _mm_xor_si128(__z, __x);
+ return _mm_xor_si128(__z, __y);
+ }
+
+ }
+
+
+ template<typename _UIntType, size_t __m,
+ size_t __pos1, size_t __sl1, size_t __sl2,
+ size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2,
+ uint32_t __msk3, uint32_t __msk4,
+ uint32_t __parity1, uint32_t __parity2,
+ uint32_t __parity3, uint32_t __parity4>
+ void simd_fast_mersenne_twister_engine<_UIntType, __m,
+ __pos1, __sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4,
+ __parity1, __parity2, __parity3,
+ __parity4>::
+ _M_gen_rand(void)
+ {
+ __m128i __r1 = _mm_load_si128(&_M_state[_M_nstate - 2]);
+ __m128i __r2 = _mm_load_si128(&_M_state[_M_nstate - 1]);
+
+ size_t __i;
+ for (__i = 0; __i < _M_nstate - __pos1; ++__i)
+ {
+ __m128i __r = __sse2_recursion<__sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4>
+ (_M_state[__i], _M_state[__i + __pos1], __r1, __r2);
+ _mm_store_si128(&_M_state[__i], __r);
+ __r1 = __r2;
+ __r2 = __r;
+ }
+ for (; __i < _M_nstate; ++__i)
+ {
+ __m128i __r = __sse2_recursion<__sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4>
+ (_M_state[__i], _M_state[__i + __pos1 - _M_nstate], __r1, __r2);
+ _mm_store_si128(&_M_state[__i], __r);
+ __r1 = __r2;
+ __r2 = __r;
+ }
+
+ _M_pos = 0;
+ }
+
+
+#else
+
+ namespace {
+
+ template<size_t __shift>
+ inline void __rshift(uint32_t *__out, const uint32_t *__in)
+ {
+ uint64_t __th = ((static_cast<uint64_t>(__in[3]) << 32)
+ | static_cast<uint64_t>(__in[2]));
+ uint64_t __tl = ((static_cast<uint64_t>(__in[1]) << 32)
+ | static_cast<uint64_t>(__in[0]));
+
+ uint64_t __oh = __th >> (__shift * 8);
+ uint64_t __ol = __tl >> (__shift * 8);
+ __ol |= __th << (64 - __shift * 8);
+ __out[1] = static_cast<uint32_t>(__ol >> 32);
+ __out[0] = static_cast<uint32_t>(__ol);
+ __out[3] = static_cast<uint32_t>(__oh >> 32);
+ __out[2] = static_cast<uint32_t>(__oh);
+ }
+
+
+ template<size_t __shift>
+ inline void __lshift(uint32_t *__out, const uint32_t *__in)
+ {
+ uint64_t __th = ((static_cast<uint64_t>(__in[3]) << 32)
+ | static_cast<uint64_t>(__in[2]));
+ uint64_t __tl = ((static_cast<uint64_t>(__in[1]) << 32)
+ | static_cast<uint64_t>(__in[0]));
+
+ uint64_t __oh = __th << (__shift * 8);
+ uint64_t __ol = __tl << (__shift * 8);
+ __oh |= __tl >> (64 - __shift * 8);
+ __out[1] = static_cast<uint32_t>(__ol >> 32);
+ __out[0] = static_cast<uint32_t>(__ol);
+ __out[3] = static_cast<uint32_t>(__oh >> 32);
+ __out[2] = static_cast<uint32_t>(__oh);
+ }
+
+
+ template<size_t __sl1, size_t __sl2, size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2, uint32_t __msk3, uint32_t __msk4>
+ inline void __recursion(uint32_t *__r,
+ const uint32_t *__a, const uint32_t *__b,
+ const uint32_t *__c, const uint32_t *__d)
+ {
+ uint32_t __x[4];
+ uint32_t __y[4];
+
+ __lshift<__sl2>(__x, __a);
+ __rshift<__sr2>(__y, __c);
+ __r[0] = (__a[0] ^ __x[0] ^ ((__b[0] >> __sr1) & __msk1)
+ ^ __y[0] ^ (__d[0] << __sl1));
+ __r[1] = (__a[1] ^ __x[1] ^ ((__b[1] >> __sr1) & __msk2)
+ ^ __y[1] ^ (__d[1] << __sl1));
+ __r[2] = (__a[2] ^ __x[2] ^ ((__b[2] >> __sr1) & __msk3)
+ ^ __y[2] ^ (__d[2] << __sl1));
+ __r[3] = (__a[3] ^ __x[3] ^ ((__b[3] >> __sr1) & __msk4)
+ ^ __y[3] ^ (__d[3] << __sl1));
+ }
+
+ }
+
+
+ template<typename _UIntType, size_t __m,
+ size_t __pos1, size_t __sl1, size_t __sl2,
+ size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2,
+ uint32_t __msk3, uint32_t __msk4,
+ uint32_t __parity1, uint32_t __parity2,
+ uint32_t __parity3, uint32_t __parity4>
+ void simd_fast_mersenne_twister_engine<_UIntType, __m,
+ __pos1, __sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4,
+ __parity1, __parity2, __parity3,
+ __parity4>::
+ _M_gen_rand(void)
+ {
+ const uint32_t *__r1 = &_M_state32[_M_nstate32 - 8];
+ const uint32_t *__r2 = &_M_state32[_M_nstate32 - 4];
+ static constexpr size_t __pos1_32 = __pos1 * 4;
+
+ size_t __i;
+ for (__i = 0; __i < _M_nstate32 - __pos1_32; __i += 4)
+ {
+ __recursion<__sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4>
+ (&_M_state32[__i], &_M_state32[__i],
+ &_M_state32[__i + __pos1_32], __r1, __r2);
+ __r1 = __r2;
+ __r2 = &_M_state32[__i];
+ }
+
+ for (; __i < _M_nstate32; __i += 4)
+ {
+ __recursion<__sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4>
+ (&_M_state32[__i], &_M_state32[__i],
+ &_M_state32[__i + __pos1_32 - _M_nstate32], __r1, __r2);
+ __r1 = __r2;
+ __r2 = &_M_state32[__i];
+ }
+
+ _M_pos = 0;
+ }
+
+#endif
+
+
+ template<typename _UIntType, size_t __m,
+ size_t __pos1, size_t __sl1, size_t __sl2,
+ size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2,
+ uint32_t __msk3, uint32_t __msk4,
+ uint32_t __parity1, uint32_t __parity2,
+ uint32_t __parity3, uint32_t __parity4,
+ typename _CharT, typename _Traits>
+ std::basic_ostream<_CharT, _Traits>&
+ operator<<(std::basic_ostream<_CharT, _Traits>& __os,
+ const __gnu_cxx::simd_fast_mersenne_twister_engine<_UIntType,
+ __m, __pos1, __sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4,
+ __parity1, __parity2, __parity3, __parity4>& __x)
+ {
+ typedef std::basic_ostream<_CharT, _Traits> __ostream_type;
+ typedef typename __ostream_type::ios_base __ios_base;
+
+ const typename __ios_base::fmtflags __flags = __os.flags();
+ const _CharT __fill = __os.fill();
+ const _CharT __space = __os.widen(' ');
+ __os.flags(__ios_base::dec | __ios_base::fixed | __ios_base::left);
+ __os.fill(__space);
+
+ for (size_t __i = 0; __i < __x._M_nstate32; ++__i)
+ __os << __x._M_state32[__i] << __space;
+ __os << __x._M_pos;
+
+ __os.flags(__flags);
+ __os.fill(__fill);
+ return __os;
+ }
+
+
+ template<typename _UIntType, size_t __m,
+ size_t __pos1, size_t __sl1, size_t __sl2,
+ size_t __sr1, size_t __sr2,
+ uint32_t __msk1, uint32_t __msk2,
+ uint32_t __msk3, uint32_t __msk4,
+ uint32_t __parity1, uint32_t __parity2,
+ uint32_t __parity3, uint32_t __parity4,
+ typename _CharT, typename _Traits>
+ std::basic_istream<_CharT, _Traits>&
+ operator>>(std::basic_istream<_CharT, _Traits>& __is,
+ __gnu_cxx::simd_fast_mersenne_twister_engine<_UIntType,
+ __m, __pos1, __sl1, __sl2, __sr1, __sr2,
+ __msk1, __msk2, __msk3, __msk4,
+ __parity1, __parity2, __parity3, __parity4>& __x)
+ {
+ typedef std::basic_istream<_CharT, _Traits> __istream_type;
+ typedef typename __istream_type::ios_base __ios_base;
+
+ const typename __ios_base::fmtflags __flags = __is.flags();
+ __is.flags(__ios_base::dec | __ios_base::skipws);
+
+ for (size_t __i = 0; __i < __x._M_nstate32; ++__i)
+ __is >> __x._M_state32[__i];
+ __is >> __x._M_pos;
+
+ __is.flags(__flags);
+ return __is;
+ }
+
+_GLIBCXX_END_NAMESPACE_VERSION
+} // namespace
+
+
+#endif // _EXT_RANDOM_TCC
--- /dev/null
+// { dg-options "-std=c++0x" }
+// { dg-require-cstdint "" }
+//
+// 2010-02-16 Paolo Carlini <paolo.carlini@oracle.com>
+// 2012-08-28 Ulrich Drepper <drepper@gmail.com>, adapted for SFMT
+//
+// Copyright (C) 2010, 2012 Free Software Foundation, Inc.
+//
+// This file is part of the GNU ISO C++ Library. This library is free
+// software; you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the
+// Free Software Foundation; either version 3, or (at your option)
+// any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License along
+// with this library; see the file COPYING3. If not see
+// <http://www.gnu.org/licenses/>.
+
+#include <ext/random>
+
+void
+test01()
+{
+ typedef uint32_t value_type;
+
+ __gnu_cxx::simd_fast_mersenne_twister_engine<value_type, 607, 2,
+ 15, 3, 13, 3,
+ 0xfdff37ffU, 0xef7f3f7dU,
+ 0xff777b7dU, 0x7ff7fb2fU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x5986f054U> e(1);
+
+ const auto f(e);
+ auto g(f);
+ g = g; // Suppress unused warning.
+}
+
+int main()
+{
+ test01();
+ return 0;
+}
--- /dev/null
+// { dg-options "-std=c++0x" }
+// { dg-require-cstdint "" }
+//
+// 2008-11-24 Edward M. Smith-Rowland <3dw4rd@verizon.net>
+// 2012-08-28 Ulrich Drepper <drepper@gmail.com>, adapted for SFMT
+//
+// Copyright (C) 2008, 2009, 2010, 2012 Free Software Foundation, Inc.
+//
+// This file is part of the GNU ISO C++ Library. This library is free
+// software; you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the
+// Free Software Foundation; either version 3, or (at your option)
+// any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License along
+// with this library; see the file COPYING3. If not see
+// <http://www.gnu.org/licenses/>.
+
+#include <iostream>
+#include <ext/random>
+#include <testsuite_hooks.h>
+
+
+template<class SFMT>
+void run_test(std::initializer_list<typename SFMT::result_type> vals)
+{
+ typedef typename SFMT::result_type result_type;
+ SFMT e;
+
+ e.seed(sizeof(result_type) == 4 ? 1234 : 4321);
+ e.discard(990);
+ bool success = true;
+ for (auto i : vals)
+ {
+ result_type r = e();
+ success &= r == i;
+ std::cout << r << " vs " << i << std::endl;
+ }
+ VERIFY( success );
+}
+
+
+void
+test01()
+{
+ bool test __attribute__((unused)) = true;
+
+ __gnu_cxx::sfmt19937 e;
+
+ VERIFY( e.min() == 0 );
+ VERIFY( e.max() == std::numeric_limits<uint32_t>::max() );
+
+ run_test<__gnu_cxx::sfmt607>({ UINT32_C(1318548553), UINT32_C(1985957974),
+ UINT32_C(1367744196), UINT32_C(3463392791), UINT32_C(2780736231),
+ UINT32_C(3894488561), UINT32_C(3157036262), UINT32_C(3491812767),
+ UINT32_C(1724574180), UINT32_C(3645035493) });
+
+ run_test<__gnu_cxx::sfmt607_64>({UINT64_C(15510024334182072935),
+ UINT64_C(5793753331747412752), UINT64_C(16198353238554625740),
+ UINT64_C(2233208824926016498), UINT64_C(3566091399820823780),
+ UINT64_C(16608268514591292798), UINT64_C(10684941689666043359),
+ UINT64_C(12463424292910456802), UINT64_C(5902567440240131366),
+ UINT64_C(7228030834036501150) });
+
+ run_test<__gnu_cxx::sfmt1279>({ UINT32_C(66657331), UINT32_C(637106837),
+ UINT32_C(406927341), UINT32_C(3964420203), UINT32_C(2127134160),
+ UINT32_C(1327235047), UINT32_C(227339400), UINT32_C(97109542),
+ UINT32_C(1814799261), UINT32_C(340888197) });
+
+ run_test<__gnu_cxx::sfmt1279_64>({ UINT64_C(16431921382083697129),
+ UINT64_C(3107599092104940900), UINT64_C(4055245506102959965),
+ UINT64_C(16096064917153424198), UINT64_C(14429331498726837109),
+ UINT64_C(9539664361920633782), UINT64_C(1435296568185387099),
+ UINT64_C(15922567183295047131), UINT64_C(641988285517426228),
+ UINT64_C(15936274870984512675) });
+
+ run_test<__gnu_cxx::sfmt2281>({ UINT32_C(2662391944), UINT32_C(1176696104),
+ UINT32_C(3587947451), UINT32_C(4098993357), UINT32_C(3140998698),
+ UINT32_C(870759742), UINT32_C(623529127), UINT32_C(3458807285),
+ UINT32_C(3341615957), UINT32_C(195614711) });
+
+ run_test<__gnu_cxx::sfmt2281_64>({ UINT64_C(16747191622237074632),
+ UINT64_C(15804170396401370381), UINT64_C(3395175278324920203),
+ UINT64_C(1541877340159274442), UINT64_C(14176322102994316687),
+ UINT64_C(5130618305074712143), UINT64_C(6769693652413407081),
+ UINT64_C(17733765687477661079), UINT64_C(5189766940360047353),
+ UINT64_C(1333654688569723389) });
+
+ run_test<__gnu_cxx::sfmt4253>({ UINT32_C(90342133), UINT32_C(1083987943),
+ UINT32_C(1785481425), UINT32_C(1921212667), UINT32_C(3164342992),
+ UINT32_C(1489324569), UINT32_C(603530523), UINT32_C(952851722),
+ UINT32_C(2380944844), UINT32_C(3335854133) });
+
+ run_test<__gnu_cxx::sfmt4253_64>({ UINT64_C(11570915401962514263),
+ UINT64_C(206693220452528225), UINT64_C(16553299974633247759),
+ UINT64_C(1069562842508952901), UINT64_C(7203975672387749585),
+ UINT64_C(7552781925224963166), UINT64_C(16865729458807008705),
+ UINT64_C(7848963629493506078), UINT64_C(9282397173969292817),
+ UINT64_C(10738488504584559289) });
+
+ run_test<__gnu_cxx::sfmt11213>({ UINT32_C(2072997009), UINT32_C(1332330347),
+ UINT32_C(179681555), UINT32_C(2315290438), UINT32_C(2429393974),
+ UINT32_C(509881964), UINT32_C(3807607878), UINT32_C(3055319970),
+ UINT32_C(671840881), UINT32_C(3477325874) });
+
+ run_test<__gnu_cxx::sfmt11213_64>({ UINT64_C(373867573626408653),
+ UINT64_C(4732829340233638861), UINT64_C(16174630176505735656),
+ UINT64_C(10063018133994900869), UINT64_C(17308645173308419196),
+ UINT64_C(11091353816581371951), UINT64_C(15078420471318089727),
+ UINT64_C(17965717592743818706), UINT64_C(12301543162252389155),
+ UINT64_C(1724943167823308511) });
+
+ run_test<__gnu_cxx::sfmt19937>({ UINT32_C(4002809368), UINT32_C(421169044),
+ UINT32_C(1112642589), UINT32_C(3076213779), UINT32_C(3387033971),
+ UINT32_C(2499610950), UINT32_C(3057240914), UINT32_C(1662679783),
+ UINT32_C(461224431), UINT32_C(1168395933) });
+
+ run_test<__gnu_cxx::sfmt19937_64>({ UINT64_C(8032857516355555296),
+ UINT64_C(14023605983059313116), UINT64_C(1032336061815461376),
+ UINT64_C(9840995337876562612), UINT64_C(9869256223029203587),
+ UINT64_C(12227975697177267636), UINT64_C(12728115115844186033),
+ UINT64_C(7752058479783205470), UINT64_C(729733219713393087),
+ UINT64_C(12954017801239007622) });
+
+ run_test<__gnu_cxx::sfmt44497>({ UINT32_C(1483092082), UINT32_C(1895679637),
+ UINT32_C(9122740), UINT32_C(635864575), UINT32_C(320732971),
+ UINT32_C(4253159584), UINT32_C(30097521), UINT32_C(839233316),
+ UINT32_C(1431693534), UINT32_C(645981752) });
+
+ run_test<__gnu_cxx::sfmt44497_64>({ UINT64_C(6246103978016445638),
+ UINT64_C(4198275826138953222), UINT64_C(12473679170573289212),
+ UINT64_C(14745709982748360209), UINT64_C(3630790792408208113),
+ UINT64_C(4195294399578350499), UINT64_C(3742595698794327253),
+ UINT64_C(17388385867517445933), UINT64_C(4261866397667814669),
+ UINT64_C(17394085161690598095) });
+
+ run_test<__gnu_cxx::sfmt86243>({ UINT32_C(3910985535), UINT32_C(100094501),
+ UINT32_C(3120362616), UINT32_C(1854432382), UINT32_C(314688154),
+ UINT32_C(522122712), UINT32_C(3026095676), UINT32_C(3681962735),
+ UINT32_C(1851548627), UINT32_C(2153846465) });
+
+ run_test<__gnu_cxx::sfmt86243_64>({ UINT64_C(250135615696586029),
+ UINT64_C(4836277875486422184), UINT64_C(12389320296183057446),
+ UINT64_C(7983028875884559442), UINT64_C(10079555227308335361),
+ UINT64_C(14829333386540244214), UINT64_C(12159744972103351172),
+ UINT64_C(4932579842314286356), UINT64_C(5200375244476537050),
+ UINT64_C(11795681221121010641) });
+
+ run_test<__gnu_cxx::sfmt132049>({ UINT32_C(1551023420), UINT32_C(1462317554),
+ UINT32_C(2882528449), UINT32_C(1136299843), UINT32_C(292840589),
+ UINT32_C(1307775247), UINT32_C(463274356), UINT32_C(1430357686),
+ UINT32_C(3907607055), UINT32_C(3462509184) });
+
+ run_test<__gnu_cxx::sfmt132049_64>({ UINT64_C(649482638765113922),
+ UINT64_C(14205859353699897918), UINT64_C(14077261854908137257),
+ UINT64_C(9564785861212212042), UINT64_C(7310747921257808846),
+ UINT64_C(13759009477111470372), UINT64_C(11942123860149328831),
+ UINT64_C(12868386070200572127), UINT64_C(18348617059674004332),
+ UINT64_C(4233208019331956061) });
+
+ run_test<__gnu_cxx::sfmt216091>({ UINT32_C(4171954654), UINT32_C(2938491210),
+ UINT32_C(1356393891), UINT32_C(3558249995), UINT32_C(3711769979),
+ UINT32_C(3434953144), UINT32_C(1601628304), UINT32_C(2187495640),
+ UINT32_C(1762169715), UINT32_C(2141213778) });
+
+ run_test<__gnu_cxx::sfmt216091_64>({ UINT64_C(11322404276387828766),
+ UINT64_C(9653391575000195546), UINT64_C(1767839622905368464),
+ UINT64_C(1690838241348740821), UINT64_C(817628268513271254),
+ UINT64_C(15111277786569319196), UINT64_C(15817118810543358764),
+ UINT64_C(5639001693408668083), UINT64_C(9959854003669400568),
+ UINT64_C(13675983279642398887) });
+}
+
+int main()
+{
+ test01();
+ return 0;
+}
--- /dev/null
+// { dg-options "-std=c++0x" }
+// { dg-require-cstdint "" }
+//
+// 2008-11-24 Edward M. Smith-Rowland <3dw4rd@verizon.net>
+// 2012-08-28 Ulrich Drepper <drepper@gmail.com>, adapted for SFMT
+//
+// Copyright (C) 2008, 2009, 2010, 2012 Free Software Foundation, Inc.
+//
+// This file is part of the GNU ISO C++ Library. This library is free
+// software; you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the
+// Free Software Foundation; either version 3, or (at your option)
+// any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License along
+// with this library; see the file COPYING3. If not see
+// <http://www.gnu.org/licenses/>.
+
+#include <ext/random>
+
+void
+test01()
+{
+ unsigned long seed = 2;
+ __gnu_cxx::simd_fast_mersenne_twister_engine<
+ uint64_t, 607, 2,
+ 15, 3, 13, 3,
+ 0xfdff37ffU, 0xef7f3f7dU,
+ 0xff777b7dU, 0x7ff7fb2fU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x5986f054U> x(seed);
+}
+
+int main()
+{
+ test01();
+ return 0;
+}
--- /dev/null
+// { dg-options "-std=c++0x" }
+// { dg-require-cstdint "" }
+//
+// 2008-11-24 Edward M. Smith-Rowland <3dw4rd@verizon.net>
+// 2012-08-28 Ulrich Drepper <drepper@gmail.com>, adapted for SFMT
+//
+// Copyright (C) 2008, 2009, 2010, 2012 Free Software Foundation, Inc.
+//
+// This file is part of the GNU ISO C++ Library. This library is free
+// software; you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the
+// Free Software Foundation; either version 3, or (at your option)
+// any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License along
+// with this library; see the file COPYING3. If not see
+// <http://www.gnu.org/licenses/>.
+
+#include <ext/random>
+
+void
+test01()
+{
+ double seed = 2.0;
+ __gnu_cxx::simd_fast_mersenne_twister_engine<
+ uint64_t, 607, 2,
+ 15, 3, 13, 3,
+ 0xfdff37ffU, 0xef7f3f7dU,
+ 0xff777b7dU, 0x7ff7fb2fU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x5986f054U> x(seed);
+}
+
+int main()
+{
+ test01();
+ return 0;
+}
--- /dev/null
+// { dg-options "-std=c++0x" }
+// { dg-require-cstdint "" }
+//
+// 2010-02-01 Paolo Carlini <paolo.carlini@oracle.com>
+// 2012-08-28 Ulrich Drepper <drepper@gmail.com>, adapted for SFMT
+//
+// Copyright (C) 2010, 2012 Free Software Foundation, Inc.
+//
+// This file is part of the GNU ISO C++ Library. This library is free
+// software; you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the
+// Free Software Foundation; either version 3, or (at your option)
+// any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License along
+// with this library; see the file COPYING3. If not see
+// <http://www.gnu.org/licenses/>.
+
+#include <ext/random>
+
+void
+test01()
+{
+ std::seed_seq seed;
+ __gnu_cxx::simd_fast_mersenne_twister_engine<
+ uint32_t, 607, 2,
+ 15, 3, 13, 3,
+ 0xfdff37ffU, 0xef7f3f7dU,
+ 0xff777b7dU, 0x7ff7fb2fU,
+ 0x00000001U, 0x00000000U,
+ 0x00000000U, 0x5986f054U> x(seed);
+}
+
+int main()
+{
+ test01();
+ return 0;
+}
--- /dev/null
+// { dg-options "-std=c++0x" }
+// { dg-require-cstdint "" }
+//
+// 2008-11-24 Edward M. Smith-Rowland <3dw4rd@verizon.net>
+// 2012-08-28 Ulrich Drepper <drepper@gmail.com>, adapted for SFMT
+//
+// Copyright (C) 2008, 2009, 2012 Free Software Foundation, Inc.
+//
+// This file is part of the GNU ISO C++ Library. This library is free
+// software; you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the
+// Free Software Foundation; either version 3, or (at your option)
+// any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License along
+// with this library; see the file COPYING3. If not see
+// <http://www.gnu.org/licenses/>.
+
+// 26.4.3.2 Class template mersenne_twister_engine [rand.eng.mers]
+// 26.4.2.2 Concept RandomNumberEngine [rand.concept.eng]
+
+#include <ext/random>
+#include <testsuite_hooks.h>
+
+void
+test01()
+{
+ bool test __attribute__((unused)) = true;
+
+ __gnu_cxx::sfmt19937 u, v;
+
+ VERIFY( u == v );
+
+ u.discard(100);
+ v.discard(100);
+
+ VERIFY( u == v );
+}
+
+int main()
+{
+ test01();
+ return 0;
+}
--- /dev/null
+// { dg-options "-std=c++0x" }
+// { dg-require-cstdint "" }
+//
+// 2010-03-16 Paolo Carlini <paolo.carlini@oracle.com>
+// 2012-08-28 Ulrich Drepper <drepper@gmail.com>, adapted for SFMT
+//
+// Copyright (C) 2010, 2012 Free Software Foundation, Inc.
+//
+// This file is part of the GNU ISO C++ Library. This library is free
+// software; you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the
+// Free Software Foundation; either version 3, or (at your option)
+// any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License along
+// with this library; see the file COPYING3. If not see
+// <http://www.gnu.org/licenses/>.
+
+// 26.5.3.2 Class template mersenne_twister_engine [rand.eng.mers]
+
+#include <ext/random>
+#include <testsuite_hooks.h>
+
+void
+test01()
+{
+ bool test __attribute__((unused)) = true;
+
+ __gnu_cxx::sfmt19937 u, v;
+
+ VERIFY( !(u != v) );
+
+ u.discard(100);
+ v.discard(100);
+
+ VERIFY( !(u != v) );
+
+ v.discard(1);
+
+ VERIFY( u != v );
+}
+
+int main()
+{
+ test01();
+ return 0;
+}
--- /dev/null
+// { dg-options "-std=c++0x" }
+// { dg-require-cstdint "" }
+//
+// 2008-11-24 Edward M. Smith-Rowland <3dw4rd@verizon.net>
+// 2012-08-28 Ulrich Drepper <drepper@gmail.com>, adapted for SFMT
+//
+// Copyright (C) 2008, 2009, 2010, 2011, 2012 Free Software Foundation, Inc.
+//
+// This file is part of the GNU ISO C++ Library. This library is free
+// software; you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the
+// Free Software Foundation; either version 3, or (at your option)
+// any later version.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License along
+// with this library; see the file COPYING3. If not see
+// <http://www.gnu.org/licenses/>.
+
+// 26.4.3.2 Class template mersenne_twister_engine [rand.eng.mers]
+// 26.4.2.2 Concept RandomNumberEngine [rand.concept.eng]
+
+#include <sstream>
+#include <ext/random>
+#include <testsuite_hooks.h>
+
+void
+test01()
+{
+ bool test __attribute__((unused)) = true;
+
+ std::stringstream str;
+ __gnu_cxx::sfmt19937 u, v;
+
+ u(); // advance
+ str << u;
+
+ VERIFY( !(u == v) );
+
+ str >> v;
+ VERIFY( u == v );
+ for (unsigned i = 0; i < 1000; ++i)
+ VERIFY( u() == v() );
+
+ str.clear();
+ str << v;
+
+ u();
+ u();
+ u();
+
+ str >> u;
+ VERIFY( u == v );
+ for (unsigned i = 0; i < 1000; ++i)
+ VERIFY( u() == v() );
+}
+
+int main()
+{
+ test01();
+ return 0;
+}
projects / thirdparty / gcc.git / commitdiff
