BUGFIX: Correct interpolation of point data on processor boundaries after parallel topological changes. Author: Zeljko Tukovic. Merge: Dominik Christ.

src/foam/fields/PointPatchFields/constraint/processor/ProcessorPointPatchField.C

@@ -511,7 +511,18 @@ evaluate
             }
 
             // Average over two sides
-            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<Type>(tpn, neiPoints)
+                );
 
             // Get internal field to insert values into
             Field<Type>& iF = const_cast<Field<Type>&>(this->internalField());