| 834 |
|
} |
| 835 |
|
|
| 836 |
|
//check whether sample time, status time, thermal time and reset time are divisble by dt |
| 837 |
< |
if (!isDivisible(globals->getSampleTime(), globals->getDt())){ |
| 837 |
> |
if (globals->haveSampleTime() && !isDivisible(globals->getSampleTime(), globals->getDt())){ |
| 838 |
|
sprintf(painCave.errMsg, |
| 839 |
|
"Sample time is not divisible by dt.\n" |
| 840 |
|
"\tThis will result in samples that are not uniformly\n" |
| 844 |
|
simError(); |
| 845 |
|
} |
| 846 |
|
|
| 847 |
< |
if (globals->haveStatusTime() && !isDivisible(globals->getSampleTime(), globals->getDt())){ |
| 847 |
> |
if (globals->haveStatusTime() && !isDivisible(globals->getStatusTime(), globals->getDt())){ |
| 848 |
|
sprintf(painCave.errMsg, |
| 849 |
|
"Status time is not divisible by dt.\n" |
| 850 |
|
"\tThis will result in status reports that are not uniformly\n" |
| 880 |
|
if (globals->haveSampleTime()){ |
| 881 |
|
info[i].sampleTime = globals->getSampleTime(); |
| 882 |
|
info[i].statusTime = info[i].sampleTime; |
| 883 |
– |
info[i].thermalTime = info[i].sampleTime; |
| 883 |
|
} |
| 884 |
|
else{ |
| 885 |
|
info[i].sampleTime = globals->getRunTime(); |
| 886 |
|
info[i].statusTime = info[i].sampleTime; |
| 888 |
– |
info[i].thermalTime = info[i].sampleTime; |
| 887 |
|
} |
| 888 |
|
|
| 889 |
|
if (globals->haveStatusTime()){ |
| 892 |
|
|
| 893 |
|
if (globals->haveThermalTime()){ |
| 894 |
|
info[i].thermalTime = globals->getThermalTime(); |
| 895 |
+ |
} else { |
| 896 |
+ |
info[i].thermalTime = globals->getRunTime(); |
| 897 |
|
} |
| 898 |
|
|
| 899 |
|
info[i].resetIntegrator = 0; |
| 1365 |
|
int localMol, allMol; |
| 1366 |
|
int local_atoms, local_bonds, local_bends, local_torsions, local_SRI; |
| 1367 |
|
int local_rigid; |
| 1368 |
– |
vector<int> globalAtomIndex; |
| 1368 |
|
vector<int> globalMolIndex; |
| 1369 |
|
|
| 1370 |
|
mpiSim = new mpiSimulation(info); |
| 1371 |
|
|
| 1372 |
|
mpiSim->divideLabor(); |
| 1373 |
|
globalAtomIndex = mpiSim->getGlobalAtomIndex(); |
| 1374 |
< |
globalMolIndex = mpiSim->getGlobalMolIndex(); |
| 1374 |
> |
//globalMolIndex = mpiSim->getGlobalMolIndex(); |
| 1375 |
|
|
| 1376 |
|
// set up the local variables |
| 1377 |
|
|