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

Bugfix: empty wedge patch in parallel. Zeljko Tukovic

Browse source
This commit is contained in:
Hrvoje Jasak 2017-05-25 14:40:47 +01:00
parent 05899eaa9e
commit 01aebe4bf9
1 changed files with 20 additions and 16 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

36
src/foam/fields/PointPatchFieldTemplates/constraint/wedge/WedgePointPatchField.C
View file

@ -202,23 +202,27 @@ WedgePointPatchField
const Pstream::commsTypes commsType
)
{
// In order to ensure that the wedge patch is always flat, take the
// normal vector from the first point
const vector& nHat = this->patch().pointNormals()[0];
tmp<Field<Type> > tvalues =
transform(I - nHat*nHat, this->patchInternalField());
const Field<Type>& values = tvalues();
// Get internal field to insert values into
Field<Type>& iF = const_cast<Field<Type>&>(this->internalField());
// Get addressing
const labelList& meshPoints = this->patch().meshPoints();
forAll (meshPoints, pointI)
// ZT, 26/02/2017: Size of the patch could be zero in parallel runs
if (!this->patch().meshPoints().empty())
{
iF[meshPoints[pointI]] = values[pointI];
// In order to ensure that the wedge patch is always flat, take the
// normal vector from the first point
const vector& nHat = this->patch().pointNormals()[0];
tmp<Field<Type> > tvalues =
transform(I - nHat*nHat, this->patchInternalField());
const Field<Type>& values = tvalues();
// Get internal field to insert values into
Field<Type>& iF = const_cast<Field<Type>&>(this->internalField());
// Get addressing
const labelList& meshPoints = this->patch().meshPoints();
forAll (meshPoints, pointI)
{
iF[meshPoints[pointI]] = values[pointI];
}
}
}