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