/usr/include/linbox/randiter/mersenne-twister.h is in liblinbox-dev 1.1.6~rc0-4.1.
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 66 67 68 69 70 71 72 73 | /* -*- mode: C++; tab-width: 8; indent-tabs-mode: t; c-basic-offset: 8 -*- */
/* linbox/randiter/mersenne-twister.h
* Copyright (C) 2002 Bradford Hovinen
*
* Written by Bradford Hovinen <hovinen@cis.udel.edu>
*
* ------------------------------------
*
* This file is part of LinBox, licensed under the GNU Lesser General Public
* License. See COPYING for more information.
*
* ------------------------------------
*
* Header file for implementation of the Mersenne twister pseudo-random number
* generator, crated by Makoto Matsumoto and Takuji Nishimura. Further
* information on the Mersenne Twister can be found at
* http://www.math.keio.ac.jp/~matsumoto/emt.html [ dead link. -bds 2007feb]
*
* This forms the basic underlying algorithm for most psuedo-random number
* generation in LinBox.
*
* N.B. This module is tested in test-modular under the random number
* generator tests.
*/
#ifndef __MERSENNE_TWISTER_H
#define __MERSENNE_TWISTER_H
#include <vector>
#include "linbox/integer.h"
namespace LinBox
{
class MersenneTwister
{
public:
MersenneTwister (uint32 seed = 0);
uint32 reload ();
uint32 randomInt ();
uint32 randomInt () const
{ return const_cast<MersenneTwister&>(*this).randomInt();}
uint32 randomIntRange (uint32 start, uint32 end);
uint32 randomIntRange (uint32 start, uint32 end) const
{ return const_cast<MersenneTwister&>(*this).randomIntRange(start,end); }
double randomDouble ();
double randomDouble () const
{ return const_cast<MersenneTwister&>(*this).randomDouble(); }
double randomDoubleRange (double start, double end)
{ return randomDouble () * (end - start) + start; }
double randomDoubleRange (double start, double end) const
{ return const_cast<MersenneTwister&>(*this).randomDoubleRange(start,end); }
void setSeed (uint32 seed);
private:
std::vector<uint32> _state;
std::vector<uint32>::iterator _next;
int _left;
};
}
#ifdef LinBoxSrcOnly
#include <linbox/randiter/mersenne-twister.C>
#endif
#endif // __MERSENNE_TWISTER_H
|