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

Comparing branches/development/src/nonbonded/MAW.cpp (file contents):
Revision 1582 by gezelter, Tue Jun 14 20:41:44 2011 UTC vs.
Revision 1583 by gezelter, Thu Jun 16 22:00:08 2011 UTC

# Line 50 | Line 50 | namespace OpenMD {
50  
51   namespace OpenMD {
52  
53 <  MAW::MAW() : name_("MAW"), initialized_(false), forceField_(NULL),
54 <                   shiftedPot_(false), shiftedFrc_(false) {}
53 >  MAW::MAW() : name_("MAW"), initialized_(false), forceField_(NULL) {}
54    
55    void MAW::initialize() {    
56  
# Line 172 | Line 171 | namespace OpenMD {
171        myPot  = D_e * (expfnc2  - 2.0 * expfnc);
172        myDeriv   = 2.0 * D_e * beta * (expfnc - expfnc2);
173        
174 <      if (MAW::shiftedPot_ || MAW::shiftedFrc_) {
174 >      if (idat.shiftedPot || idat.shiftedForce) {
175          exptC     = -beta*( *(idat.rcut)  - R_e);
176          expfncC   = exp(exptC);
177          expfnc2C  = expfncC*expfncC;
178        }
179        
180 <      if (MAW::shiftedPot_) {
180 >      if (idat.shiftedPot) {
181          myPotC = D_e * (expfnc2C - 2.0 * expfncC);
182          myDerivC = 0.0;
183 <      } else if (MAW::shiftedFrc_) {
183 >      } else if (idat.shiftedForce) {
184          myPotC = D_e * (expfnc2C - 2.0 * expfncC);
185          myDerivC  = 2.0 * D_e * beta * (expfnc2C - expfnc2C);
186          myPotC += myDerivC * ( *(idat.rij)  -  *(idat.rcut) );

Diff Legend

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