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"); |