Commit da400179 authored by thery's avatar thery

bug in Traversor2FFaV

parent 29179079
...@@ -241,6 +241,7 @@ public: ...@@ -241,6 +241,7 @@ public:
} ; } ;
// Traverse the faces adjacent to a given face through sharing a common edge // Traverse the faces adjacent to a given face through sharing a common edge
// Warning mult-incidence is not managed (some faces can be send several times)
template <typename MAP> template <typename MAP>
class Traversor2FFaE //: public Traversor<MAP> class Traversor2FFaE //: public Traversor<MAP>
{ {
......
...@@ -530,12 +530,12 @@ Traversor2FFaV<MAP>::Traversor2FFaV(const MAP& map, Dart dart) : m(map),m_QLT(NU ...@@ -530,12 +530,12 @@ Traversor2FFaV<MAP>::Traversor2FFaV(const MAP& map, Dart dart) : m(map),m_QLT(NU
else else
{ {
start = m.phi2(m.phi_1(m.phi2(m.phi_1(dart)))) ; start = m.phi2(m.phi_1(m.phi2(m.phi_1(dart)))) ;
current = start ; while (start == dart)
if(start == dart)
{ {
stop = m.phi2(m.phi_1(dart)) ; dart = m.phi1(dart);
start = next() ; start = m.phi2(m.phi_1(m.phi2(m.phi_1(dart)))) ;
} }
current = start ;
stop = dart ; stop = dart ;
if(m.isBoundaryMarked2(start)) if(m.isBoundaryMarked2(start))
start = next() ; start = next() ;
...@@ -577,12 +577,15 @@ Dart Traversor2FFaV<MAP>::next() ...@@ -577,12 +577,15 @@ Dart Traversor2FFaV<MAP>::next()
current = m.phi2(m.phi_1(m.phi2(m.phi_1(d)))) ; current = m.phi2(m.phi_1(m.phi2(m.phi_1(d)))) ;
if(current == d) if(current == d)
{ {
stop = m.phi2(m.phi_1(d)) ; stop = m.phi1(d);
current = m.phi2(d);
return next() ; return next() ;
} }
stop = d ; stop = d ;
if(m.isBoundaryMarked2(current)) }
return next() ; if(m.isBoundaryMarked2(current))
{
return next() ;
} }
if(current == start) if(current == start)
current = NIL ; current = NIL ;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment