/* * Copyright 2008-2021 NVIDIA Corporation * * Licensed under the Apache License, Version 2.0 (the "License"); * you may not use this file except in compliance with the License. * You may obtain a copy of the License at * * http://www.apache.org/licenses/LICENSE-2.0 * * Unless required by applicable law or agreed to in writing, software * distributed under the License is distributed on an "AS IS" BASIS, * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. * See the License for the specific language governing permissions and * limitations under the License. */ #pragma once #include #if defined(_CCCL_IMPLICIT_SYSTEM_HEADER_GCC) # pragma GCC system_header #elif defined(_CCCL_IMPLICIT_SYSTEM_HEADER_CLANG) # pragma clang system_header #elif defined(_CCCL_IMPLICIT_SYSTEM_HEADER_MSVC) # pragma system_header #endif // no system header #include #include #include THRUST_NAMESPACE_BEGIN namespace random { template _CCCL_HOST_DEVICE uniform_int_distribution::uniform_int_distribution(IntType a, IntType b) : m_param(a, b) {} // end uniform_int_distribution::uniform_int_distribution() template _CCCL_HOST_DEVICE uniform_int_distribution::uniform_int_distribution(const param_type& parm) : m_param(parm) {} // end uniform_int_distribution::uniform_int_distribution() template _CCCL_HOST_DEVICE void uniform_int_distribution::reset() {} // end uniform_int_distribution::reset() template template _CCCL_HOST_DEVICE typename uniform_int_distribution::result_type uniform_int_distribution::operator()(UniformRandomNumberGenerator& urng) { return operator()(urng, m_param); } // end uniform_int_distribution::operator()() template template _CCCL_HOST_DEVICE typename uniform_int_distribution::result_type uniform_int_distribution::operator()(UniformRandomNumberGenerator& urng, const param_type& parm) { // XXX this implementation is somewhat hacky and will skip // values if the range of the RNG is smaller than the range of the distribution // we should improve this implementation in a later version typedef typename thrust::detail::largest_available_float::type float_type; const float_type real_min(static_cast(parm.first)); const float_type real_max(static_cast(parm.second)); // add one to the right end of the interval because it is half-open // XXX adding 1.0 to a potentially large floating point number seems like a bad idea uniform_real_distribution real_dist(real_min, real_max + float_type(1)); return static_cast(real_dist(urng)); } // end uniform_int_distribution::operator()() template _CCCL_HOST_DEVICE typename uniform_int_distribution::result_type uniform_int_distribution::a() const { return m_param.first; } // end uniform_int_distribution::a() template _CCCL_HOST_DEVICE typename uniform_int_distribution::result_type uniform_int_distribution::b() const { return m_param.second; } // end uniform_int_distribution::b() template _CCCL_HOST_DEVICE typename uniform_int_distribution::param_type uniform_int_distribution::param() const { return m_param; } // end uniform_int_distribution::param() template _CCCL_HOST_DEVICE void uniform_int_distribution::param(const param_type& parm) { m_param = parm; } // end uniform_int_distribution::param() template _CCCL_HOST_DEVICE typename uniform_int_distribution::result_type uniform_int_distribution::min THRUST_PREVENT_MACRO_SUBSTITUTION() const { return a(); } // end uniform_int_distribution::min() template _CCCL_HOST_DEVICE typename uniform_int_distribution::result_type uniform_int_distribution::max THRUST_PREVENT_MACRO_SUBSTITUTION() const { return b(); } // end uniform_int_distribution::max() template _CCCL_HOST_DEVICE bool uniform_int_distribution::equal(const uniform_int_distribution& rhs) const { return param() == rhs.param(); } template template std::basic_ostream& uniform_int_distribution::stream_out(std::basic_ostream& os) const { typedef std::basic_ostream ostream_type; typedef typename ostream_type::ios_base ios_base; // save old flags and fill character 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); os << a() << space << b(); // restore old flags and fill character os.flags(flags); os.fill(fill); return os; } template template std::basic_istream& uniform_int_distribution::stream_in(std::basic_istream& is) { typedef std::basic_istream istream_type; typedef typename istream_type::ios_base ios_base; // save old flags const typename ios_base::fmtflags flags = is.flags(); is.flags(ios_base::skipws); is >> m_param.first >> m_param.second; // restore old flags is.flags(flags); return is; } template _CCCL_HOST_DEVICE bool operator==(const uniform_int_distribution& lhs, const uniform_int_distribution& rhs) { return thrust::random::detail::random_core_access::equal(lhs, rhs); } template _CCCL_HOST_DEVICE bool operator!=(const uniform_int_distribution& lhs, const uniform_int_distribution& rhs) { return !(lhs == rhs); } template std::basic_ostream& operator<<(std::basic_ostream& os, const uniform_int_distribution& d) { return thrust::random::detail::random_core_access::stream_out(os, d); } template std::basic_istream& operator>>(std::basic_istream& is, uniform_int_distribution& d) { return thrust::random::detail::random_core_access::stream_in(is, d); } } // namespace random THRUST_NAMESPACE_END