From 03d6e6585adb966c419c3da9cd53f8eaa8b828b1 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Wed, 4 Sep 2013 12:34:08 +0100
Subject: [PATCH] BUG: autoLayerDriver: initialisation of face merge candidates
---
.../autoHexMesh/autoHexMeshDriver/autoLayerDriver.C | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C
index 1d9ecc0bc74..cd0ebafc3a4 100644
--- a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C
+++ b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C
@@ -2450,17 +2450,17 @@ void Foam::autoLayerDriver::mergePatchFacesUndo
const dictionary& motionDict
)
{
- scalar minCos =
- Foam::cos(degToRad(layerParams.featureAngle()));
+ // Clip to 45 degrees
+ scalar planarAngle = min(45.0, layerParams.featureAngle());
+ scalar minCos = Foam::cos(degToRad(planarAngle));
- scalar concaveCos =
- Foam::cos(degToRad(layerParams.concaveAngle()));
+ scalar concaveCos = Foam::cos(degToRad(layerParams.concaveAngle()));
Info<< nl
<< "Merging all faces of a cell" << nl
<< "---------------------------" << nl
<< " - which are on the same patch" << nl
- << " - which make an angle < " << layerParams.featureAngle()
+ << " - which make an angle < " << planarAngle
<< " degrees"
<< nl
<< " (cos:" << minCos << ')' << nl
@@ -2478,7 +2478,7 @@ void Foam::autoLayerDriver::mergePatchFacesUndo
concaveCos,
meshRefiner_.meshedPatches(),
motionDict,
- labelList(mesh.nFaces() -1)
+ labelList(mesh.nFaces(), -1)
);
nChanged += meshRefiner_.mergeEdgesUndo(minCos, motionDict);
--
GitLab