35#include "CLHEP/Random/defs.h"
36#include "CLHEP/Random/RanshiEngine.h"
37#include "CLHEP/Random/engineIDulong.h"
42static const int MarkerLen = 64;
47int RanshiEngine::numEngines = 0;
51 halfBuff(0), numFlats(0)
55 buffer[i] = (
unsigned int)(numEngines+19780503L*(i+1));
58 theSeed = numEngines+19780503L*++i;
59 redSpin = (
unsigned int)(
theSeed & 0xffffffff);
61 for( i = 0; i < 10000; ++i)
flat();
66 halfBuff(0), numFlats(0)
73 halfBuff(0), numFlats(0)
75 for (
int i = 0; i < numBuff; ++i) {
76 buffer[i] = (
unsigned int)seed&0xffffffff;
79 redSpin = (
unsigned int)(
theSeed & 0xffffffff);
81 for (j = 0; j < numBuff*20; ++j) {
88 halfBuff(0), numFlats(0)
91 while( i < numBuff ) {
92 buffer[i] = (
unsigned int)((rowIndex + (i+1)*(colIndex+8))&0xffffffff);
96 redSpin = colIndex & 0xffffffff;
97 for( i = 0; i < 100; ++i)
flat();
103 unsigned int redAngle = (((numBuff/2) - 1) & redSpin) + halfBuff;
104 unsigned int blkSpin = buffer[redAngle] & 0xffffffff;
105 unsigned int boostResult = blkSpin ^ redSpin;
107 buffer[redAngle] = ((blkSpin << 17) | (blkSpin >> (32-17))) ^ redSpin;
109 redSpin = (blkSpin + numFlats++) & 0xffffffff;
110 halfBuff = numBuff/2 - halfBuff;
118 for (
int i = 0; i < size; ++i) {
130 while (seeds[i] && i < numBuff) {
131 buffer[i] = (
unsigned int)seeds[i];
134 while (i < numBuff) {
135 buffer[i] = buffer[i-1];
139 redSpin = (
unsigned int)
theSeed;
145 std::ofstream outFile(filename, std::ios::out);
146 if (!outFile.bad()) {
148 std::vector<unsigned long> v =
put();
150 std::cout <<
"Result of v = put() is:\n";
152 for (
unsigned int i=0; i<v.size(); ++i) {
153 outFile << v[i] <<
"\n";
155 std::cout << v[i] <<
" ";
156 if (i%6==0) std::cout <<
"\n";
164 if (!outFile.bad()) {
165 outFile << std::setprecision(20) <<
theSeed << std::endl;
166 for (
int i = 0; i < numBuff; ++i) {
167 outFile << buffer[i] <<
" ";
169 outFile << redSpin <<
" " << numFlats <<
" " << halfBuff << std::endl;
175 std::ifstream inFile(filename, std::ios::in);
177 std::cerr <<
" -- Engine state remains unchanged\n";
181 std::vector<unsigned long> v;
183 for (
unsigned int ivec=0; ivec < VECTOR_STATE_SIZE; ++ivec) {
186 std::cout <<
"ivec = " << ivec <<
" xin = " << xin <<
" ";
187 if (ivec%3 == 0) std::cout <<
"\n";
190 inFile.clear(std::ios::badbit | inFile.rdstate());
191 std::cerr <<
"\nRanshiEngine state (vector) description improper."
192 <<
"\nrestoreStatus has failed."
193 <<
"\nInput stream is probably mispositioned now." << std::endl;
204 for (
int i = 0; i < numBuff; ++i) {
207 inFile >> redSpin >> numFlats >> halfBuff;
212 std::cout << std::setprecision(20) << std::endl;
213 std::cout <<
"----------- Ranshi engine status ----------" << std::endl;
214 std::cout <<
"Initial seed = " <<
theSeed << std::endl;
215 std::cout <<
"Current red spin = " << redSpin << std::endl;
216 std::cout <<
"Values produced = " << numFlats << std::endl;
217 std::cout <<
"Side of buffer = " << (halfBuff ?
"upper" :
"lower")
219 std::cout <<
"Current buffer = " << std::endl;
220 for (
int i = 0; i < numBuff; i+=4) {
221 std::cout << std::setw(10) << std::setiosflags(std::ios::right)
222 << buffer[i] << std::setw(11) << buffer[i+1] << std::setw(11)
223 << buffer[i+2] << std::setw(11) << buffer[i+3] << std::endl;
225 std::cout <<
"-------------------------------------------" << std::endl;
228RanshiEngine::operator float() {
229 unsigned int redAngle = (((numBuff/2) - 1) & redSpin) + halfBuff;
230 unsigned int blkSpin = buffer[redAngle] & 0xffffffff;
232 buffer[redAngle] = ((blkSpin << 17) | (blkSpin >> (32-17))) ^ redSpin;
234 redSpin = (blkSpin + numFlats++) & 0xffffffff;
235 halfBuff = numBuff/2 - halfBuff;
237 return float(blkSpin * twoToMinus_32());
240RanshiEngine::operator
unsigned int() {
241 unsigned int redAngle = (((numBuff/2) - 1) & redSpin) + halfBuff;
242 unsigned int blkSpin = buffer[redAngle] & 0xffffffff;
244 buffer[redAngle] = ((blkSpin << 17) | (blkSpin >> (32-17))) ^ redSpin;
246 redSpin = (blkSpin + numFlats++) & 0xffffffff;
247 halfBuff = numBuff/2 - halfBuff;
253 char beginMarker[] =
"RanshiEngine-begin";
254 os << beginMarker <<
"\nUvec\n";
255 std::vector<unsigned long> v =
put();
256 for (
unsigned int i=0; i<v.size(); ++i) {
261 char endMarker[] =
"RanshiEngine-end";
262 int pr=os.precision(20);
263 os <<
" " << beginMarker <<
" ";
266 for (
int i = 0; i < numBuff; ++i) {
267 os << buffer[i] <<
"\n";
269 os << redSpin <<
" " << numFlats <<
"\n" << halfBuff;
271 os <<
" " << endMarker <<
"\n";
278 std::vector<unsigned long> v;
279 v.push_back (engineIDulong<RanshiEngine>());
280 for (
int i = 0; i < numBuff; ++i) {
281 v.push_back(
static_cast<unsigned long>(buffer[i]));
283 v.push_back(
static_cast<unsigned long>(redSpin));
284 v.push_back(
static_cast<unsigned long>(numFlats));
285 v.push_back(
static_cast<unsigned long>(halfBuff));
290 char beginMarker [MarkerLen];
296 if (strcmp(beginMarker,
"RanshiEngine-begin")) {
297 is.clear(std::ios::badbit | is.rdstate());
298 std::cerr <<
"\nInput mispositioned or"
299 <<
"\nRanshiEngine state description missing or"
300 <<
"\nwrong engine type found." << std::endl;
307 return "RanshiEngine-begin";
312 std::vector<unsigned long> v;
314 for (
unsigned int ivec=0; ivec < VECTOR_STATE_SIZE; ++ivec) {
317 is.clear(std::ios::badbit | is.rdstate());
318 std::cerr <<
"\nRanshiEngine state (vector) description improper."
319 <<
"\ngetState() has failed."
320 <<
"\nInput stream is probably mispositioned now." << std::endl;
331 char endMarker [MarkerLen];
332 for (
int i = 0; i < numBuff; ++i) {
335 is >> redSpin >> numFlats >> halfBuff;
339 if (strcmp(endMarker,
"RanshiEngine-end")) {
340 is.clear(std::ios::badbit | is.rdstate());
341 std::cerr <<
"\nRanshiEngine state description incomplete."
342 <<
"\nInput stream is probably mispositioned now." << std::endl;
349 if ((v[0] & 0xffffffffUL) != engineIDulong<RanshiEngine>()) {
351 "\nRanshiEngine get:state vector has wrong ID word - state unchanged\n";
358 if (v.size() != VECTOR_STATE_SIZE ) {
360 "\nRanshiEngine get:state vector has wrong length - state unchanged\n";
363 for (
int i = 0; i < numBuff; ++i) {
366 redSpin = v[numBuff+1];
367 numFlats = v[numBuff+2];
368 halfBuff = v[numBuff+3];
static double twoToMinus_32()
static double twoToMinus_53()
static double nearlyTwoToMinus_54()
static bool checkFile(std::istream &file, const std::string &filename, const std::string &classname, const std::string &methodname)
virtual std::istream & getState(std::istream &is)
static std::string beginTag()
void flatArray(const int size, double *vect)
virtual std::istream & get(std::istream &is)
void saveStatus(const char filename[]="RanshiEngine.conf") const
std::vector< unsigned long > put() const
static std::string engineName()
void restoreStatus(const char filename[]="RanshiEngine.conf")
void setSeeds(const long *seeds, int)
void setSeed(long seed, int)
bool possibleKeywordInput(IS &is, const std::string &key, T &t)