21 #include "CLHEP/Random/defs.h"
22 #include "CLHEP/Random/NonRandomEngine.h"
23 #include "CLHEP/Random/engineIDulong.h"
24 #include "CLHEP/Random/DoubConv.hh"
37 sequenceHasBeenSet(false),
38 intervalHasBeenSet(false) ,
41 randomInterval(0.1) { }
54 for (
int i=0; i<
n; i++) sequence.push_back(*s++);
55 assert (sequence.size() == (
unsigned int)n);
57 sequenceHasBeenSet=
true;
64 intervalHasBeenSet=
true;
70 if (sequenceHasBeenSet) {
71 double v = sequence[nInSeq++];
72 if (nInSeq >= sequence.size() ) sequenceHasBeenSet =
false;
76 if ( !nextHasBeenSet ) {
78 <<
"Attempt to use NonRandomEngine without setting next random!\n";
82 double a = nextRandom;
83 nextHasBeenSet =
false;
85 if (intervalHasBeenSet) {
86 nextRandom += randomInterval;
87 if ( nextRandom >= 1 ) nextRandom -= 1.0;
88 nextHasBeenSet =
true;
96 for (
int i = 0; i < size; ++i) {
102 std::string beginMarker =
"NonRandomEngine-begin";
103 os << beginMarker <<
"\nUvec\n";
104 std::vector<unsigned long> v =
put();
105 for (
unsigned int i=0; i<v.size(); ++i) {
110 std::string endMarker =
"NonRandomEngine-end";
111 int pr = os.precision(20);
112 os <<
" " << beginMarker <<
"\n";
113 os << nextHasBeenSet <<
" ";
114 os << sequenceHasBeenSet <<
" ";
115 os << intervalHasBeenSet <<
"\n";
116 os << nextRandom <<
" " << nInSeq <<
" " << randomInterval <<
"\n";
117 os << sequence.size() <<
"\n";
118 for (
unsigned int i = 0; i < sequence.size(); ++i) {
119 os << sequence[i] <<
"\n";
121 os << endMarker <<
"\n ";
128 std::vector<unsigned long> v;
129 v.push_back (engineIDulong<NonRandomEngine>());
130 std::vector<unsigned long> t;
131 v.push_back(static_cast<unsigned long>(nextHasBeenSet));
132 v.push_back(static_cast<unsigned long>(sequenceHasBeenSet));
133 v.push_back(static_cast<unsigned long>(intervalHasBeenSet));
135 v.push_back(t[0]); v.push_back(t[1]);
136 v.push_back(static_cast<unsigned long>(nInSeq));
138 v.push_back(t[0]); v.push_back(t[1]);
139 v.push_back(static_cast<unsigned long>(sequence.size()));
140 for (
unsigned int i=0; i<sequence.size(); ++i) {
142 v.push_back(t[0]); v.push_back(t[1]);
148 std::string beginMarker =
"NonRandomEngine-begin";
150 if (beginMarker !=
"NonRandomEngine-begin") {
151 is.clear(std::ios::badbit | is.rdstate());
152 std::cerr <<
"\nInput mispositioned or"
153 <<
"\nNonRandomEngine state description missing or"
154 <<
"\nwrong engine type found.\n";
161 return "NonRandomEngine-begin";
166 std::vector<unsigned long> v;
167 unsigned long uu = 99999;
168 unsigned long ssiz = 0;
172 for (
unsigned int istart=0; istart < 10; ++istart) {
175 is.clear(std::ios::badbit | is.rdstate());
176 std::cout <<
"istart = " << istart <<
"\n";
178 <<
"\nNonRandomEngine state (vector) description has no sequence size."
179 <<
"\ngetState() has failed."
180 <<
"\nInput stream is probably mispositioned now." << std::endl;
185 std::cout <<
"v[" << istart <<
"] = " << uu <<
"\n";
187 if (istart==9) ssiz = uu;
189 for (
unsigned int ivec=0; ivec < 2*ssiz; ++ivec) {
192 is.clear(std::ios::badbit | is.rdstate());
193 std::cerr <<
"\nNonRandomEngine state (vector) description improper."
194 <<
"\ngetState() has failed."
195 <<
"\nInput stream is probably mispositioned now." << std::endl;
200 std::cout <<
"v[" << v.size()-1 <<
"] = " << uu <<
"\n";
209 std::string endMarker =
"NonRandomEngine-end";
210 is >> sequenceHasBeenSet >> intervalHasBeenSet;
211 is >> nextRandom >> nInSeq >> randomInterval;
212 unsigned int seqSize;
216 for (
unsigned int i = 0; i < seqSize; ++i) {
218 sequence.push_back(x);
221 if (endMarker !=
"NonRandomEngine-end") {
222 is.clear(std::ios::badbit | is.rdstate());
223 std::cerr <<
"\n NonRandomEngine state description incomplete."
224 <<
"\nInput stream is probably mispositioned now." << std::endl;
231 if ((v[0] & 0xffffffffUL) != engineIDulong<NonRandomEngine>()) {
233 "\nNonRandomEngine get:state vector has wrong ID word - state unchanged\n";
240 unsigned int seqSize = v[9];
241 if (v.size() != 2*seqSize + 10 ) {
243 "\nNonRandomEngine get:state vector has wrong length - state unchanged\n";
244 std::cerr <<
" (length = " << v.size()
245 <<
"; expected " << 2*seqSize + 10 <<
")\n";
248 std::vector<unsigned long> t(2);
249 nextHasBeenSet = (v[1]!=0);
250 sequenceHasBeenSet = (v[2]!=0);
251 intervalHasBeenSet = (v[3]!=0);
256 for (
unsigned int i=0; i<seqSize; ++i) {
257 t[0] = v[2*i+10]; t[1] = v[2*i+11];
bool possibleKeywordInput(IS &is, const std::string &key, T &t)
static std::string beginTag()
void setNextRandom(double r)
static double longs2double(const std::vector< unsigned long > &v)
void flatArray(const int size, double *vect)
virtual ~NonRandomEngine()
static std::vector< unsigned long > dto2longs(double d)
void setRandomSequence(double *s, int n)
void setRandomInterval(double x)
virtual std::istream & get(std::istream &is)
virtual std::istream & getState(std::istream &is)
std::vector< unsigned long > put() const