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

Comparing branches/development/src/restraints/ThermoIntegrationForceManager.cpp (file contents):
Revision 1723 by gezelter, Thu May 24 20:59:54 2012 UTC vs.
Revision 1760 by gezelter, Thu Jun 21 19:26:46 2012 UTC

# Line 125 | Line 125 | namespace OpenMD {
125        }
126      }
127      
128 <    // set vraw to be the unmodulated potential
129 <    lrPot_ = curSnapshot->statData[Stats::LONG_RANGE_POTENTIAL];
130 <    curSnapshot->statData[Stats::VRAW] = lrPot_;
128 >    // set rawPotential to be the unmodulated potential
129 >    lrPot_ = curSnapshot->getLongRangePotential();
130 >    curSnapshot->setRawPotential(lrPot_);
131      
132      // modulate the potential and update the snapshot
133      lrPot_ *= factor_;
134 <    curSnapshot->statData[Stats::LONG_RANGE_POTENTIAL] = lrPot_;
134 >    curSnapshot->setLongRangePotential(lrPot_);
135      
136      // scale the pressure tensor
137      tempTau = curSnapshot->getStressTensor();
# Line 144 | Line 144 | namespace OpenMD {
144      
145      if (simParam->getUseRestraints()) {
146        // do restraints from RestraintForceManager:
147 –      //restPot_local = doRestraints(1.0 - factor_);    
147        restPot_local = doRestraints(1.0 - factor_);      
148        vHarm_local = getUnscaledPotential();
149      }
# Line 162 | Line 161 | namespace OpenMD {
161   #endif
162  
163      // give the final values to stats
164 <    curSnapshot->statData[Stats::LONG_RANGE_POTENTIAL] = lrPot_;
165 <    curSnapshot->statData[Stats::VHARM] = vHarm_;
164 >    curSnapshot->setLongRangePotential(lrPot_);
165 >    curSnapshot->setRestraintPotential(vHarm_);
166    }  
167   }

Diff Legend

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