| 47 |
|
* @version 1.0 |
| 48 |
|
*/ |
| 49 |
|
|
| 50 |
– |
#include <sprng.h> |
| 51 |
– |
|
| 50 |
|
#include "brains/MoleculeCreator.hpp" |
| 51 |
|
#include "brains/SimCreator.hpp" |
| 52 |
|
#include "brains/SimSnapshotManager.hpp" |
| 55 |
|
#include "UseTheForce/ForceFieldFactory.hpp" |
| 56 |
|
#include "utils/simError.h" |
| 57 |
|
#include "utils/StringUtils.hpp" |
| 58 |
< |
#include "math/MersenneTwister.hpp" |
| 58 |
> |
#include "math/SeqRandNumGen.hpp" |
| 59 |
|
#ifdef IS_MPI |
| 60 |
|
#include "io/mpiBASS.h" |
| 61 |
+ |
#include "math/ParallelRandNumGen.hpp" |
| 62 |
|
#endif |
| 63 |
|
|
| 64 |
|
namespace oopse { |
| 253 |
|
|
| 254 |
|
int seedValue; |
| 255 |
|
Globals * simParams = info->getSimParams(); |
| 256 |
< |
MTRand* myRandom; //divide labor does not need Parallel random number generator |
| 256 |
> |
SeqRandNumGen* myRandom; //divide labor does not need Parallel random number generator |
| 257 |
|
if (simParams->haveSeed()) { |
| 258 |
|
seedValue = simParams->getSeed(); |
| 259 |
< |
myRandom = new MTRand(seedValue); |
| 259 |
> |
myRandom = new SeqRandNumGen(seedValue); |
| 260 |
|
}else { |
| 261 |
< |
myRandom = new MTRand(); |
| 261 |
> |
myRandom = new SeqRandNumGen(); |
| 262 |
|
} |
| 263 |
|
|
| 264 |
|
|
