RhodeCode

    ///

    /// @{

    /// \brief Default constructor

    ///

    /// Constructor with constant seeding.

      return real<Number>() * (b - a) + a;

    }

    /// \brief Returns a random real number from the range [0, 1)

    ///

    /// It returns a random double from the range [0, 1).

      return _random_bits::IntConversion<Number, Word>::convert(core);

    }

    unsigned int uinteger() {

      return uinteger<unsigned int>();

    }

    ///\name Non-uniform distributions

    ///

    ///@{

    ///

    /// It returns a random bool with given probability of true result.

    bool boolean(double p) {

      return operator()() < p;

    }

    /// \note The Cartesian form of the Box-Muller

    /// transformation is used to generate a random normal distribution.

    double gauss()

      } while(S>=1);

      return std::sqrt(-2*std::log(S)/S)*V1;

    }

    /// \sa gauss()

    double gauss(double mean,double std_dev)

    {

    /// standard deviation. The return value can direcly be passed to

    /// lognormal().

    std::pair<double,double> lognormalParamsFromMD(double mean,

    {

      double fr=std_dev/mean;

      fr*=fr;

      return std::pair<double,double>(std::log(mean)-lg/2.0,std::sqrt(lg));

    }

    /// Lognormal distribution with given mean and standard deviation

    /// Lognormal distribution with given mean and standard deviation.

    ///

    double lognormalMD(double mean,double std_dev)

    {

      return lognormal(lognormalParamsFromMD(mean,std_dev));

    }

    /// Exponential distribution with given mean

    /// This function generates an exponential distribution random number

    ///\name Two dimensional distributions

    ///

    ///@{

    /// Uniform distribution on the full unit circle

      } while(V1*V1+V2*V2>=1);

      return dim2::Point<double>(V1,V2);

    }

    /// This function provides a turning symmetric two-dimensional distribution.

    /// Both coordinates are of standard normal distribution, but they are not