+typedef unsigned int seedType;
+typedef boost::random::mt19937 engine_type;
+typedef boost::random::uniform_01<> uniform_01_dist;
+typedef boost::random::gamma_distribution<> gamma_dist;
+typedef boost::random::variate_generator<engine_type&, uniform_01_dist> uniform_01_generator;
+typedef boost::random::variate_generator<engine_type&, gamma_dist> gamma_generator;
+
+class engineFactory {
+public:
+ static void init() { seedEngine = new engine_type(time(NULL)); }
+ static void init(seedType seed) { seedEngine = new engine_type(seed); }
+
+ static void finish() { if (seedEngine != NULL) delete seedEngine; }
+
+ static engine_type *new_engine() {
+ seedType seed;
+ static std::set<seedType> seedSet; // empty set of seeds
+ std::set<seedType>::iterator iter;
+
+ do {
+ seed = (*seedEngine)();
+ iter = seedSet.find(seed);
+ } while (iter != seedSet.end());
+ seedSet.insert(seed);
+
+ return new engine_type(seed);
+ }
+
+ private:
+ static engine_type *seedEngine;
+};
+
+engine_type* engineFactory::seedEngine = NULL;