ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

revise automatic seed generation

Browse Source
This commit is contained in:
Axel Kohlmeyer
2021-10-08 11:34:21 -04:00
parent 4b86dbd200
commit 392ebf7db7
2 changed files with 10 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
src/DPD-MESO/pair_mdpd.cpp
View File

@ -217,9 +217,12 @@ void PairMDPD::settings(int narg, char **arg)
seed = utils::inumeric(FLERR,arg[2],false,lmp);
// initialize Marsaglia RNG with processor-unique seed
// create a positive seed based on the system clock, if requested.
if (seed <= 0)
seed = (int) (platform::walltime() * 1073741824.0);
if (seed <= 0) {
constexpr double LARGE_NUM = 2<<30;
seed = int(fmod(platform::walltime() * LARGE_NUM, LARGE_NUM)) + 1;
}
delete random;
random = new RanMars(lmp,(seed + comm->me) % 900000000);

7
src/DPD-MESO/pair_tdpd.cpp
View File

@ -239,9 +239,12 @@ void PairTDPD::settings(int narg, char **arg)
seed = utils::inumeric(FLERR,arg[2],false,lmp);
// initialize Marsaglia RNG with processor-unique seed
// create a positive seed based on the system clock, if requested.
if (seed <= 0)
seed = (int) (platform::walltime() * 1073741824.0);
if (seed <= 0) {
constexpr double LARGE_NUM = 2<<30;
seed = int(fmod(platform::walltime() * LARGE_NUM, LARGE_NUM)) + 1;
}
delete random;
random = new RanMars(lmp,(seed + comm->me) % 900000000);