Commit 4807448b authored by Henry Weller's avatar Henry Weller
Browse files

cyclicACMIPolyPatch: Clear geometry when resetting

Patch contributed by Mattijs Janssens
parent 0d7bf4ff
...@@ -60,16 +60,25 @@ void Foam::cyclicACMIPolyPatch::resetAMI ...@@ -60,16 +60,25 @@ void Foam::cyclicACMIPolyPatch::resetAMI
if (boundaryMesh().mesh().hasCellCentres()) if (boundaryMesh().mesh().hasCellCentres())
{ {
WarningInFunction if (debug)
<< "The mesh already has cellCentres calculated when" {
<< " resetting ACMI " << name() << "." << endl Pout<< "cyclicACMIPolyPatch::resetAMI : clearing cellCentres"
<< "This is a problem since ACMI adapts the face areas" << " for " << name() << " and " << nonOverlapPatch.name()
<< " (to close cells) so this has" << endl << endl;
<< "to be done before cell centre calculation." << endl }
<< "This can happen if e.g. the cyclicACMI is after"
<< " any processor patches in the boundary." << endl //WarningInFunction
<< "Continuing with potential for incorrect geometry" // << "The mesh already has cellCentres calculated when"
<< " calculation and mass loss" << endl; // << " resetting ACMI " << name() << "." << endl
// << "This is a problem since ACMI adapts the face areas"
// << " (to close cells) so this has" << endl
// << "to be done before cell centre calculation." << endl
// << "This can happen if e.g. the cyclicACMI is after"
// << " any processor patches in the boundary." << endl;
const_cast<polyMesh&>
(
boundaryMesh().mesh()
).primitiveMesh::clearGeom();
} }
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment