Add face numbering

The approach is very simple. It is incorrect in parallel.

src/mesh/Connectivity.cpp

@@ -142,14 +142,13 @@ void Connectivity<3>::_computeCellFaceAndFaceNodeConnectivities()
   }
 
   {
-    int l=0;
+    FaceId l=0;
     for (const auto& face_cells_vector : face_cells_map) {
       const Face& face = face_cells_vector.first;
       m_face_number_map[face] = l;
       ++l;
     }
   }
-
 #warning check that the number of cell per faces is <=2
 }
 
@@ -279,6 +278,16 @@ Connectivity(const ConnectivityDescriptor& descriptor)
 
   if constexpr (Dimension>1) {
     this->_computeCellFaceAndFaceNodeConnectivities();
+
+    {
+#warning incorrect in parallel
+      FaceValue<int> face_number(*this);
+
+      parallel_for (face_number.size(), PASTIS_LAMBDA(const FaceId& l) {
+          face_number[l] = l;
+        });
+      m_face_number = face_number;
+    }
   }
 }
 

src/mesh/Connectivity.hpp

@@ -695,6 +695,7 @@ class Connectivity final
     return m_inv_cell_nb_nodes;
   }
 
+#warning remove this member
   unsigned int getFaceNumber(const std::vector<unsigned int>& face_nodes) const
   {
     const Face face(face_nodes);