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/uniform_on_sphere.hpp (8812B)

/* boost random/uniform_on_sphere.hpp header file
 *
 * Copyright Jens Maurer 2000-2001
 * 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$
 *
 * Revision history
 *  2001-02-18  moved to individual header files
 */

#ifndef BOOST_RANDOM_UNIFORM_ON_SPHERE_HPP
#define BOOST_RANDOM_UNIFORM_ON_SPHERE_HPP

#include <vector>
#include <algorithm>     // std::transform
#include <functional>    // std::bind2nd, std::divides
#include <boost/assert.hpp>
#include <boost/random/detail/config.hpp>
#include <boost/random/detail/operators.hpp>
#include <boost/random/normal_distribution.hpp>

namespace boost {
namespace random {

/**
 * Instantiations of class template uniform_on_sphere model a
 * \random_distribution. Such a distribution produces random
 * numbers uniformly distributed on the unit sphere of arbitrary
 * dimension @c dim. The @c Cont template parameter must be a STL-like
 * container type with begin and end operations returning non-const
 * ForwardIterators of type @c Cont::iterator. 
 */
template<class RealType = double, class Cont = std::vector<RealType> >
class uniform_on_sphere
{
public:
    typedef RealType input_type;
    typedef Cont result_type;

class param_type
    {
    public:

typedef uniform_on_sphere distribution_type;

/**
         * Constructs the parameters of a uniform_on_sphere
         * distribution, given the dimension of the sphere.
         */
        explicit param_type(int dim_arg = 2) : _dim(dim_arg)
        {
            BOOST_ASSERT(_dim >= 0);
        }

/** Returns the dimension of the sphere. */
        int dim() const { return _dim; }

/** Writes the parameters to a @c std::ostream. */
        BOOST_RANDOM_DETAIL_OSTREAM_OPERATOR(os, param_type, parm)
        {
            os << parm._dim;
            return os;
        }

/** Reads the parameters from a @c std::istream. */
        BOOST_RANDOM_DETAIL_ISTREAM_OPERATOR(is, param_type, parm)
        {
            is >> parm._dim;
            return is;
        }

/** Returns true if the two sets of parameters are equal. */
        BOOST_RANDOM_DETAIL_EQUALITY_OPERATOR(param_type, lhs, rhs)
        { return lhs._dim == rhs._dim; }

/** Returns true if the two sets of parameters are different. */
        BOOST_RANDOM_DETAIL_INEQUALITY_OPERATOR(param_type)

private:
        int _dim;
    };

/**
     * Constructs a @c uniform_on_sphere distribution.
     * @c dim is the dimension of the sphere.
     *
     * Requires: dim >= 0
     */
    explicit uniform_on_sphere(int dim_arg = 2)
      : _container(dim_arg), _dim(dim_arg) { }

/**
     * Constructs a @c uniform_on_sphere distribution from its parameters.
     */
    explicit uniform_on_sphere(const param_type& parm)
      : _container(parm.dim()), _dim(parm.dim()) { }

// compiler-generated copy ctor and assignment operator are fine

/** Returns the dimension of the sphere. */
    int dim() const { return _dim; }

/** Returns the parameters of the distribution. */
    param_type param() const { return param_type(_dim); }
    /** Sets the parameters of the distribution. */
    void param(const param_type& parm)
    {
        _dim = parm.dim();
        _container.resize(_dim);
    }

/**
     * Returns the smallest value that the distribution can produce.
     * Note that this is required to approximate the standard library's
     * requirements.  The behavior is defined according to lexicographical
     * comparison so that for a container type of std::vector,
     * dist.min() <= x <= dist.max() where x is any value produced
     * by the distribution.
     */
    result_type min BOOST_PREVENT_MACRO_SUBSTITUTION () const
    {
        result_type result(_dim);
        if(_dim != 0) {
            result.front() = RealType(-1.0);
        }
        return result;
    }
    /**
     * Returns the largest value that the distribution can produce.
     * Note that this is required to approximate the standard library's
     * requirements.  The behavior is defined according to lexicographical
     * comparison so that for a container type of std::vector,
     * dist.min() <= x <= dist.max() where x is any value produced
     * by the distribution.
     */
    result_type max BOOST_PREVENT_MACRO_SUBSTITUTION () const
    {
        result_type result(_dim);
        if(_dim != 0) {
            result.front() = RealType(1.0);
        }
        return result;
    }

/**
     * Effects: Subsequent uses of the distribution do not depend
     * on values produced by any engine prior to invoking reset.
     */
    void reset() {}

/**
     * Returns a point uniformly distributed over the surface of
     * a sphere of dimension dim().
     */
    template<class Engine>
    const result_type & operator()(Engine& eng)
    {
        using std::sqrt;
        switch(_dim)
        {
        case 0: break;
        case 1:
            {
                if(uniform_01<RealType>()(eng) < 0.5) {
                    *_container.begin() = -1;
                } else {
                    *_container.begin() = 1;
                }
                break;
            }
        case 2:
            {
                uniform_01<RealType> uniform;
                RealType sqsum;
                RealType x, y;
                do {
                    x = uniform(eng) * 2 - 1;
                    y = uniform(eng) * 2 - 1;
                    sqsum = x*x + y*y;
                } while(sqsum == 0 || sqsum > 1);
                RealType mult = 1/sqrt(sqsum);
                typename Cont::iterator iter = _container.begin();
                *iter = x * mult;
                iter++;
                *iter = y * mult;
                break;
            }
        case 3:
            {
                uniform_01<RealType> uniform;
                RealType sqsum;
                RealType x, y;
                do {
                    x = uniform(eng) * 2 - 1;
                    y = uniform(eng) * 2 - 1;
                    sqsum = x*x + y*y;
                } while(sqsum > 1);
                RealType mult = 2 * sqrt(1 - sqsum);
                typename Cont::iterator iter = _container.begin();
                *iter = x * mult;
                ++iter;
                *iter = y * mult;
                ++iter;
                *iter = 2 * sqsum - 1;
                break;
            }
        default:
            {
                detail::unit_normal_distribution<RealType> normal;
                RealType sqsum;
                do {
                    sqsum = 0;
                    for(typename Cont::iterator it = _container.begin();
                        it != _container.end();
                        ++it) {
                        RealType val = normal(eng);
                        *it = val;
                        sqsum += val * val;
                    }
                } while(sqsum == 0);
                // for all i: result[i] /= sqrt(sqsum)
                RealType inverse_distance = 1 / sqrt(sqsum);
                for(typename Cont::iterator it = _container.begin();
                    it != _container.end();
                    ++it) {
                    *it *= inverse_distance;
                }
            }
        }
        return _container;
    }

/**
     * Returns a point uniformly distributed over the surface of
     * a sphere of dimension param.dim().
     */
    template<class Engine>
    result_type operator()(Engine& eng, const param_type& parm) const
    {
        return uniform_on_sphere(parm)(eng);
    }

/** Writes the distribution to a @c std::ostream. */
    BOOST_RANDOM_DETAIL_OSTREAM_OPERATOR(os, uniform_on_sphere, sd)
    {
        os << sd._dim;
        return os;
    }

/** Reads the distribution from a @c std::istream. */
    BOOST_RANDOM_DETAIL_ISTREAM_OPERATOR(is, uniform_on_sphere, sd)
    {
        is >> sd._dim;
        sd._container.resize(sd._dim);
        return is;
    }

/**
     * Returns true if the two distributions will produce identical
     * sequences of values, given equal generators.
     */
    BOOST_RANDOM_DETAIL_EQUALITY_OPERATOR(uniform_on_sphere, lhs, rhs)
    { return lhs._dim == rhs._dim; }

/**
     * Returns true if the two distributions may produce different
     * sequences of values, given equal generators.
     */
    BOOST_RANDOM_DETAIL_INEQUALITY_OPERATOR(uniform_on_sphere)

private:
    result_type _container;
    int _dim;
};

} // namespace random

using random::uniform_on_sphere;

} // namespace boost

#endif // BOOST_RANDOM_UNIFORM_ON_SPHERE_HPP

cmd: