--- trunk/src/math/ParallelRandNumGen.cpp 2009/11/25 20:02:06 1390 +++ trunk/src/math/ParallelRandNumGen.cpp 2012/08/31 21:16:10 1793 @@ -36,7 +36,8 @@ * [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). * [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). * [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). - * [4] Vardeman & Gezelter, in progress (2009). + * [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). + * [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). */ #include "math/ParallelRandNumGen.hpp" @@ -50,10 +51,10 @@ namespace OpenMD { ParallelRandNumGen::ParallelRandNumGen(const uint32& oneSeed) { - const int masterNode = 0; unsigned long seed = oneSeed; #ifdef IS_MPI + const int masterNode = 0; MPI_Bcast(&seed, 1, MPI_UNSIGNED_LONG, masterNode, MPI_COMM_WORLD); #endif @@ -85,9 +86,9 @@ namespace OpenMD { ParallelRandNumGen::ParallelRandNumGen() { std::vector bigSeed; - const int masterNode = 0; int nProcessors; #ifdef IS_MPI + const int masterNode = 0; MPI_Comm_size( MPI_COMM_WORLD, &nProcessors); MPI_Comm_rank( MPI_COMM_WORLD, &myRank_); #else @@ -103,9 +104,9 @@ namespace OpenMD { void ParallelRandNumGen::seed( const uint32 oneSeed ) { - const int masterNode = 0; unsigned long seed = oneSeed; #ifdef IS_MPI + const int masterNode = 0; MPI_Bcast(&seed, 1, MPI_UNSIGNED_LONG, masterNode, MPI_COMM_WORLD); #endif if (seed != oneSeed) { @@ -125,8 +126,9 @@ namespace OpenMD { std::vector bigSeed; int size; - const int masterNode = 0; + #ifdef IS_MPI + const int masterNode = 0; if (worldRank == masterNode) { #endif