57 |
|
#include "UseTheForce/ForceFieldFactory.hpp" |
58 |
|
#include "utils/simError.h" |
59 |
|
#include "utils/StringUtils.hpp" |
60 |
< |
|
60 |
> |
#include "math/SeqRandNumGen.hpp" |
61 |
|
#ifdef IS_MPI |
62 |
|
#include "io/mpiBASS.h" |
63 |
|
#include "math/ParallelRandNumGen.hpp" |
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 |
|
|