|
Packit |
8dc392 |
///////////////////////////////////////////////////////////////////////////
|
|
Packit |
8dc392 |
//
|
|
Packit |
8dc392 |
// Copyright (c) 2002-2012, Industrial Light & Magic, a division of Lucas
|
|
Packit |
8dc392 |
// Digital Ltd. LLC
|
|
Packit |
8dc392 |
//
|
|
Packit |
8dc392 |
// All rights reserved.
|
|
Packit |
8dc392 |
//
|
|
Packit |
8dc392 |
// Redistribution and use in source and binary forms, with or without
|
|
Packit |
8dc392 |
// modification, are permitted provided that the following conditions are
|
|
Packit |
8dc392 |
// met:
|
|
Packit |
8dc392 |
// * Redistributions of source code must retain the above copyright
|
|
Packit |
8dc392 |
// notice, this list of conditions and the following disclaimer.
|
|
Packit |
8dc392 |
// * Redistributions in binary form must reproduce the above
|
|
Packit |
8dc392 |
// copyright notice, this list of conditions and the following disclaimer
|
|
Packit |
8dc392 |
// in the documentation and/or other materials provided with the
|
|
Packit |
8dc392 |
// distribution.
|
|
Packit |
8dc392 |
// * Neither the name of Industrial Light & Magic nor the names of
|
|
Packit |
8dc392 |
// its contributors may be used to endorse or promote products derived
|
|
Packit |
8dc392 |
// from this software without specific prior written permission.
|
|
Packit |
8dc392 |
//
|
|
Packit |
8dc392 |
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
|
|
Packit |
8dc392 |
// "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
|
|
Packit |
8dc392 |
// LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
|
|
Packit |
8dc392 |
// A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
|
|
Packit |
8dc392 |
// OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
|
|
Packit |
8dc392 |
// SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
|
|
Packit |
8dc392 |
// LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
|
|
Packit |
8dc392 |
// DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
|
|
Packit |
8dc392 |
// THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
|
|
Packit |
8dc392 |
// (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
|
|
Packit |
8dc392 |
// OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
|
Packit |
8dc392 |
//
|
|
Packit |
8dc392 |
///////////////////////////////////////////////////////////////////////////
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
#ifndef INCLUDED_IMATHRANDOM_H
|
|
Packit |
8dc392 |
#define INCLUDED_IMATHRANDOM_H
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//-----------------------------------------------------------------------------
|
|
Packit |
8dc392 |
//
|
|
Packit |
8dc392 |
// Generators for uniformly distributed pseudo-random numbers and
|
|
Packit |
8dc392 |
// functions that use those generators to generate numbers with
|
|
Packit |
8dc392 |
// non-uniform distributions:
|
|
Packit |
8dc392 |
//
|
|
Packit |
8dc392 |
// class Rand32
|
|
Packit |
8dc392 |
// class Rand48
|
|
Packit |
8dc392 |
// solidSphereRand()
|
|
Packit |
8dc392 |
// hollowSphereRand()
|
|
Packit |
8dc392 |
// gaussRand()
|
|
Packit |
8dc392 |
// gaussSphereRand()
|
|
Packit |
8dc392 |
//
|
|
Packit |
8dc392 |
// Note: class Rand48() calls erand48() and nrand48(), which are not
|
|
Packit |
8dc392 |
// available on all operating systems. For compatibility we include
|
|
Packit |
8dc392 |
// our own versions of erand48() and nrand48(). Our functions have
|
|
Packit |
8dc392 |
// been reverse-engineered from the corresponding Unix/Linux man page.
|
|
Packit |
8dc392 |
//
|
|
Packit |
8dc392 |
//-----------------------------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
#include "ImathNamespace.h"
|
|
Packit |
8dc392 |
#include "ImathExport.h"
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
#include <stdlib.h>
|
|
Packit |
8dc392 |
#include <math.h>
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
IMATH_INTERNAL_NAMESPACE_HEADER_ENTER
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//-----------------------------------------------
|
|
Packit |
8dc392 |
// Fast random-number generator that generates
|
|
Packit |
8dc392 |
// a uniformly distributed sequence with a period
|
|
Packit |
8dc392 |
// length of 2^32.
|
|
Packit |
8dc392 |
//-----------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
class IMATH_EXPORT Rand32
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
public:
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//------------
|
|
Packit |
8dc392 |
// Constructor
|
|
Packit |
8dc392 |
//------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
Rand32 (unsigned long int seed = 0);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//--------------------------------
|
|
Packit |
8dc392 |
// Re-initialize with a given seed
|
|
Packit |
8dc392 |
//--------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
void init (unsigned long int seed);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//----------------------------------------------------------
|
|
Packit |
8dc392 |
// Get the next value in the sequence (range: [false, true])
|
|
Packit |
8dc392 |
//----------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
bool nextb ();
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//---------------------------------------------------------------
|
|
Packit |
8dc392 |
// Get the next value in the sequence (range: [0 ... 0xffffffff])
|
|
Packit |
8dc392 |
//---------------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
unsigned long int nexti ();
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//------------------------------------------------------
|
|
Packit |
8dc392 |
// Get the next value in the sequence (range: [0 ... 1[)
|
|
Packit |
8dc392 |
//------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
float nextf ();
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//-------------------------------------------------------------------
|
|
Packit |
8dc392 |
// Get the next value in the sequence (range [rangeMin ... rangeMax[)
|
|
Packit |
8dc392 |
//-------------------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
float nextf (float rangeMin, float rangeMax);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
private:
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
void next ();
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
unsigned long int _state;
|
|
Packit |
8dc392 |
};
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//--------------------------------------------------------
|
|
Packit |
8dc392 |
// Random-number generator based on the C Standard Library
|
|
Packit |
8dc392 |
// functions erand48(), nrand48() & company; generates a
|
|
Packit |
8dc392 |
// uniformly distributed sequence.
|
|
Packit |
8dc392 |
//--------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
class Rand48
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
public:
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//------------
|
|
Packit |
8dc392 |
// Constructor
|
|
Packit |
8dc392 |
//------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
Rand48 (unsigned long int seed = 0);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//--------------------------------
|
|
Packit |
8dc392 |
// Re-initialize with a given seed
|
|
Packit |
8dc392 |
//--------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
void init (unsigned long int seed);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//----------------------------------------------------------
|
|
Packit |
8dc392 |
// Get the next value in the sequence (range: [false, true])
|
|
Packit |
8dc392 |
//----------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
bool nextb ();
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//---------------------------------------------------------------
|
|
Packit |
8dc392 |
// Get the next value in the sequence (range: [0 ... 0x7fffffff])
|
|
Packit |
8dc392 |
//---------------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
long int nexti ();
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//------------------------------------------------------
|
|
Packit |
8dc392 |
// Get the next value in the sequence (range: [0 ... 1[)
|
|
Packit |
8dc392 |
//------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
double nextf ();
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//-------------------------------------------------------------------
|
|
Packit |
8dc392 |
// Get the next value in the sequence (range [rangeMin ... rangeMax[)
|
|
Packit |
8dc392 |
//-------------------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
double nextf (double rangeMin, double rangeMax);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
private:
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
unsigned short int _state[3];
|
|
Packit |
8dc392 |
};
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//------------------------------------------------------------
|
|
Packit |
8dc392 |
// Return random points uniformly distributed in a sphere with
|
|
Packit |
8dc392 |
// radius 1 around the origin (distance from origin <= 1).
|
|
Packit |
8dc392 |
//------------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
template <class Vec, class Rand>
|
|
Packit |
8dc392 |
Vec
|
|
Packit |
8dc392 |
solidSphereRand (Rand &rand;;
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//-------------------------------------------------------------
|
|
Packit |
8dc392 |
// Return random points uniformly distributed on the surface of
|
|
Packit |
8dc392 |
// a sphere with radius 1 around the origin.
|
|
Packit |
8dc392 |
//-------------------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
template <class Vec, class Rand>
|
|
Packit |
8dc392 |
Vec
|
|
Packit |
8dc392 |
hollowSphereRand (Rand &rand;;
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//-----------------------------------------------
|
|
Packit |
8dc392 |
// Return random numbers with a normal (Gaussian)
|
|
Packit |
8dc392 |
// distribution with zero mean and unit variance.
|
|
Packit |
8dc392 |
//-----------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
template <class Rand>
|
|
Packit |
8dc392 |
float
|
|
Packit |
8dc392 |
gaussRand (Rand &rand;;
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//----------------------------------------------------
|
|
Packit |
8dc392 |
// Return random points whose distance from the origin
|
|
Packit |
8dc392 |
// has a normal (Gaussian) distribution with zero mean
|
|
Packit |
8dc392 |
// and unit variance.
|
|
Packit |
8dc392 |
//----------------------------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
template <class Vec, class Rand>
|
|
Packit |
8dc392 |
Vec
|
|
Packit |
8dc392 |
gaussSphereRand (Rand &rand;;
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//---------------------------------
|
|
Packit |
8dc392 |
// erand48(), nrand48() and friends
|
|
Packit |
8dc392 |
//---------------------------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
IMATH_EXPORT double erand48 (unsigned short state[3]);
|
|
Packit |
8dc392 |
IMATH_EXPORT double drand48 ();
|
|
Packit |
8dc392 |
IMATH_EXPORT long int nrand48 (unsigned short state[3]);
|
|
Packit |
8dc392 |
IMATH_EXPORT long int lrand48 ();
|
|
Packit |
8dc392 |
IMATH_EXPORT void srand48 (long int seed);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
//---------------
|
|
Packit |
8dc392 |
// Implementation
|
|
Packit |
8dc392 |
//---------------
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline void
|
|
Packit |
8dc392 |
Rand32::init (unsigned long int seed)
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
_state = (seed * 0xa5a573a5L) ^ 0x5a5a5a5aL;
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline
|
|
Packit |
8dc392 |
Rand32::Rand32 (unsigned long int seed)
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
init (seed);
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline void
|
|
Packit |
8dc392 |
Rand32::next ()
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
_state = 1664525L * _state + 1013904223L;
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline bool
|
|
Packit |
8dc392 |
Rand32::nextb ()
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
next ();
|
|
Packit |
8dc392 |
// Return the 31st (most significant) bit, by and-ing with 2 ^ 31.
|
|
Packit |
8dc392 |
return !!(_state & 2147483648UL);
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline unsigned long int
|
|
Packit |
8dc392 |
Rand32::nexti ()
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
next ();
|
|
Packit |
8dc392 |
return _state & 0xffffffff;
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline float
|
|
Packit |
8dc392 |
Rand32::nextf (float rangeMin, float rangeMax)
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
float f = nextf();
|
|
Packit |
8dc392 |
return rangeMin * (1 - f) + rangeMax * f;
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline void
|
|
Packit |
8dc392 |
Rand48::init (unsigned long int seed)
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
seed = (seed * 0xa5a573a5L) ^ 0x5a5a5a5aL;
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
_state[0] = (unsigned short int) (seed & 0xFFFF);
|
|
Packit |
8dc392 |
_state[1] = (unsigned short int) ((seed >> 16) & 0xFFFF);
|
|
Packit |
8dc392 |
_state[2] = (unsigned short int) (seed & 0xFFFF);
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline
|
|
Packit |
8dc392 |
Rand48::Rand48 (unsigned long int seed)
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
init (seed);
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline bool
|
|
Packit |
8dc392 |
Rand48::nextb ()
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
return nrand48 (_state) & 1;
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline long int
|
|
Packit |
8dc392 |
Rand48::nexti ()
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
return nrand48 (_state);
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline double
|
|
Packit |
8dc392 |
Rand48::nextf ()
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
return erand48 (_state);
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
inline double
|
|
Packit |
8dc392 |
Rand48::nextf (double rangeMin, double rangeMax)
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
double f = nextf();
|
|
Packit |
8dc392 |
return rangeMin * (1 - f) + rangeMax * f;
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
template <class Vec, class Rand>
|
|
Packit |
8dc392 |
Vec
|
|
Packit |
8dc392 |
solidSphereRand (Rand &rand)
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
Vec v;
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
do
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
for (unsigned int i = 0; i < Vec::dimensions(); i++)
|
|
Packit |
8dc392 |
v[i] = (typename Vec::BaseType) rand.nextf (-1, 1);
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
while (v.length2() > 1);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
return v;
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
template <class Vec, class Rand>
|
|
Packit |
8dc392 |
Vec
|
|
Packit |
8dc392 |
hollowSphereRand (Rand &rand)
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
Vec v;
|
|
Packit |
8dc392 |
typename Vec::BaseType length;
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
do
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
for (unsigned int i = 0; i < Vec::dimensions(); i++)
|
|
Packit |
8dc392 |
v[i] = (typename Vec::BaseType) rand.nextf (-1, 1);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
length = v.length();
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
while (length > 1 || length == 0);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
return v / length;
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
template <class Rand>
|
|
Packit |
8dc392 |
float
|
|
Packit |
8dc392 |
gaussRand (Rand &rand)
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
float x; // Note: to avoid numerical problems with very small
|
|
Packit |
8dc392 |
float y; // numbers, we make these variables singe-precision
|
|
Packit |
8dc392 |
float length2; // floats, but later we call the double-precision log()
|
|
Packit |
8dc392 |
// and sqrt() functions instead of logf() and sqrtf().
|
|
Packit |
8dc392 |
do
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
x = float (rand.nextf (-1, 1));
|
|
Packit |
8dc392 |
y = float (rand.nextf (-1, 1));
|
|
Packit |
8dc392 |
length2 = x * x + y * y;
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
while (length2 >= 1 || length2 == 0);
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
return x * sqrt (-2 * log (double (length2)) / length2);
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
template <class Vec, class Rand>
|
|
Packit |
8dc392 |
Vec
|
|
Packit |
8dc392 |
gaussSphereRand (Rand &rand)
|
|
Packit |
8dc392 |
{
|
|
Packit |
8dc392 |
return hollowSphereRand <Vec> (rand) * gaussRand (rand);
|
|
Packit |
8dc392 |
}
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
IMATH_INTERNAL_NAMESPACE_HEADER_EXIT
|
|
Packit |
8dc392 |
|
|
Packit |
8dc392 |
#endif // INCLUDED_IMATHRANDOM_H
|