From b5d1c3029964e57f4f824ae1dda68acf82ab724b Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Sun, 29 Jan 2012 11:45:46 +0000
Subject: [PATCH] Handle bridge overlap in region couple

---
 .../regionCoupling/regionCouplingFvPatchField.C     | 13 +++++++++++++
 1 file changed, 13 insertions(+)

diff --git a/src/finiteVolume/fields/fvPatchFields/constraint/regionCoupling/regionCouplingFvPatchField.C b/src/finiteVolume/fields/fvPatchFields/constraint/regionCoupling/regionCouplingFvPatchField.C
index 1a94df04e..76ed16acc 100644
--- a/src/finiteVolume/fields/fvPatchFields/constraint/regionCoupling/regionCouplingFvPatchField.C
+++ b/src/finiteVolume/fields/fvPatchFields/constraint/regionCoupling/regionCouplingFvPatchField.C
@@ -353,6 +353,19 @@ void regionCouplingFvPatchField<Type>::initEvaluate
 
     // Do interpolation
     Field<Type>::operator=(weights*fOwn + (1.0 - weights)*fNei);
+
+    if (regionCouplePatch_.bridgeOverlap())
+    {
+        // Symmetry treatment used for overlap
+        vectorField nHat = this->patch().nf();
+
+        Field<Type> pif = this->patchInternalField();
+
+        Field<Type> bridgeField =
+            0.5*(pif + transform(I - 2.0*sqr(nHat), pif));
+
+        regionCouplePatch_.bridge(bridgeField, *this);
+    }
 }