genericmap.cpp 13.8 KB
Newer Older
Pierre Kraemer's avatar
Pierre Kraemer committed
1 2 3
/*******************************************************************************
* CGoGN: Combinatorial and Geometric modeling with Generic N-dimensional Maps  *
* version 0.1                                                                  *
4
* Copyright (C) 2009-2012, IGG Team, LSIIT, University of Strasbourg           *
Pierre Kraemer's avatar
Pierre Kraemer committed
5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
*                                                                              *
* This library is free software; you can redistribute it and/or modify it      *
* under the terms of the GNU Lesser General Public License as published by the *
* Free Software Foundation; either version 2.1 of the License, or (at your     *
* option) any later version.                                                   *
*                                                                              *
* This library is distributed in the hope that it will be useful, but WITHOUT  *
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or        *
* FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License  *
* for more details.                                                            *
*                                                                              *
* You should have received a copy of the GNU Lesser General Public License     *
* along with this library; if not, write to the Free Software Foundation,      *
* Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301 USA.           *
*                                                                              *
20
* Web site: http://cgogn.unistra.fr/                                           *
Pierre Kraemer's avatar
Pierre Kraemer committed
21 22 23 24 25
* Contact information: cgogn@unistra.fr                                        *
*                                                                              *
*******************************************************************************/

#include "Topology/generic/genericmap.h"
26
#include "Topology/generic/attributeHandler.h"
27
#include "Topology/generic/traversor/traversorCell.h"
Sylvain Thery's avatar
Sylvain Thery committed
28

Pierre Kraemer's avatar
Pierre Kraemer committed
29 30
#include "Geometry/vector_gen.h"
#include "Geometry/matrix.h"
31
#include "Container/registered.h"
Pierre Kraemer's avatar
Pierre Kraemer committed
32

Sylvain Thery's avatar
Sylvain Thery committed
33 34
#include <algorithm>

