ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/trunk/src/constraints/Rattle.cpp
(Generate patch)

Comparing trunk/src/constraints/Rattle.cpp (file contents):
Revision 1982 by gezelter, Tue Apr 15 12:12:23 2014 UTC vs.
Revision 1983 by gezelter, Tue Apr 15 20:36:19 2014 UTC

# Line 45 | Line 45 | namespace OpenMD {
45   #include "utils/simError.h"
46   namespace OpenMD {
47  
48 <  Rattle::Rattle(SimInfo* info) : info_(info), maxConsIteration_(10), consTolerance_(1.0e-6), doRattle_(false) {
48 >  Rattle::Rattle(SimInfo* info) : info_(info), maxConsIteration_(10),
49 >                                  consTolerance_(1.0e-6), doRattle_(false),
50 >                                  currConstraintTime_(0.0) {
51      
52      if (info_->getNGlobalConstraints() > 0)
53        doRattle_ = true;
# Line 68 | Line 70 | namespace OpenMD {
70        constraintTime_ = simParams->getStatusTime();
71      }
72  
73 <    constraintOutputFile_ = getPrefix(info_->getFinalConfigFileName()) + ".constraintForces";
73 >    constraintOutputFile_ = getPrefix(info_->getFinalConfigFileName()) +
74 >      ".constraintForces";
75  
76      // create ConstraintWriter  
77 <    constraintWriter_ = new ConstraintWriter(info_, constraintOutputFile_.c_str());  
77 >    constraintWriter_ = new ConstraintWriter(info_,
78 >                                             constraintOutputFile_.c_str());  
79  
80      if (!constraintWriter_){
81        sprintf(painCave.errMsg, "Failed to create ConstraintWriter\n");
# Line 102 | Line 106 | namespace OpenMD {
106            constraints.push_back(consPair);
107          }
108        }
105 –      
109        constraintWriter_->writeConstraintForces(constraints);
110        currConstraintTime_ += constraintTime_;
111      }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines