changement mineurs particules a tester vertexstate et modif subdivide

include/Algo/ImplicitHierarchicalMesh/subdivision3.hpp

@@ -27,7 +27,7 @@
 #include "Algo/Modelisation/extrusion.h"
 #include "Topology/generic/dartmarker.h"
 
-#define DEBUG_affichage
+//#define DEBUG_affichage
 namespace CGoGN
 {
 
@@ -214,13 +214,15 @@ Dart subdivideVolumeClassic(typename PFP::MAP& map, Dart d, VertexAttribute<type
 
 	if(OneLevelDifference)
 	{
+
 		Traversor3WF<typename PFP::MAP> traF(map, old);
 		for(Dart dit = traF.begin(); dit != traF.end(); dit = traF.next())
 		{
 			Dart nv = map.phi3(dit);
 			if(!map.isBoundaryMarked(3, nv))
 				if(map.volumeLevel(nv) == vLevel - 1)
-					subdivideVolumeClassic<PFP>(map,nv,position,OneLevelDifference);
+
+                    subdivideVolumeClassic<PFP>(map,nv,position,OneLevelDifference);
 		}
 	}
 
@@ -354,7 +356,7 @@ Dart subdivideVolumeClassic(typename PFP::MAP& map, Dart d, VertexAttribute<type
 
 //            v.push_back(map.phi1(e));
 //            v.push_back(map.phi1(map.phi1(e)));
-			std::cout << "v size = " << v.size() << std::endl;
+//			std::cout << "v size = " << v.size() << std::endl;
             e = map.phi2(map.phi_1(e));
         }
         while(e != *edge);

include/Algo/MovingObjects/particle_cell_3D.hpp

@@ -366,7 +366,7 @@ void ParticleCell3D<PFP>::vertexState(const VEC3& current)
 #ifdef DEBUG
     std::cout << "points equal vrai :" << current << " || "<<this->m_position<<" || "<<d<<" || "<<position[d]<< std::endl;
 #endif
-//		this->m_position = som;
+        this->m_position = som;
         this->m_positionFace = som;
 		state = VERTEX;
 		return;
@@ -535,11 +535,14 @@ void ParticleCell3D<PFP>::edgeState(const VEC3& current)
 				}
                 else {  // sur le sommet de l'arete
                     this->m_position = position[2];
+                    std::cout << "edgeState va dans sommet cas 1 position :"<< this->m_position << std::endl;
+
 					vertexState(current);
 				}
             }
             else {  // sur l'autre sommet de l'arete
                 d=m.phi1(d);
+                std::cout << "edgeState va dans sommet cas 1 position :"<< this->m_position << std::endl;
 				vertexState(current);
 			}
 		}