FM / random

/usr/include/boost/random

Name	Size	Mode	Actions
detail/	-	0755	rm
additive_combine.hpp	8826	0644	edit dl rm
bernoulli_distribution.hpp	5608	0644	edit dl rm
beta_distribution.hpp	5900	0644	edit dl rm
binomial_distribution.hpp	12985	0644	edit dl rm
cauchy_distribution.hpp	6528	0644	edit dl rm
chi_squared_distribution.hpp	6322	0644	edit dl rm
discard_block.hpp	7991	0644	edit dl rm
discrete_distribution.hpp	21137	0644	edit dl rm
exponential_distribution.hpp	20895	0644	edit dl rm
extreme_value_distribution.hpp	5611	0644	edit dl rm
faure.hpp	11474	0644	edit dl rm
fisher_f_distribution.hpp	5924	0644	edit dl rm
gamma_distribution.hpp	9166	0644	edit dl rm
generate_canonical.hpp	2954	0644	edit dl rm
geometric_distribution.hpp	7339	0644	edit dl rm
hyperexponential_distribution.hpp	36980	0644	edit dl rm
independent_bits.hpp	9117	0644	edit dl rm
inversive_congruential.hpp	9550	0644	edit dl rm
lagged_fibonacci.hpp	18256	0644	edit dl rm
laplace_distribution.hpp	5743	0644	edit dl rm
linear_congruential.hpp	16214	0644	edit dl rm
linear_feedback_shift.hpp	7162	0644	edit dl rm
lognormal_distribution.hpp	8051	0644	edit dl rm
mersenne_twister.hpp	24094	0644	edit dl rm
negative_binomial_distribution.hpp	6911	0644	edit dl rm
niederreiter_base2.hpp	11736	0644	edit dl rm
non_central_chi_squared_distribution.hpp	7724	0644	edit dl rm
normal_distribution.hpp	17441	0644	edit dl rm
piecewise_constant_distribution.hpp	17186	0644	edit dl rm
piecewise_linear_distribution.hpp	18802	0644	edit dl rm
poisson_distribution.hpp	10145	0644	edit dl rm
random_device.hpp	5032	0644	edit dl rm
random_number_generator.hpp	1958	0644	edit dl rm
ranlux.hpp	3292	0644	edit dl rm
seed_seq.hpp	3792	0644	edit dl rm
shuffle_order.hpp	9181	0644	edit dl rm
shuffle_output.hpp	1352	0644	edit dl rm
sobol.hpp	7704	0644	edit dl rm
student_t_distribution.hpp	5633	0644	edit dl rm
subtract_with_carry.hpp	21590	0644	edit dl rm
taus88.hpp	1164	0644	edit dl rm
traits.hpp	4165	0644	edit dl rm
triangle_distribution.hpp	7033	0644	edit dl rm
uniform_01.hpp	7538	0644	edit dl rm
uniform_int.hpp	2909	0644	edit dl rm
uniform_int_distribution.hpp	16199	0644	edit dl rm
uniform_on_sphere.hpp	8812	0644	edit dl rm
uniform_real.hpp	2449	0644	edit dl rm
uniform_real_distribution.hpp	7726	0644	edit dl rm
uniform_smallint.hpp	12230	0644	edit dl rm
variate_generator.hpp	3708	0644	edit dl rm
weibull_distribution.hpp	5524	0644	edit dl rm
xor_combine.hpp	7049	0644	edit dl rm

Edit: /usr/include/boost/random/generate_canonical.hpp (2954B)

#ifndef BOOST_RANDOM_GENERATE_CANONICAL_HPP
#define BOOST_RANDOM_GENERATE_CANONICAL_HPP

#include <algorithm>
#include <boost/assert.hpp>
#include <boost/config/no_tr1/cmath.hpp>
#include <boost/limits.hpp>
#include <boost/type_traits/is_integral.hpp>
#include <boost/mpl/bool.hpp>
#include <boost/random/detail/signed_unsigned_tools.hpp>
#include <boost/random/detail/generator_bits.hpp>

