exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Bugfix in faceZone reconstruction procedure

Browse source 
Forgot to take into account the "turning face" index and decrement by one
This commit is contained in:
Vuko Vukcevic 2019-06-24 09:51:40 +02:00
parent 49769f0397
commit 15af49e74a
1 changed files with 13 additions and 8 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

21
src/decompositionMethods/decomposeReconstruct/reconstructTools/finiteVolume/processorMeshesRebuild.C
View file

@ -1915,18 +1915,23 @@ Foam::processorMeshesReconstructor::reconstructMesh(const Time& db)
// Get the face index in the new, reconstructed mesh
const label& newFaceI = curFaceProcAddr[oldFaceI];
// Check if the face is mapped (if the index is >= 0)
if (newFaceI > -1)
// Check if the face is mapped.
// Note:
// 1. Need to decrement by 1 because of the face turning
// 2. No need to handle negative new indices coming from
// slave processor because we'd end up with
// duplicate entries (two faces on two processors
// merged into a single one)
if (newFaceI > 0)
{
// This is a face that's been correctly
// mapped, insert the face in the new zone
zoneReconFaces[nFaces] = newFaceI;
zoneReconFaces[nFaces] = newFaceI - 1;
// Also store the flip map of the face. Note: I'm
// pretty sure that we don't need to check whether
// the flip map has been preserved because the
// combined faces are inserted from master side
// always.
// Also store the flip map of the face. We don't
// need to check whether the flip map has been
// preserved because we only get the combined faces
// that are inserted from master side.
zoneReconFaceFlips[nFaces] = zoneFlipMap[i];
// Increment the number of faces for this zone