| 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 |
|
|
