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

Comparing branches/development/src/brains/ForceManager.cpp (file contents):
Revision 1590 by gezelter, Mon Jul 11 01:39:49 2011 UTC vs.
Revision 1601 by gezelter, Thu Aug 4 20:04:35 2011 UTC

# Line 386 | Line 386 | namespace OpenMD {
386      }
387      
388      // Zero out the stress tensor
389 <    tau = Mat3x3d(0.0);
389 >    tau *= 0.0;
390      
391    }
392    
# Line 618 | Line 618 | namespace OpenMD {
618              for (vector<int>::iterator jb = atomListColumn.begin();
619                   jb != atomListColumn.end(); ++jb) {              
620                atom2 = (*jb);
621 <            
621 >
622                if (!fDecomp_->skipAtomPair(atom1, atom2)) {
623                  vpair = 0.0;
624                  workPot = 0.0;
# Line 640 | Line 640 | namespace OpenMD {
640                    idat.d = &d;
641                    idat.r2 = &r2;
642                  }
643 <                
643 >              
644                  r = sqrt( *(idat.r2) );
645                  idat.rij = &r;
646                
# Line 649 | Line 649 | namespace OpenMD {
649                  } else {
650                    interactionMan_->doPair(idat);
651                    fDecomp_->unpackInteractionData(idat, atom1, atom2);
652 +
653                    vij += vpair;
654                    fij += f1;
655                    tau -= outProduct( *(idat.d), f1);

Diff Legend

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