I have a Monte Carlo simulation program (Ising model) in which I wanted to
use Boost.Random number generators instead of one from *Numerical Recipes. *It
works when the generator is called only in main, but when I use it inside
function, I have "identifier not found". I don't want to initialize them
each time I call a function because this will repeat the sequences of
random number. Is there a way to overcome it?
I don't know much about typedefs and templates, so my initialization of the
generator is just copied from some tutorial and changed in, I hope,
appropriate places.
Here is the example code (I wouldn't paste original program, it's too long,
but here is the same problem)
// random_test.cpp : Defines the entry point for the console application.
//
#include "stdafx.h"
#include "iostream"
#include <ctime>
#include
using namespace std;
namespace boost {
namespace random {
template
class lagged_fibonacci_engine;
template
class lagged_fibonacci_01_engine;
typedef lagged_fibonacci_01_engine< double, 48, 607, 273 >
lagged_fibonacci607;
typedef lagged_fibonacci_01_engine< double, 48, 1279, 418 >
lagged_fibonacci1279;
typedef lagged_fibonacci_01_engine< double, 48, 2281, 1252 >
lagged_fibonacci2281;
typedef lagged_fibonacci_01_engine< double, 48, 3217, 576 >
lagged_fibonacci3217;
typedef lagged_fibonacci_01_engine< double, 48, 4423, 2098 >
lagged_fibonacci4423;
typedef lagged_fibonacci_01_engine< double, 48, 9689, 5502 >
lagged_fibonacci9689;
typedef lagged_fibonacci_01_engine< double, 48, 19937, 9842 >
lagged_fibonacci19937;
typedef lagged_fibonacci_01_engine< double, 48, 23209, 13470 >
lagged_fibonacci23209;
typedef lagged_fibonacci_01_engine< double, 48, 44497, 21034 >
lagged_fibonacci44497;
}
}
void generate();
int _tmain(int argc, _TCHAR* argv[])
{const int rangeMin = 1;
const int rangeMax = 10;
typedef boost::uniform_real<> NumberDistribution;
typedef boost::random::lagged_fibonacci9689 RandomNumberGenerator;
typedef boost::variate_generator Generator;
NumberDistribution distribution(rangeMin, rangeMax);
RandomNumberGenerator generator;
Generator numberGenerator(generator, distribution);
generator.seed(time(0)); // seed with the current time
generate();
return 0;
}
void generate(){cout << numberGenerator() << endl;}
