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 1808 by gezelter, Mon Oct 22 20:42:10 2012 UTC vs.
Revision 1847 by gezelter, Thu Jan 31 17:57:07 2013 UTC

# Line 745 | Line 745 | namespace OpenMD {
745            if (iLoop == PAIR_LOOP) {
746              vij = 0.0;
747              fij = V3Zero;
748 +            eField1 = V3Zero;
749 +            eField2 = V3Zero;
750            }
751            
752            in_switching_region = switcher_->getSwitch(rgrpsq, sw, dswdr,
# Line 819 | Line 821 | namespace OpenMD {
821                fij += fg;
822  
823                if (atomListRow.size() == 1 && atomListColumn.size() == 1) {
824 <                stressTensor -= outProduct( *(idat.d), fg);
825 <                if (doHeatFlux_)
826 <                  fDecomp_->addToHeatFlux(*(idat.d) * dot(fg, vel2));
827 <                
824 >                if (!fDecomp_->skipAtomPair(atomListRow[0],
825 >                                            atomListColumn[0],
826 >                                            cg1, cg2)) {
827 >                  stressTensor -= outProduct( *(idat.d), fg);
828 >                  if (doHeatFlux_)
829 >                    fDecomp_->addToHeatFlux(*(idat.d) * dot(fg, vel2));
830 >                }                
831                }
832            
833                for (ia = atomListRow.begin();

Diff Legend

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