From 0a8167203feb85f4c840fbd71bde7c7b64d79460 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 12 Mar 2012 10:32:58 +0000
Subject: [PATCH] BUG: decompositionMethod.C: multiple faces inbetween cells
---
.../decompositionMethod/decompositionMethod.C | 19 ++++++++++++++++---
1 file changed, 16 insertions(+), 3 deletions(-)
diff --git a/src/parallel/decompose/decompositionMethods/decompositionMethod/decompositionMethod.C b/src/parallel/decompose/decompositionMethods/decompositionMethod/decompositionMethod.C
index 573b6e52066..ab18fad6e27 100644
--- a/src/parallel/decompose/decompositionMethods/decompositionMethod/decompositionMethod.C
+++ b/src/parallel/decompose/decompositionMethods/decompositionMethod/decompositionMethod.C
@@ -307,12 +307,20 @@ void Foam::decompositionMethod::calcCellCells
// Done as postprocessing step since we now have cellCells.
label newIndex = 0;
labelHashSet nbrCells;
+
+
+ if (cellCells.size() == 0)
+ {
+ return;
+ }
+
+ label startIndex = cellCells.offsets()[0];
+
forAll(cellCells, cellI)
{
nbrCells.clear();
nbrCells.insert(globalAgglom.toGlobal(cellI));
- label startIndex = cellCells.offsets()[cellI];
label endIndex = cellCells.offsets()[cellI+1];
for (label i = startIndex; i < endIndex; i++)
@@ -322,16 +330,21 @@ void Foam::decompositionMethod::calcCellCells
cellCells.m()[newIndex++] = cellCells.m()[i];
}
}
+ startIndex = endIndex;
cellCells.offsets()[cellI+1] = newIndex;
}
cellCells.m().setSize(newIndex);
-
//forAll(cellCells, cellI)
//{
- // const labelUList cCells = cellCells[cellI];
+ // Pout<< "Original: Coarse cell " << cellI << endl;
+ // forAll(mesh.cellCells()[cellI], i)
+ // {
+ // Pout<< " nbr:" << mesh.cellCells()[cellI][i] << endl;
+ // }
// Pout<< "Compacted: Coarse cell " << cellI << endl;
+ // const labelUList cCells = cellCells[cellI];
// forAll(cCells, i)
// {
// Pout<< " nbr:" << cCells[i] << endl;
--
GitLab