test_mt19937.cpp 2.1 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081
  1. /* test_mt19937.cpp
  2. *
  3. * Copyright Steven Watanabe 2011
  4. * Distributed under the Boost Software License, Version 1.0. (See
  5. * accompanying file LICENSE_1_0.txt or copy at
  6. * http://www.boost.org/LICENSE_1_0.txt)
  7. *
  8. * $Id$
  9. *
  10. */
  11. #include <boost/random/mersenne_twister.hpp>
  12. #include <algorithm>
  13. #include <vector>
  14. #include <boost/cstdint.hpp>
  15. #define BOOST_RANDOM_URNG boost::random::mt19937
  16. #define BOOST_RANDOM_SEED_WORDS 624
  17. // validation by experiment from mt19937.c
  18. #define BOOST_RANDOM_VALIDATION_VALUE 4123659995U
  19. #define BOOST_RANDOM_SEED_SEQ_VALIDATION_VALUE 666528879U
  20. #define BOOST_RANDOM_ITERATOR_VALIDATION_VALUE 3408548740U
  21. #define BOOST_RANDOM_DISCARD_COUNT1 9307
  22. #define BOOST_RANDOM_DISCARD_COUNT2 20000000
  23. #define BOOST_RANDOM_DISCARD_MAX
  24. #define BOOST_RANDOM_GENERATE_VALUES { 0xD091BB5CU, 0x22AE9EF6U, 0xE7E1FAEEU, 0xD5C31F79U }
  25. #include "test_generator.ipp"
  26. struct seed_seq_0 {
  27. template<class It>
  28. void generate(It begin, It end) const {
  29. std::fill(begin, end, boost::uint32_t(0));
  30. }
  31. };
  32. struct seed_seq_1 {
  33. template<class It>
  34. void generate(It begin, It end) const {
  35. std::fill(begin, end, boost::uint32_t(0));
  36. *(end - 1) = 1;
  37. }
  38. };
  39. BOOST_AUTO_TEST_CASE(test_special_seed) {
  40. {
  41. seed_seq_1 seed;
  42. std::vector<boost::uint32_t> vec(624);
  43. seed.generate(vec.begin(), vec.end());
  44. std::vector<boost::uint32_t>::iterator it = vec.begin();
  45. boost::mt19937 gen1(it, vec.end());
  46. BOOST_CHECK_EQUAL(gen1(), 0u);
  47. BOOST_CHECK_EQUAL(gen1(), 0u);
  48. boost::mt19937 gen2(seed);
  49. BOOST_CHECK_EQUAL(gen2(), 0u);
  50. BOOST_CHECK_EQUAL(gen2(), 0u);
  51. BOOST_CHECK_EQUAL(gen1, gen2);
  52. }
  53. {
  54. seed_seq_0 seed;
  55. std::vector<boost::uint32_t> vec(624);
  56. seed.generate(vec.begin(), vec.end());
  57. std::vector<boost::uint32_t>::iterator it = vec.begin();
  58. boost::mt19937 gen1(it, vec.end());
  59. BOOST_CHECK_EQUAL(gen1(), 1141379330u);
  60. BOOST_CHECK_EQUAL(gen1(), 0u);
  61. boost::mt19937 gen2(seed);
  62. BOOST_CHECK_EQUAL(gen2(), 1141379330u);
  63. BOOST_CHECK_EQUAL(gen2(), 0u);
  64. BOOST_CHECK_EQUAL(gen1, gen2);
  65. }
  66. }