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

Comparing trunk/src/restraints/ThermoIntegrationForceManager.cpp (file contents):
Revision 417 by chrisfen, Thu Mar 10 15:10:24 2005 UTC vs.
Revision 419 by chrisfen, Thu Mar 10 16:15:13 2005 UTC

# Line 90 | Line 90 | namespace oopse {
90      
91      // build the scaling factor used to modulate the forces and torques
92      factor_ = pow(tIntLambda_, tIntK_);
93 +
94 +    printf("%f is the factor\n",factor_);
95      
96    }
97    
# Line 106 | Line 108 | namespace oopse {
108      Vector3d frc;
109      Vector3d trq;
110      
109 –    curSnapshot = info_->getSnapshotManager()->getCurrentSnapshot();
110 –    
111      // perform the standard calcForces first
112      ForceManager::calcForces(needPotential, needStress);
113      
114 +    curSnapshot = info_->getSnapshotManager()->getCurrentSnapshot();
115 +
116      // now scale forces and torques of all the integrableObjects
117        
118      for (mol = info_->beginMolecule(mi); mol != NULL;
# Line 130 | Line 132 | namespace oopse {
132        }
133        
134        // set vraw to be the unmodulated potential
135 <      lrPot_ = curSnapshot->statData[Stats::POTENTIAL_ENERGY];
135 >      lrPot_ = curSnapshot->statData[Stats::LONG_RANGE_POTENTIAL];
136        curSnapshot->statData[Stats::VRAW] = lrPot_;
137        
138        // modulate the potential and update the snapshot
139        lrPot_ *= factor_;
140 <      curSnapshot->statData[Stats::POTENTIAL_ENERGY] = lrPot_;
140 >      curSnapshot->statData[Stats::LONG_RANGE_POTENTIAL] = lrPot_;
141      }
142      
143      // do crystal restraint forces for thermodynamic integration
144      if (simParam->getUseSolidThermInt()) {
145        
146        lrPot_ += restraint_->Calc_Restraint_Forces();
147 <      curSnapshot->statData[Stats::POTENTIAL_ENERGY] = lrPot_;
147 >      curSnapshot->statData[Stats::LONG_RANGE_POTENTIAL] = lrPot_;
148        
149        vHarm_ = restraint_->getVharm();
150        curSnapshot->statData[Stats::VHARM] = vHarm_;

Diff Legend

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