| 45 |
|
#include "utils/simError.h" |
| 46 |
|
namespace OpenMD { |
| 47 |
|
|
| 48 |
< |
Shake::Shake(SimInfo* info) : info_(info), maxConsIteration_(10), consTolerance_(1.0e-6), doShake_(false) { |
| 48 |
> |
Shake::Shake(SimInfo* info) : info_(info), maxConsIteration_(10), |
| 49 |
> |
consTolerance_(1.0e-6), doShake_(false), |
| 50 |
> |
currConstraintTime_(0.0) { |
| 51 |
|
|
| 52 |
|
if (info_->getNGlobalConstraints() > 0) |
| 53 |
|
doShake_ = true; |
| 61 |
|
constraintTime_ = simParams->getStatusTime(); |
| 62 |
|
} |
| 63 |
|
|
| 64 |
< |
constraintOutputFile_ = getPrefix(info_->getFinalConfigFileName()) + ".constraintForces"; |
| 64 |
> |
constraintOutputFile_ = getPrefix(info_->getFinalConfigFileName()) + |
| 65 |
> |
".constraintForces"; |
| 66 |
|
|
| 67 |
|
// create ConstraintWriter |
| 68 |
< |
constraintWriter_ = new ConstraintWriter(info_, constraintOutputFile_.c_str()); |
| 68 |
> |
constraintWriter_ = new ConstraintWriter(info_, |
| 69 |
> |
constraintOutputFile_.c_str()); |
| 70 |
|
|
| 71 |
|
if (!constraintWriter_){ |
| 72 |
|
sprintf(painCave.errMsg, "Failed to create ConstraintWriter\n"); |