From 9afbd6805901a248d424503b1240689b94373e51 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Fri, 3 Dec 2010 12:03:37 +0000
Subject: [PATCH] ENH: cyclicPolyPatch : disable warning about faces not
perpendicular to rotation axis
---
.../constraint/cyclic/cyclicPolyPatch.C | 50 +++++++++----------
1 file changed, 25 insertions(+), 25 deletions(-)
diff --git a/src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C b/src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C
index bd06a2755b8..e7d1991fb95 100644
--- a/src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C
+++ b/src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C
@@ -194,8 +194,9 @@ void Foam::cyclicPolyPatch::calcTransforms
vectorField half0Normals(half0Areas.size());
vectorField half1Normals(half1Areas.size());
- scalar maxCos = -GREAT;
- label maxFacei = -1;
+ //- Additional warning about faces non-aligned with rotation axis
+ //scalar maxCos = -GREAT;
+ //label maxFacei = -1;
forAll(half0, facei)
{
@@ -237,32 +238,31 @@ void Foam::cyclicPolyPatch::calcTransforms
half0Normals[facei] = half0Areas[facei] / magSf;
half1Normals[facei] = half1Areas[facei] / nbrMagSf;
- if (transform_ == ROTATIONAL)
- {
- scalar cos = mag(half0Normals[facei] & rotationAxis_);
- if (cos > maxCos)
- {
- maxCos = cos;
- maxFacei = facei;
- }
- }
+ //if (transform_ == ROTATIONAL)
+ //{
+ // scalar cos = mag(half0Normals[facei] & rotationAxis_);
+ // if (cos > maxCos)
+ // {
+ // maxCos = cos;
+ // maxFacei = facei;
+ // }
+ //}
}
}
- if (maxCos > sqrt(SMALL))
- {
- WarningIn
- (
- "cyclicPolyPatch::calcTransforms()"
- ) << "on patch " << name()
- << " face:" << maxFacei << " fc:" << half0Ctrs[maxFacei]
- << " is not perpendicular to the rotationAxis." << endl
- << "This might cause conservation problems"
- << " or problems with geometry or topology changes." << endl
- << "rotation axis : " << rotationAxis_ << endl
- << "face normal : " << half0Normals[maxFacei] << endl
- << "cosine of angle : " << maxCos << endl;
- }
+ //if (maxCos > sqrt(SMALL))
+ //{
+ // WarningIn
+ // (
+ // "cyclicPolyPatch::calcTransforms()"
+ // ) << "on patch " << name()
+ // << " face:" << maxFacei << " fc:" << half0Ctrs[maxFacei]
+ // << " is not perpendicular to the rotationAxis." << endl
+ // << "This will cause problems with topology changes." << endl
+ // << "rotation axis : " << rotationAxis_ << endl
+ // << "face normal : " << half0Normals[maxFacei] << endl
+ // << "cosine of angle : " << maxCos << endl;
+ //}
// Calculate transformation tensors
calcTransformTensors
--
GitLab