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

BUG: correct the EnSight point value patch conversion (#2010)

Browse Source 
- patch point ordering was affected by the EnSight face sorting
  (tri, quad, nsided). This did not affect other point fields,
  since they use the internal field values.
This commit is contained in:
Mark Olesen
2021-03-10 18:58:15 +01:00
parent 9dc70ab856
commit 273741c052
1 changed files with 31 additions and 19 deletions
Download Patch File Download Diff File Expand all files Collapse all files

50
src/conversion/ensight/output/ensightOutputVolFieldTemplates.C
View File

@ -224,14 +224,14 @@ bool Foam::ensightOutput::writePointField
os.beginPart(part.index());
}
labelList uniqueMeshPointLabels;
part.uniqueMeshPoints(mesh, uniqueMeshPointLabels, parallel);
labelList uniquePointLabels;
part.uniqueMeshPoints(mesh, uniquePointLabels, parallel);
ensightOutput::Detail::writeFieldComponents
(
os,
ensightFile::coordinates,
UIndirectList<Type>(pf.internalField(), uniqueMeshPointLabels),
UIndirectList<Type>(pf.internalField(), uniquePointLabels),
parallel
);
}
@ -249,30 +249,42 @@ bool Foam::ensightOutput::writePointField
os.beginPart(part.index());
}
labelList uniquePointLabels;
part.uniqueMeshPoints(mesh, uniquePointLabels, parallel);
const auto& bfld = pf.boundaryField()[patchId];
// Only valuePointPatchField is actually derived from Field
const auto* vpp = isA<Field<Type>>(bfld);
if (vpp)
{
ensightOutput::Detail::writeFieldComponents
(
os,
ensightFile::coordinates,
*vpp,
parallel
);
}
else
{
labelList uniqueMeshPointLabels;
part.uniqueMeshPoints(mesh, uniqueMeshPointLabels, parallel);
// Require patch local indices, not mesh point labels.
// But need to use polyPatch meshPointMap() to recover the
// local indices since the ensight output will have jumbled
// the face output order
const polyPatch& pp = mesh.boundaryMesh()[patchId];
for (label& pointi : uniquePointLabels)
{
pointi = pp.meshPointMap()[pointi];
}
ensightOutput::Detail::writeFieldComponents
(
os,
ensightFile::coordinates,
UIndirectList<Type>(pf.internalField(), uniqueMeshPointLabels),
UIndirectList<Type>(*vpp, uniquePointLabels),
parallel
);
}
else
{
ensightOutput::Detail::writeFieldComponents
(
os,
ensightFile::coordinates,
UIndirectList<Type>(pf.internalField(), uniquePointLabels),
parallel
);
}
@ -294,14 +306,14 @@ bool Foam::ensightOutput::writePointField
// uses internalField only
{
labelList uniqueMeshPointLabels;
part.uniqueMeshPoints(mesh, uniqueMeshPointLabels, parallel);
labelList uniquePointLabels;
part.uniqueMeshPoints(mesh, uniquePointLabels, parallel);
ensightOutput::Detail::writeFieldComponents
(
os,
ensightFile::coordinates,
UIndirectList<Type>(pf.internalField(), uniqueMeshPointLabels),
UIndirectList<Type>(pf.internalField(), uniquePointLabels),
parallel
);
}