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

Comparing trunk/src/applications/staticProps/P2OrderParameter.cpp (file contents):
Revision 544 by tim, Fri May 27 04:41:34 2005 UTC vs.
Revision 545 by tim, Fri May 27 21:11:24 2005 UTC

# Line 125 | Line 125 | void P2OrderParameter::process() {
125        Mat3x3d orderTensor(0.0);
126        for (std::vector<std::pair<StuntDouble*, StuntDouble*> >::iterator j = sdPairs_.begin(); j != sdPairs_.end(); ++j) {
127            Vector3d vec = j->first->getPos() - j->second->getPos();
128 +          currentSnapshot_->wrapVector(vec);
129            vec.normalize();
130            orderTensor +=outProduct(vec, vec);
131        }
# Line 155 | Line 156 | void P2OrderParameter::process() {
156        double angle = 0.0;
157        for (std::vector<std::pair<StuntDouble*, StuntDouble*> >::iterator j = sdPairs_.begin(); j != sdPairs_.end(); ++j) {
158            Vector3d vec = j->first->getPos() - j->second->getPos();
159 +          currentSnapshot_->wrapVector(vec);
160            vec.normalize();
161  
162            angle += acos(dot(vec, director)) ;

Diff Legend

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