From 806bdec1786946d42cfdf9e5e5af189d50dc3db0 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Tue, 13 Mar 2012 13:03:31 +0000
Subject: [PATCH] BUG: axisRotationMotion: fix angle specification

---
 .../axisRotationMotion/axisRotationMotion.C                | 7 +------
 1 file changed, 1 insertion(+), 6 deletions(-)

diff --git a/src/dynamicFvMesh/solidBodyMotionFvMesh/solidBodyMotionFunctions/axisRotationMotion/axisRotationMotion.C b/src/dynamicFvMesh/solidBodyMotionFvMesh/solidBodyMotionFunctions/axisRotationMotion/axisRotationMotion.C
index 456441026f..2de3fab7e0 100644
--- a/src/dynamicFvMesh/solidBodyMotionFvMesh/solidBodyMotionFunctions/axisRotationMotion/axisRotationMotion.C
+++ b/src/dynamicFvMesh/solidBodyMotionFvMesh/solidBodyMotionFunctions/axisRotationMotion/axisRotationMotion.C
@@ -73,8 +73,6 @@ Foam::solidBodyMotionFunctions::axisRotationMotion::transformation() const
 {
     scalar t = time_.value();
 
-    // Motion around a centre of gravity
-
     // Rotation around centre of gravity (in radians)
     vector omega
     (
@@ -84,10 +82,7 @@ Foam::solidBodyMotionFunctions::axisRotationMotion::transformation() const
     );
 
     scalar magOmega = mag(omega);
-
-    scalar cosHalfTheta = cos(0.5*magOmega);
-
-    quaternion R(cosHalfTheta, omega/magOmega);
+    quaternion R(omega/magOmega, magOmega);
     septernion TR(septernion(CofG_)*R*septernion(-CofG_));
 
     Info<< "solidBodyMotionFunctions::axisRotationMotion::transformation(): "