From f1e7cae6469c753042200d12a6c853646c7cfb4b Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 12 Apr 2010 11:48:04 +0100
Subject: [PATCH] ENH: Have edge addressing from coupledPatch to mesh

---
 .../polyMesh/globalMeshData/globalMeshData.C  | 42 +++++++++++++++++++
 .../polyMesh/globalMeshData/globalMeshData.H  | 11 +++++
 2 files changed, 53 insertions(+)

diff --git a/src/OpenFOAM/meshes/polyMesh/globalMeshData/globalMeshData.C b/src/OpenFOAM/meshes/polyMesh/globalMeshData/globalMeshData.C
index 88d02457de..17c8254387 100644
--- a/src/OpenFOAM/meshes/polyMesh/globalMeshData/globalMeshData.C
+++ b/src/OpenFOAM/meshes/polyMesh/globalMeshData/globalMeshData.C
@@ -1122,6 +1122,9 @@ void Foam::globalMeshData::clearOut()
     deleteDemandDrivenData(sharedEdgeAddrPtr_);
 
     coupledPatchPtr_.clear();
+    coupledPatchMeshEdgesPtr_.clear();
+    coupledPatchMeshEdgeMapPtr_.clear();
+
     // Point
     globalPointNumberingPtr_.clear();
     globalPointSlavesPtr_.clear();
@@ -1409,6 +1412,45 @@ const Foam::indirectPrimitivePatch& Foam::globalMeshData::coupledPatch() const
 }
 
 
+const Foam::labelList& Foam::globalMeshData::coupledPatchMeshEdges() const
+{
+    if (!coupledPatchMeshEdgesPtr_.valid())
+    {
+        coupledPatchMeshEdgesPtr_.reset
+        (
+            new labelList
+            (
+                coupledPatch().meshEdges
+                (
+                    mesh_.edges(),
+                    mesh_.pointEdges()
+                )
+            )
+        );
+    }
+    return coupledPatchMeshEdgesPtr_();
+}
+
+
+const Foam::Map<Foam::label>& Foam::globalMeshData::coupledPatchMeshEdgeMap()
+const
+{
+    if (!coupledPatchMeshEdgeMapPtr_.valid())
+    {
+        const labelList& me = coupledPatchMeshEdges();
+
+        coupledPatchMeshEdgeMapPtr_.reset(new Map<label>(2*me.size()));
+        Map<label>& em = coupledPatchMeshEdgeMapPtr_();
+
+        forAll(me, i)
+        {
+            em.insert(me[i], i);
+        }
+    }
+    return coupledPatchMeshEdgeMapPtr_();
+}
+
+
 const Foam::globalIndex& Foam::globalMeshData::globalPointNumbering() const
 {
     if (!globalPointNumberingPtr_.valid())
diff --git a/src/OpenFOAM/meshes/polyMesh/globalMeshData/globalMeshData.H b/src/OpenFOAM/meshes/polyMesh/globalMeshData/globalMeshData.H
index fe2e04731b..79e086fd23 100644
--- a/src/OpenFOAM/meshes/polyMesh/globalMeshData/globalMeshData.H
+++ b/src/OpenFOAM/meshes/polyMesh/globalMeshData/globalMeshData.H
@@ -208,6 +208,10 @@ class globalMeshData
 
             mutable autoPtr<indirectPrimitivePatch> coupledPatchPtr_;
 
+            mutable autoPtr<labelList> coupledPatchMeshEdgesPtr_;
+
+            mutable autoPtr<Map<label> > coupledPatchMeshEdgeMapPtr_;
+
             // Collocated
 
                 // Coupled point to collocated coupled points
@@ -506,6 +510,13 @@ public:
             //- Return patch of all coupled faces
             const indirectPrimitivePatch& coupledPatch() const;
 
+            //- Return map from coupledPatch edges to mesh edges
+            const labelList& coupledPatchMeshEdges() const;
+
+            //- Return map from mesh edges to coupledPatch edges
+            const Map<label>& coupledPatchMeshEdgeMap() const;
+
+
             // Coupled point to collocated coupled points. Coupled points are
             // points on any coupled patch.