From 7bf524013921b8b5bd58de72719a7d661e14bced Mon Sep 17 00:00:00 2001
From: jfa <jfa@opencascade.com>
Date: Wed, 21 Apr 2010 12:46:40 +0000
Subject: [PATCH] Mantis issue 0020422: EDF 1070 GEOM: Partition anomaly -
 REGRESSION from 4.1.x.

---
 .../StdMeshers_CompositeSegment_1D.cxx        | 29 +++++++++++++++----
 1 file changed, 23 insertions(+), 6 deletions(-)

diff --git a/src/StdMeshers/StdMeshers_CompositeSegment_1D.cxx b/src/StdMeshers/StdMeshers_CompositeSegment_1D.cxx
index eee3ef012..a3a78e6ec 100644
--- a/src/StdMeshers/StdMeshers_CompositeSegment_1D.cxx
+++ b/src/StdMeshers/StdMeshers_CompositeSegment_1D.cxx
@@ -81,7 +81,16 @@ namespace {
     TopTools_MapOfShape edgeCounter;
     edgeCounter.Add( edge );
     TopoDS_Vertex v;
-    v = forward ? TopExp::LastVertex( edge,1 ) : TopExp::FirstVertex( edge,1 );
+    if (edge.Orientation() <= TopAbs_REVERSED) // FORWARD || REVERSED
+      v = forward ? TopExp::LastVertex(edge,true) : TopExp::FirstVertex(edge,true);
+    else {
+      TopoDS_Iterator vIt (edge);
+      if (forward)
+        for (; vIt.More(); vIt.Next())
+          v = TopoDS::Vertex(vIt.Value()); // last
+      else if (vIt.More())
+        v = TopoDS::Vertex(vIt.Value()); // first
+    }
     TopTools_ListIteratorOfListOfShape ancestIt = aMesh.GetAncestors( v );
     for ( ; ancestIt.More(); ancestIt.Next() )
     {
@@ -92,11 +101,19 @@ namespace {
     if ( edgeCounter.Extent() < 3 && !eNext.IsNull() ) {
       if ( SMESH_Algo::IsContinuous( edge, eNext )) {
         // care of orientation
-        bool reverse;
-        if ( forward )
-          reverse = ( !v.IsSame( TopExp::FirstVertex( eNext, true )));
-        else
-          reverse = ( !v.IsSame( TopExp::LastVertex( eNext, true )));
+        TopoDS_Vertex vn;
+        if (eNext.Orientation() <= TopAbs_REVERSED) // FORWARD || REVERSED
+          vn = forward ? TopExp::FirstVertex(eNext,true) : TopExp::LastVertex(eNext,true);
+        else {
+          TopoDS_Iterator vIt (eNext);
+          if (forward)
+            vn = TopoDS::Vertex(vIt.Value()); // first
+          else if (vIt.More())
+            for (; vIt.More(); vIt.Next())
+              vn = TopoDS::Vertex(vIt.Value()); // last
+        }
+
+        bool reverse = (!v.IsSame(vn));
         if ( reverse )
           eNext.Reverse();
         return eNext;
-- 
2.30.2