From 957f8d9d70b4c8644da44890f17029a36cf18804 Mon Sep 17 00:00:00 2001 From: mattijs Date: Wed, 15 Dec 2010 15:35:50 +0000 Subject: [PATCH] COMP: transform : split off operations --- src/OpenFOAM/primitives/transform/transform.H | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/src/OpenFOAM/primitives/transform/transform.H b/src/OpenFOAM/primitives/transform/transform.H index 70e6821612..c880e05888 100644 --- a/src/OpenFOAM/primitives/transform/transform.H +++ b/src/OpenFOAM/primitives/transform/transform.H @@ -48,9 +48,11 @@ inline tensor rotationTensor const vector& n2 ) { + const scalar s = n1 & n2; + const vector n3 = n1 ^ n2; return - (n1 & n2)*I - + (1 - (n1 & n2))*sqr(n1 ^ n2)/(magSqr(n1 ^ n2) + VSMALL) + s*I + + (1 - s)*sqr(n3)/(magSqr(n3) + VSMALL) + (n2*n1 - n1*n2); }