57 |
|
#include "UseTheForce/ForceFieldFactory.hpp" |
58 |
|
#include "utils/simError.h" |
59 |
|
#include "utils/StringUtils.hpp" |
60 |
< |
#include "math/MersenneTwister.hpp" |
60 |
> |
#include "math/SeqRandNumGen.hpp" |
61 |
|
#ifdef IS_MPI |
62 |
|
#include "io/mpiBASS.h" |
63 |
+ |
#include "math/ParallelRandNumGen.hpp" |
64 |
|
#endif |
65 |
|
|
66 |
|
namespace oopse { |
255 |
|
|
256 |
|
int seedValue; |
257 |
|
Globals * simParams = info->getSimParams(); |
258 |
< |
MTRand* myRandom; //divide labor does not need Parallel random number generator |
258 |
> |
SeqRandNumGen* myRandom; //divide labor does not need Parallel random number generator |
259 |
|
if (simParams->haveSeed()) { |
260 |
|
seedValue = simParams->getSeed(); |
261 |
< |
myRandom = new MTRand(seedValue); |
261 |
> |
myRandom = new SeqRandNumGen(seedValue); |
262 |
|
}else { |
263 |
< |
myRandom = new MTRand(); |
263 |
> |
myRandom = new SeqRandNumGen(); |
264 |
|
} |
265 |
|
|
266 |
|
|