From 14dce8cef92362a7b86cc43f7eaab70e1ea03edf Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Tue, 29 May 2012 11:04:49 +0100
Subject: [PATCH] BUG: extendedFeatureEdgeMesh: updated indexing of edges

---
 .../extendedFeatureEdgeMesh.C                 | 20 +++++++++++++++++--
 1 file changed, 18 insertions(+), 2 deletions(-)

diff --git a/src/edgeMesh/extendedFeatureEdgeMesh/extendedFeatureEdgeMesh.C b/src/edgeMesh/extendedFeatureEdgeMesh/extendedFeatureEdgeMesh.C
index 7899c7743a..76d336c8ea 100644
--- a/src/edgeMesh/extendedFeatureEdgeMesh/extendedFeatureEdgeMesh.C
+++ b/src/edgeMesh/extendedFeatureEdgeMesh/extendedFeatureEdgeMesh.C
@@ -1189,8 +1189,24 @@ void Foam::extendedFeatureEdgeMesh::add(const extendedFeatureEdgeMesh& fem)
     }
 
     edgeList newEdges(newEdgeI);
-    UIndirectList<edge>(newEdges, reverseEdgeMap) = edges();
-    UIndirectList<edge>(newEdges, reverseFemEdgeMap) = fem.edges();
+    forAll(edges(), i)
+    {
+        const edge& e = edges()[i];
+        newEdges[reverseEdgeMap[i]] = edge
+        (
+            reversePointMap[e[0]],
+            reversePointMap[e[1]]
+        );
+    }
+    forAll(fem.edges(), i)
+    {
+        const edge& e = fem.edges()[i];
+        newEdges[reverseFemEdgeMap[i]] = edge
+        (
+            reverseFemPointMap[e[0]],
+            reverseFemPointMap[e[1]]
+        );
+    }
 
     pointField newEdgeDirections(newEdgeI);
     newEdgeDirections.rmap(edgeDirections(), reverseEdgeMap);