namespace boost {
namespace random {

namespace detail {

template<class RealType, std::size_t bits, class URNG>
RealType generate_canonical_impl(URNG& g, boost::mpl::true_ /*is_integral*/)
{
    using std::pow;
    typedef typename URNG::result_type base_result;
    std::size_t digits = std::numeric_limits<RealType>::digits;
    RealType R = RealType((g.max)()) - RealType((g.min)()) + 1;
    RealType mult = R;
    RealType limit =
        pow(RealType(2),
            RealType((std::min)(static_cast<std::size_t>(bits), digits)));
    RealType S = RealType(detail::subtract<base_result>()(g(), (g.min)()));
    while(mult < limit) {
        RealType inc = RealType(detail::subtract<base_result>()(g(), (g.min)()));
        S += inc * mult;
        mult *= R;
    }
    return S / mult;
}

template<class RealType, std::size_t bits, class URNG>
RealType generate_canonical_impl(URNG& g, boost::mpl::false_ /*is_integral*/)
{
    using std::pow;
    using std::floor;
    BOOST_ASSERT((g.min)() == 0);
    BOOST_ASSERT((g.max)() == 1);
    std::size_t digits = std::numeric_limits<RealType>::digits;
    std::size_t engine_bits = detail::generator_bits<URNG>::value();
    std::size_t b = (std::min)(bits, digits);
    RealType R = pow(RealType(2), RealType(engine_bits));
    RealType mult = R;
    RealType limit = pow(RealType(2), RealType(b));
    RealType S = RealType(g() - (g.min)());
    while(mult < limit) {
        RealType inc(floor((RealType(g()) - RealType((g.min)())) * R));
        S += inc * mult;
        mult *= R;
    }
    return S / mult;
}

}

/**
 * Returns a value uniformly distributed in the range [0, 1)
 * with at least @c bits random bits.
 */
template<class RealType, std::size_t bits, class URNG>
RealType generate_canonical(URNG& g)
{
    RealType result = detail::generate_canonical_impl<RealType, bits>(
        g, boost::random::traits::is_integral<typename URNG::result_type>());
    BOOST_ASSERT(result >= 0);
    BOOST_ASSERT(result <= 1);
    if(result == 1) {
        result -= std::numeric_limits<RealType>::epsilon() / 2;
        BOOST_ASSERT(result != 1);
    }
    return result;
}

} // namespace random
} // namespace boost

#endif // BOOST_RANDOM_GENERATE_CANONICAL_HPP

/* boost random/generate_canonical.hpp header file * * Copyright Steven Watanabe 2011 * Distributed under the Boost Software License, Version 1.0. (See * accompanying file LICENSE_1_0.txt or copy at * http://www.boost.org/LICENSE_1_0.txt) * * See http://www.boost.org for most recent version including documentation. * * $Id$ * */ #ifndef BOOST_RANDOM_GENERATE_CANONICAL_HPP #define BOOST_RANDOM_GENERATE_CANONICAL_HPP #include #include #include #include #include #include #include #include namespace boost { namespace random { namespace detail { template RealType generate_canonical_impl(URNG& g, boost::mpl::true_ /*is_integral*/) { using std::pow; typedef typename URNG::result_type base_result; std::size_t digits = std::numeric_limits::digits; RealType R = RealType((g.max)()) - RealType((g.min)()) + 1; RealType mult = R; RealType limit = pow(RealType(2), RealType((std::min)(static_cast(bits), digits))); RealType S = RealType(detail::subtract()(g(), (g.min)())); while(mult < limit) { RealType inc = RealType(detail::subtract()(g(), (g.min)())); S += inc * mult; mult *= R; } return S / mult; } template RealType generate_canonical_impl(URNG& g, boost::mpl::false_ /*is_integral*/) { using std::pow; using std::floor; BOOST_ASSERT((g.min)() == 0); BOOST_ASSERT((g.max)() == 1); std::size_t digits = std::numeric_limits::digits; std::size_t engine_bits = detail::generator_bits::value(); std::size_t b = (std::min)(bits, digits); RealType R = pow(RealType(2), RealType(engine_bits)); RealType mult = R; RealType limit = pow(RealType(2), RealType(b)); RealType S = RealType(g() - (g.min)()); while(mult < limit) { RealType inc(floor((RealType(g()) - RealType((g.min)())) * R)); S += inc * mult; mult *= R; } return S / mult; } } /** * Returns a value uniformly distributed in the range [0, 1) * with at least @c bits random bits. */ template RealType generate_canonical(URNG& g) { RealType result = detail::generate_canonical_impl( g, boost::random::traits::is_integral()); BOOST_ASSERT(result >= 0); BOOST_ASSERT(result <= 1); if(result == 1) { result -= std::numeric_limits::epsilon() / 2; BOOST_ASSERT(result != 1); } return result; } } // namespace random } // namespace boost #endif // BOOST_RANDOM_GENERATE_CANONICAL_HPP