123 |
|
|
124 |
|
initialize(); |
125 |
|
|
126 |
< |
while (snap->getTime() < runTime) { |
126 |
> |
while (snap->getTime() < runTime) { |
127 |
|
preStep(); |
128 |
|
integrateStep(); |
129 |
|
postStep(); |
143 |
|
|
144 |
|
void VelocityVerletIntegrator::postStep() { |
145 |
|
|
146 |
– |
//save snapshot |
147 |
– |
info_->getSnapshotManager()->advance(); |
148 |
– |
|
149 |
– |
//increase time |
150 |
– |
snap->increaseTime(dt); |
146 |
|
|
147 |
|
if (needVelocityScaling) { |
148 |
|
if (snap->getTime() >= currThermal) { |
166 |
|
|
167 |
|
if (snap->getTime() >= currStatus) { |
168 |
|
//save statistics, before writeStat, we must save statistics |
174 |
– |
stats->collectStats(); |
169 |
|
saveConservedQuantity(); |
170 |
+ |
stats->collectStats(); |
171 |
|
|
172 |
|
if (simParams->getRNEMDParameters()->getUseRNEMD()) { |
173 |
|
rnemd_->writeOutputFile(); |
187 |
|
resetIntegrator(); |
188 |
|
currReset += resetTime; |
189 |
|
} |
190 |
+ |
//save snapshot |
191 |
+ |
info_->getSnapshotManager()->advance(); |
192 |
+ |
|
193 |
+ |
//increase time |
194 |
+ |
snap->increaseTime(dt); |
195 |
+ |
|
196 |
|
} |
197 |
|
|
198 |
|
|