#ifndef SPROUT_RANDOM_UNIFORM_SMALLINT_HPP #define SPROUT_RANDOM_UNIFORM_SMALLINT_HPP #include #include #include #include #include #include #include namespace sprout { namespace random { // // uniform_smallint // template class uniform_smallint { public: typedef IntType input_type; typedef IntType result_type; private: static SPROUT_CONSTEXPR IntType arg_check(IntType min_arg, IntType max_arg) { return min_arg <= max_arg ? min_arg : throw "assert(min_arg <= max_arg)" ; } public: // // param_type // class param_type { public: typedef uniform_smallint distribution_type; private: IntType min_; IntType max_; public: SPROUT_CONSTEXPR param_type() : min_(0) , max_(9) {} SPROUT_CONSTEXPR explicit param_type(IntType min_arg, IntType max_arg = 9) : min_(arg_check(min_arg, max_arg)) , max_(max_arg) {} SPROUT_CONSTEXPR IntType a() const { return min_; } SPROUT_CONSTEXPR IntType b() const { return max_; } template friend std::basic_ostream& operator>>( std::basic_istream& lhs, param_type const& rhs ) { return lhs >> rhs.min_ >> std::ws >> rhs.max_; } template friend std::basic_ostream& operator<<( std::basic_ostream& lhs, param_type const& rhs ) { return lhs << rhs.min_ << " " << rhs.max_; } SPROUT_CONSTEXPR friend bool operator==(param_type const& lhs, param_type const& rhs) { return lhs.min_ == rhs.min_ && lhs.max_ == rhs.max_; } SPROUT_CONSTEXPR friend bool operator!=(param_type const& lhs, param_type const& rhs) { return !(lhs == rhs); } }; private: result_type min_; result_type max_; private: template SPROUT_CONSTEXPR sprout::random::random_result generate_true_2( Engine const& eng, sprout::random::random_result const& rnd, RangeType range, BaseUnsigned base_range, BaseUnsigned val ) const { return range >= base_range ? sprout::random::random_result( sprout::random::detail::add()(static_cast(val), min_), rnd.engine(), *this ) : sprout::random::random_result( sprout::random::detail::add()(static_cast(val % (static_cast(range) + 1)), min_), rnd.engine(), *this ) ; } template SPROUT_CONSTEXPR sprout::random::random_result generate_true_1( Engine const& eng, sprout::random::random_result const& rnd, RangeType range, BaseUnsigned base_range ) const { typedef typename Engine::result_type base_result; return generate_true_2( eng, rnd, range, base_range, BaseUnsigned(sprout::random::detail::subtract()(rnd.result(), eng.min())) ); } template SPROUT_CONSTEXPR sprout::random::random_result generate( Engine const& eng, std::true_type ) const { typedef typename Engine::result_type base_result; typedef typename std::make_unsigned::type base_unsigned; typedef typename std::make_unsigned::type range_type; return generate_true_1( eng, eng(), range_type(sprout::random::detail::subtract()(max_, min_)), base_unsigned(sprout::random::detail::subtract()(eng.max(), eng.min())) ); } template SPROUT_CONSTEXPR sprout::random::random_result generate_false_2( Engine const& eng, sprout::random::random_result > const& rnd, RangeType range, RangeType offset ) const { return offset > range ? sprout::random::random_result( max_, rnd.engine(), *this ) : sprout::random::random_result( sprout::random::detail::add()(offset , min_), rnd.engine(), *this ) ; } template SPROUT_CONSTEXPR sprout::random::random_result generate_false_1( Engine const& eng, sprout::random::random_result > const& rnd, RangeType range ) const { typedef typename Engine::result_type base_result; return generate_false_2( eng, rnd, RangeType(sprout::random::detail::subtract()(max_, min_)), RangeType(static_cast(rnd.result() * (static_cast(range) + 1))) ); } template SPROUT_CONSTEXPR sprout::random::random_result generate( Engine const& eng, std::false_type ) const { typedef typename Engine::result_type base_result; typedef typename std::make_unsigned::type range_type; return generate_false_1( eng, sprout::random::uniform_01()(eng), range_type(sprout::random::detail::subtract()(max_, min_)) ); } public: SPROUT_CONSTEXPR uniform_smallint() : min_(0) , max_(9) {} SPROUT_CONSTEXPR explicit uniform_smallint(IntType min_arg, IntType max_arg = 9) : min_(arg_check(min_arg, max_arg)) , max_(max_arg) {} explicit uniform_smallint(param_type const& parm) : min_(parm.a()) , max_(parm.b()) {} SPROUT_CONSTEXPR result_type a() const { return min_; } SPROUT_CONSTEXPR result_type b() const { return max_; } SPROUT_CONSTEXPR result_type min() const { return min_; } SPROUT_CONSTEXPR result_type max() const { return max_; } SPROUT_CONSTEXPR param_type param() const { return param_type(min_, max_); } void param(param_type const& parm) { min_ = parm.a(); max_ = parm.b(); } template SPROUT_CONSTEXPR sprout::random::random_result operator()(Engine const& eng) const { typedef typename Engine::result_type base_result; return generate(eng, typename std::is_integral::type()); } template friend std::basic_ostream& operator>>( std::basic_istream& lhs, uniform_smallint const& rhs ) { param_type parm; return lhs >> parm; param(parm); return lhs; } template friend std::basic_ostream& operator<<( std::basic_ostream& lhs, uniform_smallint const& rhs ) { return lhs << param(); } SPROUT_CONSTEXPR friend bool operator==(uniform_smallint const& lhs, uniform_smallint const& rhs) { return lhs.param() == rhs.param(); } SPROUT_CONSTEXPR friend bool operator!=(uniform_smallint const& lhs, uniform_smallint const& rhs) { return !(lhs == rhs); } }; } // namespace random using sprout::random::uniform_smallint; } // namespace sprout #endif // #ifndef SPROUT_RANDOM_UNIFORM_SMALLINT_HPP