restore oldi-school embnull testing algos

include/Algo/MC/marchingcube.hpp

@@ -26,8 +26,6 @@
 #include "Topology/generic/dartmarker.h"
 #include <vector>
 
-#include "Algo/Modelisation/boundEmb.h"
-
 namespace CGoGN
 {
 
@@ -119,9 +117,7 @@ void MarchingCube<DataType, Windowing, PFP>::deleteMesh()
 template< typename  DataType, template < typename D2 > class Windowing, typename PFP >
 Dart  MarchingCube<DataType, Windowing, PFP>::createTriEmb(unsigned int e1, unsigned int e2, unsigned int e3)
 {
-//	L_DART d = m_map->newFace(3,false);
-
-	L_DART d = Modelisation::newFaceEmb<PFP>(*m_map,3);
+	L_DART d = m_map->newFace(3,false);
 		
 	FunctorSetEmb<GenericMap, VERTEX> fsetemb(*m_map, e1);
 	m_map->template foreach_dart_of_orbit<PFP::MAP::VERTEX_OF_PARENT>(d, fsetemb);
@@ -926,9 +922,7 @@ void MarchingCube<DataType, Windowing, PFP>::setNeighbourSimple(L_DART d1, L_DAR
 {
 	if (m_map->phi2(d1) != d2)
 	{
-//		m_map->sewFaces(d1,d2,false);
-		Modelisation::sewFaceEmb<PFP>(*m_map,d1,d2);
-
+		m_map->sewFaces(d1,d2,false);
 		#ifdef MC_WIDTH_EDGE_Z_EMBEDED
 		if (m_zslice!=NULL)
 		{
@@ -945,9 +939,7 @@ void MarchingCube<DataType, Windowing, PFP>::setNeighbour(L_DART d1, L_DART d2)
 {
 	if (m_map->phi2(d1) != d2)
 	{
-//		m_map->sewFaces(d1,d2,false);
-		Modelisation::sewFaceEmb<PFP>(*m_map,d1,d2);
-
+		m_map->sewFaces(d1,d2,false);
 		#ifdef MC_WIDTH_EDGE_Z_EMBEDED
 		if (m_zslice!=NULL)
 		{

include/Algo/Modelisation/polyhedron.h

@@ -200,7 +200,8 @@ protected:
 
 	void computeCenter();
 
-
+	Dart grid_topo_open(unsigned int x, unsigned int y);
+	Dart cylinder_topo_open(unsigned int n, unsigned int z);
 public:
 	/**
 	* Constructor
@@ -239,6 +240,8 @@ public:
 	*/
 	std::vector<Dart>& getVertexDarts() { return m_tableVertDarts; }
 
+
+
 	/**
 	* Create a 2D grid
 	* quads are oriented counter-clockwise and the returned dart