From 0e113e1a5de95f9e91eed5f32e882b9d6b157517 Mon Sep 17 00:00:00 2001 From: Dominik Christ Date: Wed, 22 Apr 2015 18:20:22 +0100 Subject: [PATCH] Correct averaging in ProcessorPointPatchField --- ListOfContributors | 1 + .../processor/ProcessorPointPatchField.C | 15 +++++++++------ 2 files changed, 10 insertions(+), 6 deletions(-) diff --git a/ListOfContributors b/ListOfContributors index 7f813b0ac..c9cf38af8 100644 --- a/ListOfContributors +++ b/ListOfContributors @@ -73,3 +73,4 @@ Tian Tang Zeljko Tukovic Niklas Wikstrom + Vanja Skuric diff --git a/src/foam/fields/PointPatchFields/constraint/processor/ProcessorPointPatchField.C b/src/foam/fields/PointPatchFields/constraint/processor/ProcessorPointPatchField.C index e119c8dc2..b76b37208 100644 --- a/src/foam/fields/PointPatchFields/constraint/processor/ProcessorPointPatchField.C +++ b/src/foam/fields/PointPatchFields/constraint/processor/ProcessorPointPatchField.C @@ -511,18 +511,21 @@ evaluate } // Average over two sides + // VS, 2015-04-12 - This works in parallel! + tpn = 0.5*(this->patchInternalField(this->internalField()) + tpn); // ZT, 22-07-2014 - point ordering is not same // at master and slave side after topology change const labelList& neiPoints = procPatch_.procPolyPatch().neighbPoints(); - tpn = - 0.5* - ( - this->patchInternalField(this->internalField()) - + Field(tpn, neiPoints) - ); + // VS, 2015-04-12 - This doesn't work in parallel! +// tpn = +// 0.5* +// ( +// this->patchInternalField(this->internalField()) +// + Field(tpn, neiPoints) +// ); // Get internal field to insert values into Field& iF = const_cast&>(this->internalField());