zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

createNonConformalCouples: Fix writing of non-conformal processor cyclic fields

Browse Source
This commit is contained in:
Will Bainbridge
2024-06-27 13:02:57 +01:00
parent 928db2bdce
commit 0dd55ebad0
1 changed files with 60 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

60
applications/utilities/mesh/manipulation/createNonConformalCouples/createNonConformalCouples.C
View File

@ -240,6 +240,51 @@ struct nonConformalCouple
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
template<class Type>
void evaluateNonConformalProcessorCyclics(const fvMesh& mesh)
{
UPtrList<VolField<Type>> fields(mesh.fields<VolField<Type>>());
forAll(fields, i)
{
const label nReq = Pstream::nRequests();
forAll(mesh.boundary(), patchi)
{
typename VolField<Type>::Patch& pf =
fields[i].boundaryFieldRef()[patchi];
if (isA<nonConformalProcessorCyclicPolyPatch>(pf.patch().patch()))
{
pf.initEvaluate(Pstream::defaultCommsType);
}
}
if
(
Pstream::parRun()
&& Pstream::defaultCommsType == Pstream::commsTypes::nonBlocking
)
{
Pstream::waitRequests(nReq);
}
forAll(mesh.boundary(), patchi)
{
typename VolField<Type>::Patch& pf =
fields[i].boundaryFieldRef()[patchi];
if (isA<nonConformalProcessorCyclicPolyPatch>(pf.patch().patch()))
{
pf.evaluate(Pstream::defaultCommsType);
}
}
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
int main(int argc, char *argv[])
@ -754,6 +799,21 @@ int main(int argc, char *argv[])
}
}
// Communicate values across non-conformal processor cyclics so that they
// contain valid values that can be written to disk
if (Pstream::parRun())
{
forAll(regionMeshes, regioni)
{
const fvMesh& mesh = regionMeshes[regioni];
#define EVALUATE_NON_CONFORMAL_PROCESSOR_CYCLICS(Type, nullArg) \
evaluateNonConformalProcessorCyclics<Type>(mesh);
FOR_ALL_FIELD_TYPES(EVALUATE_NON_CONFORMAL_PROCESSOR_CYCLICS)
#undef EVALUATE_NON_CONFORMAL_PROCESSOR_CYCLICS
}
}
// Set the precision of the points data to 10
IOstream::defaultPrecision(max(10u, IOstream::defaultPrecision()));