44 |
|
#include "primitives/Molecule.hpp" |
45 |
|
#include "primitives/StuntDouble.hpp" |
46 |
|
|
47 |
+ |
#ifndef IS_MPI |
48 |
+ |
#include "math/SeqRandNumGen.hpp" |
49 |
+ |
#else |
50 |
+ |
#include "math/ParallelRandNumGen.hpp" |
51 |
+ |
#endif |
52 |
+ |
|
53 |
|
namespace oopse { |
54 |
|
|
55 |
< |
Velocitizer::Velocitizer(SimInfo* info) { |
55 |
> |
Velocitizer::Velocitizer(SimInfo* info) : info_(info) { |
56 |
|
|
57 |
|
int seedValue; |
58 |
|
Globals * simParams = info->getSimParams(); |
60 |
|
#ifndef IS_MPI |
61 |
|
if (simParams->haveSeed()) { |
62 |
|
seedValue = simParams->getSeed(); |
63 |
< |
randNumGen_ = new MTRand(seedValue); |
63 |
> |
randNumGen_ = new SeqRandNumGen(seedValue); |
64 |
|
}else { |
65 |
< |
randNumGen_ = new MTRand(); |
65 |
> |
randNumGen_ = new SeqRandNumGen(); |
66 |
|
} |
67 |
|
#else |
68 |
|
if (simParams->haveSeed()) { |