From 14d9cf32bed13af1f6d798e793257b56d2d83399 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Wed, 1 Jun 2011 18:08:05 +0100
Subject: [PATCH] BUG: decompositionMethod: agglomeration across coupled
patches
---
.../decompositionMethod/decompositionMethod.C | 10 +++++++---
1 file changed, 7 insertions(+), 3 deletions(-)
diff --git a/src/parallel/decompose/decompositionMethods/decompositionMethod/decompositionMethod.C b/src/parallel/decompose/decompositionMethods/decompositionMethod/decompositionMethod.C
index 1cfb22887b7..44d3cc027b7 100644
--- a/src/parallel/decompose/decompositionMethods/decompositionMethod/decompositionMethod.C
+++ b/src/parallel/decompose/decompositionMethods/decompositionMethod/decompositionMethod.C
@@ -448,8 +448,8 @@ void Foam::decompositionMethod::calcCellCells
// Count number of faces (internal + coupled)
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- // Number of faces per cell
- labelList nFacesPerCell(mesh.nCells(), 0);
+ // Number of faces per coarse cell
+ labelList nFacesPerCell(nCoarse, 0);
for (label faceI = 0; faceI < mesh.nInternalFaces(); faceI++)
{
@@ -481,7 +481,11 @@ void Foam::decompositionMethod::calcCellCells
{
label own = agglom[faceOwner[faceI]];
label globalNei = globalNeighbour[bFaceI];
- if (cellPair.insert(labelPair(own, globalNei)))
+ if
+ (
+ globalAgglom.toGlobal(own) != globalNei
+ && cellPair.insert(labelPair(own, globalNei))
+ )
{
nFacesPerCell[own]++;
}
--
GitLab