zhyang-dev/OpenFOAM-6
1
0
Fork 0
mirror of https://github.com/OpenFOAM/OpenFOAM-6.git synced 2025-12-08 06:57:46 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

cyclicRepeatAMI: Added error message for zero-sized transform patch

Browse Source 
Sometimes decomposition can remove a cyclic patch entirely, converting
it into a processor cyclic. If this patch is being used to specify the
transform for a cyclicRepeatAMI patch then the calculation will fail.
This change adds a check for this situation, and errors with the
suggestion that the transform patch be preserved during decomposition.
This commit is contained in:
Will Bainbridge
2018-06-04 10:47:49 +01:00
parent 7a93d91a80
commit 1223da98f7
1 changed files with 14 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
src/meshTools/AMIInterpolation/patches/cyclicRepeatAMI/cyclicRepeatAMIPolyPatch/cyclicRepeatAMIPolyPatch.C
View File

@ -82,8 +82,9 @@ void Foam::cyclicRepeatAMIPolyPatch::resetAMI() const
} }
Info<< indent << typeName <<" : Creating addressing and weights between " Info<< indent << typeName <<" : Creating addressing and weights between "
<< this->size() << " source faces and " << neighbPatch().size() << returnReduce(this->size(), sumOp<label>()) << " source faces and "
<< " target faces" << endl; << returnReduce(neighbPatch().size(), sumOp<label>()) << " target faces"
<< endl;
// Get the transform associated with the transform patch // Get the transform associated with the transform patch
vectorTensorTransform t; vectorTensorTransform t;
@ -129,6 +130,17 @@ void Foam::cyclicRepeatAMIPolyPatch::resetAMI() const
reduce(bt, keepIfTrueOp<vectorTensorTransform>()); reduce(bt, keepIfTrueOp<vectorTensorTransform>());
if (!bt.first())
{
FatalErrorInFunction
<< "Transform patch " << transformPatch.name() << " for "
<< typeName << " patch " << name() << " has zero faces. It may "
<< "have been converted to a processor cyclic during "
<< "decomposition. Consider adding " << transformPatch.name()
<< " and it's neighbour to the list of preserved patches."
<< exit(FatalError);
}
t = bt.second(); t = bt.second();
} }
const vectorTensorTransform tInv(inv(t)); const vectorTensorTransform tInv(inv(t));