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 1592 by gezelter, Tue Jul 12 20:33:14 2011 UTC vs.
Revision 1601 by gezelter, Thu Aug 4 20:04:35 2011 UTC

# 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;
654 –                  cerr << "ats = " << atom1 << " " << atom2 << "d = " << *(idat.d) << " f1 = " << f1 << "\n";
655                    tau -= outProduct( *(idat.d), f1);
656                  }
657                }

Diff Legend

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