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(newEdges, reverseEdgeMap) = edges(); - UIndirectList(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);