Pierre Kraemer's avatar
Pierre Kraemer committed
35 36
namespace CGoGN
{
37

Sylvain Thery's avatar
Sylvain Thery committed
38 39 40 41 42 43
namespace Parallel
{
//int NumberOfThreads=1;
int NumberOfThreads = getSystemNumberOfCores();
}

44

Pierre Kraemer's avatar
Pierre Kraemer committed
45
std::map<std::string, RegisteredBaseAttribute*>* GenericMap::m_attributes_registry_map = NULL;
46

47
int GenericMap::m_nbInstances = 0;
48

49 50 51
std::vector< std::vector<Dart>* >* GenericMap::s_vdartsBuffers = NULL;

std::vector< std::vector<unsigned int>* >* GenericMap::s_vintsBuffers = NULL;
Sylvain Thery's avatar
Sylvain Thery committed
52

53
std::vector<GenericMap*>*  GenericMap::s_instances=NULL;
Sylvain Thery's avatar
Sylvain Thery committed
54 55


56

Sylvain Thery's avatar
Sylvain Thery committed
57
GenericMap::GenericMap()
Pierre Kraemer's avatar
Pierre Kraemer committed
58
{
59
	if(m_attributes_registry_map == NULL)
Pierre Kraemer's avatar
Pierre Kraemer committed
60
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
61
		m_attributes_registry_map = new std::map<std::string, RegisteredBaseAttribute*>;
62

Pierre Kraemer's avatar
Pierre Kraemer committed
63 64 65
		// register all known types
		registerAttribute<Dart>("Dart");
		registerAttribute<Mark>("Mark");
Pierre Kraemer's avatar
Pierre Kraemer committed
66

Pierre Kraemer's avatar
Pierre Kraemer committed
67
		registerAttribute<char>("char");
68 69 70
		registerAttribute<short>("short");
		registerAttribute<int>("int");
		registerAttribute<long>("long");
Pierre Kraemer's avatar
Pierre Kraemer committed
71

Pierre Kraemer's avatar
Pierre Kraemer committed
72
		registerAttribute<unsigned char>("unsigned char");
73 74 75 76 77 78
		registerAttribute<unsigned short>("unsigned short");
		registerAttribute<unsigned int>("unsigned int");
		registerAttribute<unsigned long>("unsigned long");

		registerAttribute<float>("float");
		registerAttribute<double>("double");
Pierre Kraemer's avatar
Pierre Kraemer committed
79

Pierre Kraemer's avatar
Pierre Kraemer committed
80 81 82
		registerAttribute<Geom::Vec2f>(Geom::Vec2f::CGoGNnameOfType());
		registerAttribute<Geom::Vec3f>(Geom::Vec3f::CGoGNnameOfType());
		registerAttribute<Geom::Vec4f>(Geom::Vec4f::CGoGNnameOfType());
Pierre Kraemer's avatar
Pierre Kraemer committed
83

Pierre Kraemer's avatar
Pierre Kraemer committed
84 85 86
		registerAttribute<Geom::Vec2d>(Geom::Vec2d::CGoGNnameOfType());
		registerAttribute<Geom::Vec3d>(Geom::Vec3d::CGoGNnameOfType());
		registerAttribute<Geom::Vec4d>(Geom::Vec4d::CGoGNnameOfType());
Pierre Kraemer's avatar
Pierre Kraemer committed
87

Pierre Kraemer's avatar
Pierre Kraemer committed
88 89
		registerAttribute<Geom::Matrix33f>(Geom::Matrix33f::CGoGNnameOfType());
		registerAttribute<Geom::Matrix44f>(Geom::Matrix44f::CGoGNnameOfType());
Pierre Kraemer's avatar
Pierre Kraemer committed
90

Pierre Kraemer's avatar
Pierre Kraemer committed
91 92 93 94
		registerAttribute<Geom::Matrix33d>(Geom::Matrix33d::CGoGNnameOfType());
		registerAttribute<Geom::Matrix44d>(Geom::Matrix44d::CGoGNnameOfType());
	}

Sylvain Thery's avatar
Sylvain Thery committed
95 96


Pierre Kraemer's avatar
Pierre Kraemer committed
97
	m_nbInstances++;
98 99
	if (s_instances==NULL)
		s_instances= new std::vector<GenericMap*>;
Pierre Kraemer's avatar
Pierre Kraemer committed
100

101
	s_instances->push_back(this);
102

103
	if (s_vdartsBuffers == NULL)
Sylvain Thery's avatar
Sylvain Thery committed
104
	{
105 106 107
		s_vdartsBuffers = new std::vector< std::vector<Dart>* >[NB_THREAD];
		s_vintsBuffers = new std::vector< std::vector<unsigned int>* >[NB_THREAD];
		for(unsigned int i = 0; i < NB_THREAD; ++i)
Sylvain Thery's avatar
Sylvain Thery committed
108 109 110
		{
			s_vdartsBuffers[i].reserve(8);
			s_vintsBuffers[i].reserve(8);
111
				// prealloc ?
Sylvain Thery's avatar
Sylvain Thery committed
112 113
		}
	}
114 115 116 117 118 119 120 121


	for(unsigned int i = 0; i < NB_ORBITS; ++i)
	{
		m_attribs[i].setOrbit(i) ;
		m_attribs[i].setRegistry(m_attributes_registry_map) ;
	}

Pierre Kraemer's avatar
Pierre Kraemer committed
122
	init();
Pierre Kraemer's avatar
Pierre Kraemer committed
123 124 125 126
}

