/usr/include/random/chisquare.h is in libblitz0-dev 1:0.10-3.3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 | // -*- C++ -*-
// $Id$
/*
* Chi^2 distribution
*
* This code has been adapted from RANDLIB.C 1.3, by
* Barry W. Brown, James Lovato, Kathy Russell, and John Venier.
* Code was originally by Ahrens and Dieter (see above).
*
* Adapter's notes:
*/
#ifndef BZ_RANDOM_CHISQUARE
#define BZ_RANDOM_CHISQUARE
#ifndef BZ_RANDOM_GAMMA
#include <random/gamma.h>
#endif
BZ_NAMESPACE(ranlib)
template<typename T = double, typename IRNG = defaultIRNG,
typename stateTag = defaultState>
class ChiSquare : public Gamma<T,IRNG,stateTag>
{
public:
typedef T T_numtype;
ChiSquare(T df)
: Gamma<T,IRNG,stateTag>(df/2.0) // isn't this redundant with setDF call?
{
setDF(df);
}
ChiSquare(T df, unsigned int i)
: Gamma<T,IRNG,stateTag>(df/2.0, i)
{
setDF(df);
}
void setDF(T _df)
{
BZPRECONDITION(_df > 0.0);
df = _df;
Gamma<T,IRNG,stateTag>::setMean(df/2.0);
}
T random()
{
return 2.0 * sgamma();
}
protected:
T sgamma()
{
return Gamma<T,IRNG,stateTag>::random();
}
T df;
};
BZ_NAMESPACE_END
#endif // BZ_RANDOM_CHISQUARE
|