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 272 by gezelter, Wed Jan 12 22:41:40 2005 UTC vs.
Revision 273 by tim, Tue Jan 25 17:45:23 2005 UTC

# Line 109 | Line 109 | void DLM::rotateStep(int axes1, int axes2, double angl
109  
110      // use a small angle aproximation for sin and cosine
111  
112 <    angleSqr = angle * angle;
113 <    angleSqrOver4 = angleSqr / 4.0;
114 <    top = 1.0 - angleSqrOver4;
115 <    bottom = 1.0 + angleSqrOver4;
112 >    //angleSqr = angle * angle;
113 >    //angleSqrOver4 = angleSqr / 4.0;
114 >    //top = 1.0 - angleSqrOver4;
115 >    //bottom = 1.0 + angleSqrOver4;
116  
117 <    cosAngle = top / bottom;
118 <    sinAngle = angle / bottom;
119 <
117 >    //cosAngle = top / bottom;
118 >    //sinAngle = angle / bottom;
119 >    cosAngle = cos(angle);
120 >    sinAngle = sin(angle);
121      rot(axes1, axes1) = cosAngle;
122      rot(axes2, axes2) = cosAngle;
123  

Diff Legend

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