ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE-4/src/applications/staticProps/GofAngle2.cpp
(Generate patch)

Comparing trunk/OOPSE-4/src/applications/staticProps/GofAngle2.cpp (file contents):
Revision 2759 by tim, Wed May 17 21:51:42 2006 UTC vs.
Revision 3054 by gezelter, Wed Oct 18 21:58:48 2006 UTC

# Line 93 | Line 93 | namespace oopse {
93      Vector3d pos1 = sd1->getPos();
94      Vector3d pos2 = sd2->getPos();
95      Vector3d r12 = pos1 - pos2;
96 <    currentSnapshot_->wrapVector(r12);
96 >    if (usePeriodicBoundaryConditions_)
97 >      currentSnapshot_->wrapVector(r12);
98      Vector3d dipole1 = sd1->getElectroFrame().getColumn(2);
99      Vector3d dipole2 = sd2->getElectroFrame().getColumn(2);
100      
# Line 107 | Line 108 | namespace oopse {
108  
109      RealType halfBin = (nAngleBins_ - 1) * 0.5;
110      int angleBin1 = halfBin * (cosAngle1 + 1.0);
111 <    int angleBin2 = halfBin * (cosAngle1 + 1.0);
111 >    int angleBin2 = halfBin * (cosAngle2 + 1.0);
112  
113 <    ++histogram_[angleBin1][angleBin1];    
113 >    ++histogram_[angleBin1][angleBin2];    
114      ++npairs_;
115    }
116  

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines