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

Comparing trunk/src/integrators/DLM.cpp (file contents):
Revision 1879 by gezelter, Sun Jun 16 15:15:42 2013 UTC vs.
Revision 1880 by gezelter, Mon Jun 17 18:28:30 2013 UTC

# Line 108 | Line 108 | namespace OpenMD {
108      RotMat3x3d rot = RotMat3x3d::identity(); // initalize rot as a unit matrix
109  
110      // use a small angle aproximation for sin and cosine
111 <    
111 >
112      angleSqr = angle * angle;
113      angleSqrOver4 = angleSqr / 4.0;
114      top = 1.0 - angleSqrOver4;
# Line 116 | Line 116 | namespace OpenMD {
116  
117      cosAngle = top / bottom;
118      sinAngle = angle / bottom;
119 <    
119 >
120      // or don't use the small angle approximation:
121      //cosAngle = cos(angle);
122      //sinAngle = sin(angle);
123 –
123      rot(axes1, axes1) = cosAngle;
124      rot(axes2, axes2) = cosAngle;
125  

Diff Legend

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