From 6e32ca010ac022c2e1f4f2a1beae3a0455ea99e6 Mon Sep 17 00:00:00 2001 From: Will Bainbridge Date: Mon, 27 Jan 2020 17:06:39 +0000 Subject: [PATCH] globalIndexAndTransform: Correction to handle rotation around a non-zero centre --- .../globalIndexAndTransform/globalIndexAndTransform.C | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/src/OpenFOAM/primitives/globalIndexAndTransform/globalIndexAndTransform.C b/src/OpenFOAM/primitives/globalIndexAndTransform/globalIndexAndTransform.C index c33e27ae50..8fd56044c8 100644 --- a/src/OpenFOAM/primitives/globalIndexAndTransform/globalIndexAndTransform.C +++ b/src/OpenFOAM/primitives/globalIndexAndTransform/globalIndexAndTransform.C @@ -96,8 +96,10 @@ Foam::label Foam::globalIndexAndTransform::matchTransform { // Test the inverse transform differences too + const transformer testInvTransform = inv(testTransform); + vectorDiff = - mag(refTransform.t() + testTransform.t()) + mag(refTransform.t() - testInvTransform.t()) /(maxVectorMag + vSmall) /tolerance; @@ -106,7 +108,7 @@ Foam::label Foam::globalIndexAndTransform::matchTransform if (refTransform.transforms() || testTransform.transforms()) { tensorDiff = - mag(refTransform.T() - testTransform.T().T()) + mag(refTransform.T() - testInvTransform.T()) /sqrt(3.0) /tolerance; }