GenericMap::~GenericMap()
{
127

Pierre Kraemer's avatar
Pierre Kraemer committed
128 129 130 131 132
	for(unsigned int i = 0; i < NB_ORBITS; ++i)
	{
		if(isOrbitEmbedded(i))
			m_attribs[i].clear(true) ;
	}
133 134 135 136 137

	for(std::multimap<AttributeMultiVectorGen*, AttributeHandlerGen*>::iterator it = attributeHandlers.begin(); it != attributeHandlers.end(); ++it)
		(*it).second->setInvalid() ;
	attributeHandlers.clear() ;

138 139 140

	// clean type registry if necessary
	m_nbInstances--;
Pierre Kraemer's avatar
Pierre Kraemer committed
141
	if (m_nbInstances <= 0)
142
	{
143 144 145
		for (std::map<std::string, RegisteredBaseAttribute*>::iterator it =  m_attributes_registry_map->begin(); it != m_attributes_registry_map->end(); ++it)
			delete it->second;

146
		delete m_attributes_registry_map;
147
		m_attributes_registry_map = NULL;
Sylvain Thery's avatar
Sylvain Thery committed
148 149 150 151 152 153 154 155 156 157

		for(unsigned int i = 0; i < NB_THREAD; ++i)
		{
			for (auto it =s_vdartsBuffers[i].begin(); it != s_vdartsBuffers[i].end(); ++it)
				delete *it;
			for (auto it =s_vintsBuffers[i].begin(); it != s_vintsBuffers[i].end(); ++it)
				delete *it;
		}


158
	}
159

160 161 162 163
	// remove instance of table
	auto it = std::find(s_instances->begin(), s_instances->end(), this);
	*it = s_instances->back();
	s_instances->pop_back();
Pierre Kraemer's avatar
Pierre Kraemer committed
164 165
}

Pierre Kraemer's avatar
Pierre Kraemer committed
166
void GenericMap::init()
167
{
Pierre Kraemer's avatar
Pierre Kraemer committed
168
	for(unsigned int i = 0; i < NB_ORBITS; ++i)
169
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
170 171 172 173 174
		m_attribs[i].clear(true) ;
		m_embeddings[i] = NULL ;
		m_quickTraversal[i] = NULL;

		for(unsigned int j = 0; j < NB_ORBITS; ++j)
175
		{
Pierre Kraemer's avatar
Pierre Kraemer committed
176 177
			m_quickLocalIncidentTraversal[i][j] = NULL ;
			m_quickLocalAdjacentTraversal[i][j] = NULL ;
178
		}
Pierre Kraemer's avatar
Pierre Kraemer committed
179 180
	}

Sylvain Thery's avatar
Sylvain Thery committed
181 182
	m_boundaryMarkers[0] = m_attribs[DART].addAttribute<MarkerBool>("BoundaryMark0") ;
	m_boundaryMarkers[1] = m_attribs[DART].addAttribute<MarkerBool>("BoundaryMark1") ;
Pierre Kraemer's avatar
Pierre Kraemer committed
183

Sylvain Thery's avatar
Sylvain Thery committed
184
	 for(std::multimap<AttributeMultiVectorGen*, AttributeHandlerGen*>::iterator it = attributeHandlers.begin(); it != attributeHandlers.end(); ++it)
Pierre Kraemer's avatar
Pierre Kraemer committed
185 186 187 188 189 190 191 192
		(*it).second->setInvalid() ;
	attributeHandlers.clear() ;
}

void GenericMap::clear(bool removeAttrib)
{
	if (removeAttrib)
		init();
Sylvain Thery's avatar
Sylvain Thery committed
193 194 195 196 197
	else
	{
		for(unsigned int i = 0; i < NB_ORBITS; ++i)
			m_attribs[i].clear(false) ;
	}
untereiner's avatar
untereiner committed
198 199
}

Pierre Kraemer's avatar
Pierre Kraemer committed
200 201 202 203
/****************************************
 *        ATTRIBUTES MANAGEMENT         *
 ****************************************/

