From 390843cfea2625c1a37c9a26b6f015c5442cbcbd Mon Sep 17 00:00:00 2001
From: henry <henry>
Date: Fri, 27 Mar 2009 09:08:54 +0000
Subject: [PATCH] Changed non-orthogonality check to handle both SP and DP.

---
 .../coordinateRotation/coordinateRotation.C            | 18 ++----------------
 1 file changed, 2 insertions(+), 16 deletions(-)

diff --git a/src/meshTools/coordinateSystems/coordinateRotation/coordinateRotation.C b/src/meshTools/coordinateSystems/coordinateRotation/coordinateRotation.C
index 649013f..30becd5 100644
--- a/src/meshTools/coordinateSystems/coordinateRotation/coordinateRotation.C
+++ b/src/meshTools/coordinateSystems/coordinateRotation/coordinateRotation.C
@@ -45,15 +45,13 @@ void Foam::coordinateRotation::calcTransform
     const axisOrder& order
 )
 {
-    const scalar orthogonalityError = 1.0e-8;
-
     vector a = axis1 / mag(axis1);
     vector b = axis2;
 
-    // Absorb minor nonorthogonality into axis2
+    // Absorb minor non-orthogonality into axis2
     b = b - (b & a)*a;
 
-    if (mag(b) < orthogonalityError)
+    if (mag(b) < SMALL)
     {
         FatalErrorIn
         (
@@ -64,18 +62,6 @@ void Foam::coordinateRotation::calcTransform
             << abort(FatalError);
     }
 
-    // this error check is probably no longer required
-    if (mag(a & b)/(mag(a)*mag(b)) >= orthogonalityError)
-    {
-        FatalErrorIn
-        (
-            "coordinateRotation::calcTransform()"
-        )
-            << "coordinate system nonorthogonality " << nl
-            << "mag(axis1 & axis2) = " << mag(a & b)
-            << abort(FatalError);
-    }
-
     b = b / mag(b);
     vector c = a ^ b;
 
-- 
2.11.4.GIT