zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: polyTopoChange: efficient sorting

Browse Source
This commit is contained in:
mattijs
2011-10-28 11:31:23 +01:00
parent d666fda130
commit a0352b0ae0
1 changed files with 20 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

26
src/dynamicMesh/polyTopoChange/polyTopoChange/polyTopoChange.C
View File

@ -680,13 +680,17 @@ void Foam::polyTopoChange::getFaceOrder
// First unassigned face
label newFaceI = 0;
labelList nbr;
labelList order;
forAll(cellMap_, cellI)
{
label startOfCell = cellFaceOffsets[cellI];
label nFaces = cellFaceOffsets[cellI+1] - startOfCell;
// Neighbouring cells
SortableList<label> nbr(nFaces);
//SortableList<label> nbr(nFaces);
nbr.setSize(nFaces);
for (label i = 0; i < nFaces; i++)
{
@ -725,14 +729,24 @@ void Foam::polyTopoChange::getFaceOrder
}
}
nbr.sort();
//nbr.sort();
order.setSize(nFaces);
sortedOrder(nbr, order);
forAll(nbr, i)
//forAll(nbr, i)
//{
// if (nbr[i] != -1)
// {
// oldToNew[cellFaces[startOfCell + nbr.indices()[i]]] =
// newFaceI++;
// }
//}
forAll(order, i)
{
if (nbr[i] != -1)
label index = order[i];
if (nbr[index] != -1)
{
oldToNew[cellFaces[startOfCell + nbr.indices()[i]]] =
newFaceI++;
oldToNew[cellFaces[startOfCell + index]] = newFaceI++;
}
}
}