untereiner's avatar
untereiner committed
204 205 206 207 208 209 210 211 212 213 214 215
void GenericMap::swapEmbeddingContainers(unsigned int orbit1, unsigned int orbit2)
{
	assert(orbit1 != orbit2 || !"Cannot swap a container with itself") ;
	assert((orbit1 != DART && orbit2 != DART) || !"Cannot swap the darts container") ;

	m_attribs[orbit1].swap(m_attribs[orbit2]) ;
	m_attribs[orbit1].setOrbit(orbit1) ;	// to update the orbit information
	m_attribs[orbit2].setOrbit(orbit2) ;	// in the contained AttributeMultiVectors

	m_embeddings[orbit1]->swap(m_embeddings[orbit2]) ;


Sylvain Thery's avatar
Sylvain Thery committed
216
	// NE MARCHE PLUS PAS POSSIBLE DE CHANGER LES CONTAINER
untereiner's avatar
untereiner committed
217

Sylvain Thery's avatar
Sylvain Thery committed
218 219 220 221 222 223 224 225 226 227
//	for (unsigned int i = 0; i < NB_THREAD; ++i)
//	{
//		for(std::vector<CellMarkerGen*>::iterator it = cellMarkers[i].begin(); it != cellMarkers[i].end(); ++it)
//		{
//			if((*it)->m_cell == orbit1)
//				(*it)->m_cell = orbit2 ;
//			else if((*it)->m_cell == orbit2)
//				(*it)->m_cell = orbit1 ;
//		}
//	}
untereiner's avatar
untereiner committed
228
}
Pierre Kraemer's avatar
Pierre Kraemer committed
229

Pierre Kraemer's avatar
Pierre Kraemer committed
230 231 232 233 234 235 236 237 238 239 240 241
void GenericMap::viewAttributesTables()
{
	std::cout << "======================="<< std::endl ;
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
	{
		std::cout << "ATTRIBUTE_CONTAINER " << i << std::endl ;
		AttributeContainer& cont = m_attribs[i] ;

		// get the list of attributes
		std::vector<std::string> listeNames ;
		cont.getAttributesNames(listeNames) ;
		for (std::vector<std::string>::iterator it = listeNames.begin(); it != listeNames.end(); ++it)
242 243 244 245 246
		{
			unsigned int id = cont.getAttributeIndex(*it);
			std::cout << "    " << *it << " ("<<id<<")"<<std::endl ;
		}

Pierre Kraemer's avatar
Pierre Kraemer committed
247 248 249 250 251 252 253
		std::cout << "-------------------------" << std::endl ;
	}
	std::cout << "m_embeddings: " << std::hex ;
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
		std::cout << (long)(m_embeddings[i]) << " / " ;
	std::cout << std::endl << "-------------------------" << std::endl ;

Sylvain Thery's avatar
Sylvain Thery committed
254 255 256 257
//	std::cout << "m_markTables: " ;
//	for (unsigned int i = 0; i < NB_ORBITS; ++i)
//		std::cout << (long)(m_markTables[i][0]) << " / " ;
//	std::cout << std::endl << "-------------------------" << std::endl << std::dec ;
Pierre Kraemer's avatar
Pierre Kraemer committed
258 259
}

260 261 262 263 264 265 266
void GenericMap::printDartsTable()
{
	std::cout << "======================="<< std::endl ;

	//m_attribs[DART]
}

267 268 269 270
/****************************************
 *          THREAD MANAGEMENT           *
 ****************************************/

Sylvain Thery's avatar
Sylvain Thery committed
271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312
//void GenericMap::addThreadMarker(unsigned int nb)
//{
//	unsigned int th ;

//	for (unsigned int j = 0; j < nb; ++j)
//	{
//		th = m_nbThreadMarkers ;
//		m_nbThreadMarkers++ ;

//		for (unsigned int i = 0; i < NB_ORBITS; ++i)
//		{
//			std::stringstream ss ;
//			ss << "Mark_"<< th ;
//			AttributeContainer& cellCont = m_attribs[i] ;
//			AttributeMultiVector<Mark>* amvMark = cellCont.addAttribute<Mark>(ss.str()) ;
//			m_markTables[i][th] = amvMark ;
//		}
//	}
//}

//unsigned int GenericMap::getNbThreadMarkers() const
//{
//	return m_nbThreadMarkers;
//}

//void GenericMap::removeThreadMarker(unsigned int nb)
//{
//	unsigned int th = 0;
//	while ((m_nbThreadMarkers > 1) && (nb > 0))
//	{
//		th = --m_nbThreadMarkers ;
//		--nb;
//		for (unsigned int i = 0; i < NB_ORBITS; ++i)
//		{
//			std::stringstream ss ;
//			ss << "Mark_"<< th ;
//			AttributeContainer& cellCont = m_attribs[i] ;
//			cellCont.removeAttribute<Mark>(ss.str()) ;
//			m_markTables[i][th] = NULL ;
//		}
//	}
//}
Pierre Kraemer's avatar
Pierre Kraemer committed
313 314 315 316 317

