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

Comparing trunk/src/brains/ForceManager.cpp (file contents):
Revision 2020 by gezelter, Mon Sep 22 19:18:35 2014 UTC vs.
Revision 2031 by jmichalk, Fri Oct 31 18:40:40 2014 UTC

# Line 843 | Line 843 | namespace OpenMD {
843                  } else {
844                    interactionMan_->doPair(idat);
845                    fDecomp_->unpackInteractionData(idat, atom1, atom2);
846 +                  /*
847 +                  if(atom1 == 971){
848 +                    cerr << "972: " << *(idat.dVdFQ1) << "\n";
849 +                  }
850 +                  if(atom2 == 971){
851 +                    cerr << "972: " << *(idat.dVdFQ2) << "\n";
852 +                  }
853 +                  if(atom1 == 766){
854 +                    cerr << "767: " << *(idat.dVdFQ1) << "\n";
855 +                  }
856 +                  if(atom2 == 766){
857 +                    cerr << "767: " << *(idat.dVdFQ2) << "\n";
858 +                  }
859 +                  */
860                    vij += vpair;
861                    fij += f1;
862                    stressTensor -= outProduct( *(idat.d), f1);
# Line 945 | Line 959 | namespace OpenMD {
959      if (info_->requiresSelfCorrection()) {
960        for (unsigned int atom1 = 0; atom1 < info_->getNAtoms(); atom1++) {
961          fDecomp_->fillSelfData(sdat, atom1);
962 +        /*      
963 +        if(atom1 == 971 || atom1 == 766 || atom1 == 898){
964 +          cout << "atom id: " << atom1 << "\tsdat.dVdFQ: " << *(sdat.dVdFQ) << "\n";
965 +        }
966 +        */
967          interactionMan_->doSelfCorrection(sdat);
968 +        /*
969 +        if(atom1 == 971 || atom1 == 766 || atom1 == 898){
970 +          cout << "\tsdat.dVdFQ: " << *(sdat.dVdFQ) << "\n";
971 +        }
972 +        */
973        }
974      }
975  

Diff Legend

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