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