/****************************************
 *             SAVE & LOAD              *
 ****************************************/

Pierre Kraemer's avatar
Pierre Kraemer committed
318
void GenericMap::restore_shortcuts()
Sylvain Thery's avatar
Sylvain Thery committed
319
{
Pierre Kraemer's avatar
Pierre Kraemer committed
320 321
	// EMBEDDING

Sylvain Thery's avatar
Sylvain Thery committed
322 323 324 325 326 327 328 329 330 331
	// get container of dart orbit
	AttributeContainer& cont = m_attribs[DART] ;

	// get the list of attributes
	std::vector<std::string> listeNames;
	cont.getAttributesNames(listeNames);

	// check if there are EMB_X attributes
	for (unsigned int i = 0;  i < listeNames.size(); ++i)
	{
Pierre Kraemer's avatar
merges  
Pierre Kraemer committed
332
		std::string sub = listeNames[i].substr(0, listeNames[i].size() - 1);
Sylvain Thery's avatar
Sylvain Thery committed
333 334 335 336 337 338 339 340
		if (sub == "EMB_")
		{
			unsigned int orb = listeNames[i][4]-'0'; // easy atoi computation for one char;
			AttributeMultiVector<unsigned int>* amv = cont.getDataVector<unsigned int>(i);
			m_embeddings[orb] = amv ;
		}
	}

Sylvain Thery's avatar
Sylvain Thery committed
341 342
	// MARKERS
	m_attribs[DART].getAttributesNames(listeNames);
Pierre Kraemer's avatar
Pierre Kraemer committed
343

Sylvain Thery's avatar
Sylvain Thery committed
344
	for (unsigned int i = 0;  i < listeNames.size(); ++i)
Sylvain Thery's avatar
Sylvain Thery committed
345
	{
Sylvain Thery's avatar
Sylvain Thery committed
346 347 348 349 350 351
		if (listeNames[i] == "BoundaryMark0")
			m_boundaryMarkers[0] = cont.getDataVector<MarkerBool>(i);

		if (listeNames[i] == "BoundaryMark1")
			m_boundaryMarkers[1] = cont.getDataVector<MarkerBool>(i);
	}
Sylvain Thery's avatar
Sylvain Thery committed
352

Sylvain Thery's avatar
Sylvain Thery committed
353
	// QUICK TRAVERSAL
Pierre Kraemer's avatar
Pierre Kraemer committed
354

Sylvain Thery's avatar
Sylvain Thery committed
355 356 357
	for(unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
	{
		AttributeContainer& cont = m_attribs[orbit];
Pierre Kraemer's avatar
Pierre Kraemer committed
358 359 360 361 362 363 364
		m_quickTraversal[orbit] = cont.getDataVector<Dart>("quick_traversal") ;
		for(unsigned int j = 0; j < NB_ORBITS; ++j)
		{
			std::stringstream ss;
			ss << "quickLocalIncidentTraversal_" << j;
			m_quickLocalIncidentTraversal[orbit][j] = cont.getDataVector< NoTypeNameAttribute<std::vector<Dart> > >(ss.str()) ;
			std::stringstream ss2;
Pierre Kraemer's avatar
Pierre Kraemer committed
365
			ss2 << "quickLocalAdjacentTraversal_" << j;
Pierre Kraemer's avatar
Pierre Kraemer committed
366 367
			m_quickLocalAdjacentTraversal[orbit][j] = cont.getDataVector< NoTypeNameAttribute<std::vector<Dart> > >(ss2.str()) ;
		}
368
	}
369

370 371 372 373
	// set Attribute handlers invalid
	for(std::multimap<AttributeMultiVectorGen*, AttributeHandlerGen*>::iterator it = attributeHandlers.begin(); it != attributeHandlers.end(); ++it)
		(*it).second->setInvalid() ;
	attributeHandlers.clear() ;
Sylvain Thery's avatar
Sylvain Thery committed
374 375 376 377
}

void GenericMap::dumpAttributesAndMarkers()
{
Pierre Kraemer's avatar
merges  
Pierre Kraemer committed
378
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
Sylvain Thery's avatar
Sylvain Thery committed
379 380
	{
		std::vector<std::string> names;
Pierre Kraemer's avatar
merges  
Pierre Kraemer committed
381
		names.reserve(32); 				//just to limit reallocation
Sylvain Thery's avatar
Sylvain Thery committed
382 383
		m_attribs[i].getAttributesNames(names);
		unsigned int nb = names.size();
Pierre Kraemer's avatar
merges  
Pierre Kraemer committed
384
		if (nb > 0)
Sylvain Thery's avatar
Sylvain Thery committed
385
		{
Pierre Kraemer's avatar
merges  
Pierre Kraemer committed
386
			CGoGNout << "ORBIT "<< i << CGoGNendl;
Sylvain Thery's avatar
Sylvain Thery committed
387 388 389
			std::vector<std::string> types;
			types.reserve(nb);
			m_attribs[i].getAttributesTypes(types);
Pierre Kraemer's avatar
merges  
Pierre Kraemer committed
390 391
			for (unsigned int j = 0; j < nb; ++j)
				CGoGNout << "    " << j << " : " << types[j] << " " << names[j] << CGoGNendl;
Sylvain Thery's avatar
Sylvain Thery committed
392 393
		}
	}
Sylvain Thery's avatar
Sylvain Thery committed
394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413
//	CGoGNout << "RESERVED MARKERS "<< CGoGNendl;
//	for (unsigned int i = 0; i < NB_ORBITS; ++i)
//	{
//		for (unsigned int j = 0; j < NB_THREAD; ++j)
//		{
//			MarkSet ms = m_marksets[i][j];
//			if (!ms.isClear())
//			{
//				CGoGNout << "Orbit " << i << "  thread " << j << " : ";
//				Mark m(1);
//				for (unsigned i = 0; i < Mark::getNbMarks(); ++i)
//				{
//					if (ms.testMark(m))
//						CGoGNout << m.getMarkVal() << ", ";
//					m.setMarkVal(m.getMarkVal()<<1);
//				}
//				CGoGNout << CGoGNendl;
//			}
//		}
//	}
Sylvain Thery's avatar
Sylvain Thery committed
414 415
}

416 417
void GenericMap::compact()
{
418 419 420
	// compact embedding attribs
	std::vector< std::vector<unsigned int>* > oldnews;
	oldnews.resize(NB_ORBITS);
421 422 423 424
	for (unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
	{
		if ((orbit != DART) && (isOrbitEmbedded(orbit)))
		{
425 426 427 428
			oldnews[orbit] = new std::vector<unsigned int>;
			m_attribs[orbit].compact(*(oldnews[orbit]));
		}
	}
429

430 431 432 433 434 435
	// update embedding indices of topo
	for (unsigned int i = m_attribs[DART].begin(); i != m_attribs[DART].end(); m_attribs[DART].next(i))
	{
		for (unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
		{
			if ((orbit != DART) && (isOrbitEmbedded(orbit)))
436
			{
437 438
				unsigned int& idx = m_embeddings[orbit]->operator[](i);
				unsigned int jdx = oldnews[orbit]->operator[](idx);
Pierre Kraemer's avatar
merges  
Pierre Kraemer committed
439
				if ((jdx != 0xffffffff) && (jdx != idx))
440 441 442 443
					idx = jdx;
			}
		}
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
444

445 446 447 448 449
	// delete allocated vectors
	for (unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
		if ((orbit != DART) && (isOrbitEmbedded(orbit)))
			delete[] oldnews[orbit];

450 451
	// compact topo (depends on map implementation)
	compactTopo();
452 453
}

454 455 456 457 458 459 460 461 462
void GenericMap::dumpCSV() const
{
	for (unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
	{
		m_attribs[orbit].dumpCSV();
	}
	CGoGNout << CGoGNendl;
}

Pierre Kraemer's avatar
Pierre Kraemer committed
463
} // namespace CGoGN