diff --git a/include/Algo/ImplicitHierarchicalMesh/subdivision.hpp b/include/Algo/ImplicitHierarchicalMesh/subdivision.hpp
index c0b35f0158e6d918559343dc00b2dc5955ef5cba..488a92b49543a2191c6b4a831dab52a4db11206c 100644
--- a/include/Algo/ImplicitHierarchicalMesh/subdivision.hpp
+++ b/include/Algo/ImplicitHierarchicalMesh/subdivision.hpp
@@ -92,38 +92,64 @@ void subdivideFace(typename PFP::MAP& map, Dart d, VertexAttribute<typename PFP:
 		Dart dd = map.phi1(old) ;
 		Dart e = map.phi1(map.phi1(dd)) ;
 		map.splitFace(dd, e) ;					// insert a new edge
-		unsigned int id = map.getNewEdgeId() ;
-		map.setEdgeId(map.phi_1(dd), id) ;		// set the edge id of the inserted
-		map.setEdgeId(map.phi_1(e), id) ;		// edge to the next available id
+		//unsigned int id = map.getNewEdgeId() ;
+		//map.setEdgeId(map.phi_1(dd), id) ;		// set the edge id of the inserted
+		//map.setEdgeId(map.phi_1(e), id) ;		// edge to the next available id
 
 		dd = e ;
 		e = map.phi1(map.phi1(dd)) ;
 		map.splitFace(dd, e) ;
-		id = map.getNewEdgeId() ;
-		map.setEdgeId(map.phi_1(dd), id) ;
-		map.setEdgeId(map.phi_1(e), id) ;
+		//id = map.getNewEdgeId() ;
+		//map.setEdgeId(map.phi_1(dd), id) ;
+		//map.setEdgeId(map.phi_1(e), id) ;
 
 		dd = e ;
 		e = map.phi1(map.phi1(dd)) ;
 		map.splitFace(dd, e) ;
-		id = map.getNewEdgeId() ;
-		map.setEdgeId(map.phi_1(dd), id) ;
-		map.setEdgeId(map.phi_1(e), id) ;
+		//id = map.getNewEdgeId() ;
+		//map.setEdgeId(map.phi_1(dd), id) ;
+		//map.setEdgeId(map.phi_1(e), id) ;
+
+
+		Dart dit = map.phi2(map.phi_1(dd));
+		do
+		{
+			unsigned int dId = map.getEdgeId(map.phi_1(map.phi2(dit)));
+			unsigned int eId = map.getEdgeId(map.phi1(map.phi2(dit)));
+
+			unsigned int t = dId + eId;
+
+			if(t == 1)
+			{
+				map.setEdgeId(dit, 2);
+			}
+			else if(t == 2)
+			{
+				map.setEdgeId(dit, 1);
+			}
+			else if(t == 3)
+			{
+				map.setEdgeId(dit, 0);
+			}
+
+			dit = map.phi1(dit);
+		}while(dit != map.phi2(map.phi_1(dd)));
+
 	}
 	else											// if subdividing a polygonal face
 	{
 		Dart dd = map.phi1(old) ;
 		map.splitFace(dd, map.phi1(map.phi1(dd))) ;	// insert a first edge
 		Dart ne = map.alpha1(dd) ;
-		Dart ne2 = map.phi2(ne) ;
+		//Dart ne2 = map.phi2(ne) ;
 
 		map.cutEdge(ne) ;							// cut the new edge to insert the central vertex
-		unsigned int id = map.getNewEdgeId() ;
-		map.setEdgeId(ne, id) ;
-		map.setEdgeId(map.phi2(ne), id) ;			// set the edge id of the inserted
-		id = map.getNewEdgeId() ;
-		map.setEdgeId(ne2, id) ;					// edges to the next available ids
-		map.setEdgeId(map.phi2(ne2), id) ;
+		//unsigned int id = map.getNewEdgeId() ;
+		//map.setEdgeId(ne, id) ;
+		//map.setEdgeId(map.phi2(ne), id) ;			// set the edge id of the inserted
+		//id = map.getNewEdgeId() ;
+		//map.setEdgeId(ne2, id) ;					// edges to the next available ids
+		//map.setEdgeId(map.phi2(ne2), id) ;
 
 		position[map.phi2(ne)] = p ;
 
@@ -133,11 +159,24 @@ void subdivideFace(typename PFP::MAP& map, Dart d, VertexAttribute<typename PFP:
 			Dart next = map.phi1(map.phi1(dd)) ;
 			map.splitFace(map.phi1(ne), dd) ;
 			Dart nne = map.alpha1(dd) ;
-			id = map.getNewEdgeId() ;
-			map.setEdgeId(nne, id) ;
-			map.setEdgeId(map.phi2(nne), id) ;
+			//id = map.getNewEdgeId() ;
+			//map.setEdgeId(nne, id) ;
+			//map.setEdgeId(map.phi2(nne), id) ;
 			dd = next ;
 		}
+
+		Dart dit = map.phi2(ne);
+		do
+		{
+			unsigned int eId = map.getEdgeId(map.phi1(dit));
+			if(eId == 0)
+				map.setEdgeId(dit, 1);
+			else if(eId == 1)
+				map.setEdgeId(dit, 0);
+
+			dit = map.phi2(map.phi_1(ne));
+		}
+		while(dit != map.phi2(ne));
 	}
 
 	map.setCurrentLevel(cur) ;
diff --git a/src/Algo/ImplicitHierarchicalMesh/ihm.cpp b/src/Algo/ImplicitHierarchicalMesh/ihm.cpp
index 649fe108f6dc442b80f62e456eb48ee54ad529ff..3c4392e48e6a37817182bf30ab9660800c2413ec 100644
--- a/src/Algo/ImplicitHierarchicalMesh/ihm.cpp
+++ b/src/Algo/ImplicitHierarchicalMesh/ihm.cpp
@@ -53,7 +53,14 @@ ImplicitHierarchicalMap::~ImplicitHierarchicalMap()
 
 void ImplicitHierarchicalMap::init()
 {
-	initEdgeId() ;
+	//initEdgeId() ;
+
+	//init each edge Id at 0
+	TraversorE<ImplicitHierarchicalMap> te(*this);
+	for(Dart dit = te.begin() ; dit != te.next() ; dit = te.next())
+	{
+		m_edgeId[dit] = 0;
+	}
 
 	for(unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
 	{