genericmap.cpp 31 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"
Sylvain Thery's avatar
Sylvain Thery committed
27
28
#include "Topology/generic/traversorCell.h"

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
33
34

namespace CGoGN
{
35

Pierre Kraemer's avatar
Pierre Kraemer committed
36
std::map<std::string, RegisteredBaseAttribute*>* GenericMap::m_attributes_registry_map = NULL;
37
int GenericMap::m_nbInstances = 0;
38

39
GenericMap::GenericMap() : m_nbThreads(1)
Pierre Kraemer's avatar
Pierre Kraemer committed
40
{
41
	if(m_attributes_registry_map == NULL)
Pierre Kraemer's avatar
Pierre Kraemer committed
42
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
43
		m_attributes_registry_map = new std::map<std::string, RegisteredBaseAttribute*>;
44

Pierre Kraemer's avatar
Pierre Kraemer committed
45
46
47
		// register all known types
		registerAttribute<Dart>("Dart");
		registerAttribute<Mark>("Mark");
Pierre Kraemer's avatar
Pierre Kraemer committed
48

Pierre Kraemer's avatar
Pierre Kraemer committed
49
50
51
52
		registerAttribute<long>("long");
		registerAttribute<int>("int");
		registerAttribute<short>("short");
		registerAttribute<char>("char");
Pierre Kraemer's avatar
Pierre Kraemer committed
53

Pierre Kraemer's avatar
Pierre Kraemer committed
54
55
56
57
		registerAttribute<unsigned long>("unsigned long");
		registerAttribute<unsigned int>("unsigned int");
		registerAttribute<unsigned short>("unsigned short");
		registerAttribute<unsigned char>("unsigned char");
Pierre Kraemer's avatar
Pierre Kraemer committed
58

Pierre Kraemer's avatar
Pierre Kraemer committed
59
60
61
		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
62

Pierre Kraemer's avatar
Pierre Kraemer committed
63
64
65
		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
66

Pierre Kraemer's avatar
Pierre Kraemer committed
67
68
		registerAttribute<Geom::Matrix33f>(Geom::Matrix33f::CGoGNnameOfType());
		registerAttribute<Geom::Matrix44f>(Geom::Matrix44f::CGoGNnameOfType());
Pierre Kraemer's avatar
Pierre Kraemer committed
69

Pierre Kraemer's avatar
Pierre Kraemer committed
70
71
72
73
74
		registerAttribute<Geom::Matrix33d>(Geom::Matrix33d::CGoGNnameOfType());
		registerAttribute<Geom::Matrix44d>(Geom::Matrix44d::CGoGNnameOfType());
	}

	m_nbInstances++;
Pierre Kraemer's avatar
Pierre Kraemer committed
75

76
	for(unsigned int i = 0; i < NB_ORBITS; ++i)
Pierre Kraemer's avatar
Pierre Kraemer committed
77
	{
78
		m_attribs[i].setOrbit(i) ;
79
		m_attribs[i].setRegistry(m_attributes_registry_map) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
80
		m_embeddings[i] = NULL ;
Pierre Kraemer's avatar
Pierre Kraemer committed
81
		m_quickTraversal[i] = NULL ;
82
83
84
		for(unsigned int j = 0; j < NB_THREAD; ++j)
		{
			m_marksets[i][j].clear() ;
Pierre Kraemer's avatar
Pierre Kraemer committed
85
			m_markTables[i][j] = NULL ;
86
		}
Pierre Kraemer's avatar
Pierre Kraemer committed
87
	}
88

89
90
91
92
93
	for (unsigned int i=0; i<NB_THREAD; ++i)
	{
		dartMarkers[i].reserve(16) ;
		cellMarkers[i].reserve(16) ;
	}
94

95
	// get & lock marker for boundary
Thery Sylvain's avatar
Thery Sylvain committed
96
97
	m_boundaryMarkers[0] =  m_marksets[DART][0].getNewMark();
	m_boundaryMarkers[1] =  m_marksets[DART][0].getNewMark();
98

untereiner's avatar
untereiner committed
99
100
	if (m_isMultiRes)
		initMR() ;
Pierre Kraemer's avatar
Pierre Kraemer committed
101
102
103
104
}

GenericMap::~GenericMap()
{
105
	// release marker for boundary
Thery Sylvain's avatar
Thery Sylvain committed
106
107
	m_marksets[DART][0].releaseMark(m_boundaryMarkers[0]);
	m_marksets[DART][0].releaseMark(m_boundaryMarkers[1]);
108

Pierre Kraemer's avatar
Pierre Kraemer committed
109
110
111
112
113
	for(unsigned int i = 0; i < NB_ORBITS; ++i)
	{
		if(isOrbitEmbedded(i))
			m_attribs[i].clear(true) ;
	}
114
115
116
117
118

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

119
120
121
122
123
	for (unsigned int i=0; i<NB_THREAD; ++i)
	{
		for(std::vector<DartMarkerGen*>::iterator it = dartMarkers[i].begin(); it != dartMarkers[i].end(); ++it)
			(*it)->setReleaseOnDestruct(false) ;
		dartMarkers[i].clear() ;
124

125
126
127
128
		for(std::vector<CellMarkerGen*>::iterator it = cellMarkers[i].begin(); it != cellMarkers[i].end(); ++it)
			(*it)->setReleaseOnDestruct(false) ;
		cellMarkers[i].clear() ;
	}
129
130
131

	// clean type registry if necessary
	m_nbInstances--;
Pierre Kraemer's avatar
Pierre Kraemer committed
132
	if (m_nbInstances <= 0)
133
	{
134
135
136
		for (std::map<std::string, RegisteredBaseAttribute*>::iterator it =  m_attributes_registry_map->begin(); it != m_attributes_registry_map->end(); ++it)
			delete it->second;

137
		delete m_attributes_registry_map;
138
139
		m_attributes_registry_map = NULL;
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
140
141
}

142
143
void GenericMap::clear(bool removeAttrib)
{
Sylvain Thery's avatar
Sylvain Thery committed
144
	if (removeAttrib)
145
	{
Sylvain Thery's avatar
Sylvain Thery committed
146
		for(unsigned int i = 0; i < NB_ORBITS; ++i)
147
		{
Sylvain Thery's avatar
Sylvain Thery committed
148
			m_attribs[i].clear(true) ;
149
			m_embeddings[i] = NULL ;
150
			m_quickTraversal[i] = NULL;
151
		}
152
153
		for(std::multimap<AttributeMultiVectorGen*, AttributeHandlerGen*>::iterator it = attributeHandlers.begin(); it != attributeHandlers.end(); ++it)
			(*it).second->setInvalid() ;
154
		attributeHandlers.clear() ;
155
	}
Sylvain Thery's avatar
Sylvain Thery committed
156
157
158
159
160
	else
	{
		for(unsigned int i = 0; i < NB_ORBITS; ++i)
			m_attribs[i].clear(false) ;
	}
161
162

	if (m_isMultiRes)
163
		initMR() ;
164
165
}

Pierre Kraemer's avatar
Pierre Kraemer committed
166
167
168
169
/****************************************
 *           MULTIRES                   *
 ****************************************/

Pierre Kraemer's avatar
Pierre Kraemer committed
170
171
172
173
174
175
176
177
178
179
void GenericMap::printMR()
{
	std::cout << std::endl ;

	for(unsigned int j = 0; j < m_mrNbDarts.size(); ++j)
		std::cout << m_mrNbDarts[j] << " / " ;
	std::cout << std::endl << "==========" << std::endl ;

	for(unsigned int i = m_mrattribs.begin(); i != m_mrattribs.end(); m_mrattribs.next(i))
	{
180
		std::cout << i << " : " << (*m_mrLevels)[i] << " / " ;
Pierre Kraemer's avatar
Pierre Kraemer committed
181
182
183
184
185
186
187
188
189
		for(unsigned int j = 0; j < m_mrDarts.size(); ++j)
			std::cout << (*m_mrDarts[j])[i] << " ; " ;
		std::cout << std::endl ;
	}
}

void GenericMap::initMR()
{
	m_mrattribs.clear(true) ;
Sylvain Thery's avatar
Sylvain Thery committed
190
191
	m_mrattribs.setRegistry(m_attributes_registry_map) ;

Pierre Kraemer's avatar
Pierre Kraemer committed
192
193
194
195
196
197
198
199
	m_mrDarts.clear() ;
	m_mrDarts.reserve(16) ;
	m_mrNbDarts.clear();
	m_mrNbDarts.reserve(16);
	m_mrLevelStack.clear() ;
	m_mrLevelStack.reserve(16) ;

	m_mrLevels = m_mrattribs.addAttribute<unsigned int>("MRLevel") ;
untereiner's avatar
untereiner committed
200
201
202
203
204

	AttributeMultiVector<unsigned int>* newAttrib = m_mrattribs.addAttribute<unsigned int>("MRdart_0") ;
	m_mrDarts.push_back(newAttrib) ;
	m_mrNbDarts.push_back(0) ;

Pierre Kraemer's avatar
Pierre Kraemer committed
205
206
207
	setCurrentLevel(0) ;
}

208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
//AttributeMultiVector<unsigned int>* GenericMap::addLevel()
//{
//	unsigned int newLevel = m_mrDarts.size() ;
//	std::stringstream ss ;
//	ss << "MRdart_"<< newLevel ;
//	AttributeMultiVector<unsigned int>* newAttrib = m_mrattribs.addAttribute<unsigned int>(ss.str()) ;
//	AttributeMultiVector<unsigned int>* prevAttrib = m_mrDarts[newLevel - 1];
//
//	// copy the indices of previous level into new level
//	m_mrattribs.copyAttribute(newAttrib->getIndex(), prevAttrib->getIndex()) ;
//
//	return newAttrib;
//}

void GenericMap::addLevelBack()
Pierre Kraemer's avatar
Pierre Kraemer committed
223
{
224
225
	//AttributeMultiVector<unsigned int>* newAttrib = addLevel();

Pierre Kraemer's avatar
Pierre Kraemer committed
226
	unsigned int newLevel = m_mrDarts.size() ;
Pierre Kraemer's avatar
Pierre Kraemer committed
227
	std::stringstream ss ;
Pierre Kraemer's avatar
Pierre Kraemer committed
228
	ss << "MRdart_"<< newLevel ;
229
	AttributeMultiVector<unsigned int>* newAttrib = m_mrattribs.addAttribute<unsigned int>(ss.str()) ;
untereiner's avatar
untereiner committed
230
	AttributeMultiVector<unsigned int>* prevAttrib = m_mrDarts[newLevel - 1];
Pierre Kraemer's avatar
Pierre Kraemer committed
231

untereiner's avatar
untereiner committed
232
233
	// copy the indices of previous level into new level
	m_mrattribs.copyAttribute(newAttrib->getIndex(), prevAttrib->getIndex()) ;
234

untereiner's avatar
untereiner committed
235
236
237
	m_mrDarts.push_back(newAttrib) ;
	m_mrNbDarts.push_back(0) ;
}
238

untereiner's avatar
untereiner committed
239
240
void GenericMap::addLevelFront()
{
241
242
243
244
245
246
247
248
249
250
251
	//AttributeMultiVector<unsigned int>* newAttrib = addLevel();

	unsigned int newLevel = m_mrDarts.size() ;
	std::stringstream ss ;
	ss << "MRdart_"<< newLevel ;
	AttributeMultiVector<unsigned int>* newAttrib = m_mrattribs.addAttribute<unsigned int>(ss.str()) ;
	AttributeMultiVector<unsigned int>* prevAttrib = m_mrDarts[0];

	// copy the indices of previous level into new level
	m_mrattribs.copyAttribute(newAttrib->getIndex(), prevAttrib->getIndex()) ;

252
253
254
255
	m_mrDarts.insert(m_mrDarts.begin(), newAttrib) ;
	m_mrNbDarts.insert(m_mrNbDarts.begin(), 0) ;
}

untereiner's avatar
untereiner committed
256
void GenericMap::removeLevelBack()
Pierre Kraemer's avatar
Pierre Kraemer committed
257
{
Pierre Kraemer's avatar
Pierre Kraemer committed
258
259
	unsigned int maxL = getMaxLevel() ;
	if(maxL > 0)
Pierre Kraemer's avatar
Pierre Kraemer committed
260
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
261
262
		AttributeMultiVector<unsigned int>* maxMR = m_mrDarts[maxL] ;
		AttributeMultiVector<unsigned int>* prevMR = m_mrDarts[maxL - 1] ;
Pierre Kraemer's avatar
Pierre Kraemer committed
263
264
		for(unsigned int i = m_mrattribs.begin(); i != m_mrattribs.end(); m_mrattribs.next(i))
		{
Pierre Kraemer's avatar
Pierre Kraemer committed
265
266
267
268
269
270
271
272
273
274
275
			unsigned int idx = (*maxMR)[i] ;
			if((*m_mrLevels)[i] == maxL)	// if the MRdart was introduced on the level we're removing
			{
				deleteDartLine(idx) ;		// delete the pointed dart line
				m_mrattribs.removeLine(i) ;	// delete the MRdart line
			}
			else							// if the dart was introduced on a previous level
			{
				if(idx != (*prevMR)[i])		// delete the pointed dart line only if
					deleteDartLine(idx) ;	// it is not shared with previous level
			}
Pierre Kraemer's avatar
Pierre Kraemer committed
276
277
		}

Pierre Kraemer's avatar
Pierre Kraemer committed
278
		m_mrattribs.removeAttribute<unsigned int>(maxMR->getIndex()) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
279
280
281
282
283
		m_mrDarts.pop_back() ;
		m_mrNbDarts.pop_back() ;

		if(m_mrCurrentLevel == maxL)
			--m_mrCurrentLevel ;
Pierre Kraemer's avatar
Pierre Kraemer committed
284
285
286
	}
}

untereiner's avatar
untereiner committed
287
void GenericMap::removeLevelFront()
288
{
289
290
291
292
	unsigned int maxL = getMaxLevel() ;
	if(maxL > 0) //must have at min 2 levels (0 and 1) to remove the front one
	{
		AttributeMultiVector<unsigned int>* minMR = m_mrDarts[0] ;
Thery Sylvain's avatar
Thery Sylvain committed
293
//		AttributeMultiVector<unsigned int>* firstMR = m_mrDarts[1] ;
294
295
		for(unsigned int i = m_mrattribs.begin(); i != m_mrattribs.end(); m_mrattribs.next(i))
		{
Thery Sylvain's avatar
Thery Sylvain committed
296
//			unsigned int idx = (*minMR)[i] ;
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
			if((*m_mrLevels)[i] != 0)	// if the MRdart was introduced after the level we're removing
			{
				--(*m_mrLevels)[i]; //decrement his level of insertion
			}
			else							// if the dart was introduced on a this level and not used after
			{
//				if(idx != (*firstMR)[i])		// delete the pointed dart line only if
//					deleteDartLine(idx) ;	// it is not shared with next level
			}
		}

		m_mrNbDarts[1] += m_mrNbDarts[0];

		m_mrattribs.removeAttribute<unsigned int>(minMR->getIndex()) ;
		m_mrDarts.erase(m_mrDarts.begin()) ;
		m_mrNbDarts.erase(m_mrNbDarts.begin()) ;

		--m_mrCurrentLevel ;
	}
316
317
}

untereiner's avatar
untereiner committed
318
319
320
321
322
323
324
325
326
327
328
void GenericMap::copyLevel(unsigned int level)
{
	AttributeMultiVector<unsigned int>* newAttrib = m_mrDarts[level] ;
	AttributeMultiVector<unsigned int>* prevAttrib = m_mrDarts[level - 1];

	// copy the indices of previous level into new level
	m_mrattribs.copyAttribute(newAttrib->getIndex(), prevAttrib->getIndex()) ;
}

void GenericMap::duplicateDarts(unsigned int level)
{
329
	AttributeMultiVector<unsigned int>* attrib = m_mrDarts[level] ;
untereiner's avatar
untereiner committed
330
331
332

	for(unsigned int i = m_mrattribs.begin(); i != m_mrattribs.end(); m_mrattribs.next(i))
	{
333
334
		unsigned int oldi = (*attrib)[i] ;	// get the index of the dart in previous level
		(*attrib)[i] = copyDartLine(oldi) ;	// copy the dart and affect it to the new level
untereiner's avatar
untereiner committed
335
336
337
	}
}

Pierre Kraemer's avatar
Pierre Kraemer committed
338
339
340
341
/****************************************
 *        ATTRIBUTES MANAGEMENT         *
 ****************************************/

untereiner's avatar
untereiner committed
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
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]) ;

	for(unsigned int t = 0; t < m_nbThreads; ++t)
	{
		AttributeMultiVector<Mark>* m = m_markTables[orbit1][t] ;
		m_markTables[orbit1][t] = m_markTables[orbit2][t] ;
		m_markTables[orbit2][t] = m ;

		MarkSet ms = m_marksets[orbit1][t] ;
		m_marksets[orbit1][t] = m_marksets[orbit2][t] ;
		m_marksets[orbit2][t] = ms ;
	}

	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 ;
		}
	}
}
Pierre Kraemer's avatar
Pierre Kraemer committed
375

Pierre Kraemer's avatar
Pierre Kraemer committed
376
377
378
379
380
381
382
383
384
385
386
387
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)
388
389
390
391
392
		{
			unsigned int id = cont.getAttributeIndex(*it);
			std::cout << "    " << *it << " ("<<id<<")"<<std::endl ;
		}

Pierre Kraemer's avatar
Pierre Kraemer committed
393
394
395
396
397
398
399
400
401
402
403
404
405
		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 ;

	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 ;
}

406
407
408
409
410
411
412
void GenericMap::printDartsTable()
{
	std::cout << "======================="<< std::endl ;

	//m_attribs[DART]
}

413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
/****************************************
 *          THREAD MANAGEMENT           *
 ****************************************/

void GenericMap::addThreadMarker(unsigned int nb)
{
	unsigned int th ;

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

		for (unsigned int i = 0; i < NB_ORBITS; ++i)
		{
428
429
430
431
432
			std::stringstream ss ;
			ss << "Mark_"<< th ;
			AttributeContainer& cellCont = m_attribs[i] ;
			AttributeMultiVector<Mark>* amvMark = cellCont.addAttribute<Mark>(ss.str()) ;
			m_markTables[i][th] = amvMark ;
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
		}
	}
}

unsigned int GenericMap::getNbThreadMarkers()
{
	return m_nbThreads;
}

void GenericMap::removeThreadMarker(unsigned int nb)
{
	unsigned int th = 0;
	while ((m_nbThreads > 1) && (nb > 0))
	{
		th = --m_nbThreads ;
		--nb;
		for (unsigned int i = 0; i < NB_ORBITS; ++i)
		{
451
452
453
454
455
			std::stringstream ss ;
			ss << "Mark_"<< th ;
			AttributeContainer& cellCont = m_attribs[i] ;
			cellCont.removeAttribute<Mark>(ss.str()) ;
			m_markTables[i][th] = NULL ;
456
457
458
		}
	}
}
Pierre Kraemer's avatar
Pierre Kraemer committed
459
460
461
462
463
464
465
466
467
468

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

bool GenericMap::saveMapBin(const std::string& filename)
{
	CGoGNostream fs(filename.c_str(), std::ios::out|std::ios::binary);
	if (!fs)
	{
469
		CGoGNerr << "Unable to open file for writing: " << filename << CGoGNendl;
Pierre Kraemer's avatar
Pierre Kraemer committed
470
471
472
473
474
		return false;
	}

	// Entete
	char* buff = new char[256];
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
475
	for (int i = 0; i < 256; ++i)
Pierre Kraemer's avatar
Pierre Kraemer committed
476
		buff[i] = char(255);
477
478
479
480
481
482
483
484
485
486
	if (m_isMultiRes)
	{
		const char* cgogn = "CGoGN_MRMap";
		memcpy(buff, cgogn, 12);
	}
	else
	{
		const char* cgogn = "CGoGN_Map";
		memcpy(buff, cgogn, 10);
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
487
488
	std::string mt = mapTypeName();
	const char* mtc = mt.c_str();
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
489
490
	memcpy(buff+32, mtc, mt.size()+1);
	unsigned int *buffi = reinterpret_cast<unsigned int*>(buff + 64);
Pierre Kraemer's avatar
Pierre Kraemer committed
491
492
493
494
	*buffi = NB_ORBITS;
	fs.write(reinterpret_cast<const char*>(buff), 256);
	delete buff;

Sylvain Thery's avatar
Sylvain Thery committed
495
	// save all attribs
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
496
497
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
		m_attribs[i].saveBin(fs, i);
Pierre Kraemer's avatar
Pierre Kraemer committed
498

499
	if (m_isMultiRes)
Sylvain Thery's avatar
Sylvain Thery committed
500
	{
501
502
		m_mrattribs.saveBin(fs, 00);

Sylvain Thery's avatar
Sylvain Thery committed
503
504
505
506
507
508
509
		fs.write(reinterpret_cast<const char*>(&m_mrCurrentLevel), sizeof(unsigned int));

		unsigned int nb = m_mrNbDarts.size();
		fs.write(reinterpret_cast<const char*>(&nb), sizeof(unsigned int));
		fs.write(reinterpret_cast<const char*>(&(m_mrNbDarts[0])), nb *sizeof(unsigned int));
	}

Pierre Kraemer's avatar
Pierre Kraemer committed
510
511
512
513
514
	return true;
}

bool GenericMap::loadMapBin(const std::string& filename)
{
Sylvain Thery's avatar
Sylvain Thery committed
515

Pierre Kraemer's avatar
Pierre Kraemer committed
516
517
518
	CGoGNistream fs(filename.c_str(), std::ios::in|std::ios::binary);
	if (!fs)
	{
519
		CGoGNerr << "Unable to open file for loading" << CGoGNendl;
Pierre Kraemer's avatar
Pierre Kraemer committed
520
521
522
		return false;
	}

Sylvain Thery's avatar
Sylvain Thery committed
523
524
525
	GenericMap::clear(true);

	// read info
Pierre Kraemer's avatar
Pierre Kraemer committed
526
527
528
529
530
	char* buff = new char[256];
	fs.read(reinterpret_cast<char*>(buff), 256);

	std::string buff_str(buff);
	// Check file type
531
	if (m_isMultiRes)
Pierre Kraemer's avatar
Pierre Kraemer committed
532
	{
533
534
535
536
537
538
539
540
541
542
		if (buff_str == "CGoGN_Map")
		{
			CGoGNerr<< "Wrong binary file format, file is not a MR-Map"<< CGoGNendl;
			return false;
		}
		if (buff_str != "CGoGN_MRMap")
		{
			CGoGNerr<< "Wrong binary file format"<< CGoGNendl;
			return false;
		}
Pierre Kraemer's avatar
Pierre Kraemer committed
543
	}
544
545
546
547
548
549
550
551
552
553
554
555
556
557
	else
	{
		if (buff_str == "CGoGN_MRMap")
		{
			CGoGNerr<< "Wrong binary file format, file is a MR-Map"<< CGoGNendl;
			return false;
		}
		if (buff_str != "CGoGN_Map")
		{
			CGoGNerr<< "Wrong binary file format"<< CGoGNendl;
			return false;
		}
	}

Pierre Kraemer's avatar
Pierre Kraemer committed
558
559

	// Check map type
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
560
	buff_str = std::string(buff + 32);
Pierre Kraemer's avatar
Pierre Kraemer committed
561
562

	std::string localType = this->mapTypeName();
Sylvain Thery's avatar
Sylvain Thery committed
563
564

	std::string fileType = buff_str;
Pierre Kraemer's avatar
Pierre Kraemer committed
565
566
567

	if (fileType != localType)
	{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
568
		CGoGNerr << "Not possible to load "<< fileType << " into " << localType << " object" << CGoGNendl;
Pierre Kraemer's avatar
Pierre Kraemer committed
569
570
571
572
		return false;
	}

	// Check max nb orbit
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
573
	unsigned int *ptr_nbo = reinterpret_cast<unsigned int*>(buff + 64);
Pierre Kraemer's avatar
Pierre Kraemer committed
574
575
576
	unsigned int nbo = *ptr_nbo;
	if (nbo != NB_ORBITS)
	{
577
		CGoGNerr << "Wrond max orbit number in file" << CGoGNendl;
Pierre Kraemer's avatar
Pierre Kraemer committed
578
579
580
581
		return  false;
	}

	// load attrib container
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
582
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
Pierre Kraemer's avatar
Pierre Kraemer committed
583
	{
584
		unsigned int id = AttributeContainer::loadBinId(fs);
Pierre Kraemer's avatar
Pierre Kraemer committed
585
586
587
		m_attribs[id].loadBin(fs);
	}

588
	if (m_isMultiRes)
Sylvain Thery's avatar
Sylvain Thery committed
589
590
	{
		AttributeContainer::loadBinId(fs); // not used but need to read to skip
591
592
		m_mrattribs.loadBin(fs);

Sylvain Thery's avatar
Sylvain Thery committed
593
594
595
596
597
598
599
		fs.read(reinterpret_cast<char*>(&m_mrCurrentLevel), sizeof(unsigned int));
		unsigned int nb;
		fs.read(reinterpret_cast<char*>(&nb), sizeof(unsigned int));
		m_mrNbDarts.resize(nb);
		fs.read(reinterpret_cast<char*>(&(m_mrNbDarts[0])), nb *sizeof(unsigned int));
	}

600
	// retrieve m_embeddings (from m_attribs)
Pierre Kraemer's avatar
Pierre Kraemer committed
601
	update_m_emb_afterLoad();
602

Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
603
	// recursive call from real type of map (for topo relation attributes pointers) down to GenericMap (for Marker_cleaning & pointers)
Sylvain Thery's avatar
Sylvain Thery committed
604
	update_topo_shortcuts();
Sylvain Thery's avatar
Sylvain Thery committed
605

606
607
608
609
610
611
612
613
614
615
616
617
618
	// restore nbThreads
	std::vector<std::string> typeMark;
	unsigned int nbatt0 = m_attribs[0].getAttributesTypes(typeMark);
	m_nbThreads = 0;
	for (unsigned int i = 0; i < nbatt0; ++i)
	{
		if (typeMark[i] == "Mark")
			++m_nbThreads;
	}




Pierre Kraemer's avatar
Pierre Kraemer committed
619
620
621
	return true;
}

Sylvain Thery's avatar
Sylvain Thery committed
622
623
624
625
bool GenericMap::copyFrom(const GenericMap& map)
{
	if (mapTypeName() != map.mapTypeName())
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
626
		CGoGNerr << "try to copy from incompatible type map" << CGoGNendl;
Sylvain Thery's avatar
Sylvain Thery committed
627
628
629
630
631
632
633
634
635
636
637
638
639
640
		return false;
	}

	GenericMap::clear(true);

	// load attrib container
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
		m_attribs[i].copyFrom(map.m_attribs[i]);

	if (m_isMultiRes)
	{
		m_mrattribs.copyFrom(map.m_mrattribs);
		m_mrCurrentLevel = map.m_mrCurrentLevel;

Pierre Kraemer's avatar
Pierre Kraemer committed
641
		unsigned int nb = map.m_mrNbDarts.size();
Sylvain Thery's avatar
Sylvain Thery committed
642
		m_mrNbDarts.resize(nb);
Pierre Kraemer's avatar
Pierre Kraemer committed
643
		for (unsigned int i = 0; i < nb; ++i)
Sylvain Thery's avatar
Sylvain Thery committed
644
645
646
647
648
649
650
651
652
653
654
655
			m_mrNbDarts[i] = map.m_mrNbDarts[i];
	}

	// retrieve m_embeddings (from m_attribs)
	update_m_emb_afterLoad();

	// recursive call from real type of map (for topo relation attributes pointers) down to GenericMap (for Marker_cleaning & pointers)
	update_topo_shortcuts();

	return true;
}

Sylvain Thery's avatar
Sylvain Thery committed
656
657
658
659
660
661
662
663
664
665
666
667
void GenericMap::update_m_emb_afterLoad()
{
	// 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
668
		std::string sub = listeNames[i].substr(0, listeNames[i].size() - 1);
Sylvain Thery's avatar
Sylvain Thery committed
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
		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 ;
		}
	}
}

void GenericMap::update_topo_shortcuts()
{
	for(unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
	{
		AttributeContainer& cont = m_attribs[orbit];

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

		for (unsigned int i = 0;  i < listeNames.size(); ++i)
		{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
690
			std::string sub = listeNames[i].substr(0, 5);
Sylvain Thery's avatar
Sylvain Thery committed
691
692
693
694
			if (sub == "Mark_")
			{
				// get thread number
				unsigned int thread = listeNames[i][5]-'0';
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
695
				if (listeNames[i].size() > 6) 					// thread number is >9
Pierre Kraemer's avatar
Pierre Kraemer committed
696
					thread = 10*thread + (listeNames[i][6]-'0');
Sylvain Thery's avatar
Sylvain Thery committed
697
698
699
700

				AttributeMultiVector<Mark>* amvMark = cont.getDataVector<Mark>(i);
				m_markTables[orbit][thread] = amvMark ;

Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
701
				if ((orbit == DART) && (thread == 0))	// for Marker of dart of thread O keep the boundary marker
Sylvain Thery's avatar
Sylvain Thery committed
702
				{
Thery Sylvain's avatar
Thery Sylvain committed
703
704
705
			// TODO Verifier ce qu fait exactement ce unsetMark sur m.invert ??
//					Mark m(m_boundaryMarker);
					Mark m(m_boundaryMarkers[0]+m_boundaryMarkers[1]);
Sylvain Thery's avatar
Sylvain Thery committed
706
					m.invert();
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
707
					for (unsigned int i = cont.begin(); i != cont.end(); cont.next(i))
Sylvain Thery's avatar
Sylvain Thery committed
708
709
710
711
						amvMark->operator[](i).unsetMark(m);
				}
				else								// for others clear all
				{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
712
					for (unsigned int i = cont.begin(); i != cont.end(); cont.next(i))
Sylvain Thery's avatar
Sylvain Thery committed
713
714
715
716
717
						amvMark->operator[](i).clear();
				}
			}
		}
	}
718
719
720
721
722

	if (m_isMultiRes)
	{
		std::vector<std::string> names;
		m_mrattribs.getAttributesNames(names);
Pierre Kraemer's avatar
Pierre Kraemer committed
723
724
		m_mrDarts.resize(names.size() - 1);
		for (unsigned int i = 0; i < m_mrDarts.size(); ++i)
725
726
727
728
729
730
			m_mrDarts[i] = NULL;

		for (unsigned int i = 0;  i < names.size(); ++i)
		{
			std::string sub = names[i].substr(0, 7);

Pierre Kraemer's avatar
Pierre Kraemer committed
731
			if (sub == "MRLevel")
732
				m_mrLevels = m_mrattribs.getDataVector<unsigned int>(i);
733

Pierre Kraemer's avatar
Pierre Kraemer committed
734
			if (sub == "MRdart_")
735
736
			{
				sub = names[i].substr(7);	// compute number following MT_Dart_
Pierre Kraemer's avatar
Pierre Kraemer committed
737
738
739
740
				unsigned int idx = 0;
				for (unsigned int j = 0; j < sub.length(); j++)
					idx = 10*idx + (sub[j]-'0');
				if (idx < names.size() - 1)
741
					m_mrDarts[idx] = m_mrattribs.getDataVector<unsigned int>(i);
742
				else
Pierre Kraemer's avatar
Pierre Kraemer committed
743
					CGoGNerr << "Warning problem updating MR_DARTS" << CGoGNendl;
744
745
746
			}
		}
		// check if all pointers are != NULL
Pierre Kraemer's avatar
Pierre Kraemer committed
747
		for (unsigned int i = 0; i < m_mrDarts.size(); ++i)
748
749
		{
			if (m_mrDarts[i] == NULL)
Pierre Kraemer's avatar
Pierre Kraemer committed
750
				CGoGNerr << "Warning problem MR_DARTS = NULL" << CGoGNendl;
751
752
		}
	}
Sylvain Thery's avatar
Sylvain Thery committed
753
754
755
756
}

void GenericMap::dumpAttributesAndMarkers()
{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
757
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
Sylvain Thery's avatar
Sylvain Thery committed
758
759
	{
		std::vector<std::string> names;
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
760
		names.reserve(32); 				//just to limit reallocation
Sylvain Thery's avatar
Sylvain Thery committed
761
762
		m_attribs[i].getAttributesNames(names);
		unsigned int nb = names.size();
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
763
		if (nb > 0)
Sylvain Thery's avatar
Sylvain Thery committed
764
		{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
765
			CGoGNout << "ORBIT "<< i << CGoGNendl;
Sylvain Thery's avatar
Sylvain Thery committed
766
767
768
			std::vector<std::string> types;
			types.reserve(nb);
			m_attribs[i].getAttributesTypes(types);
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
769
770
			for (unsigned int j = 0; j < nb; ++j)
				CGoGNout << "    " << j << " : " << types[j] << " " << names[j] << CGoGNendl;
Sylvain Thery's avatar
Sylvain Thery committed
771
772
773
		}
	}
	CGoGNout << "RESERVED MARKERS "<< CGoGNendl;
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
774
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
Sylvain Thery's avatar
Sylvain Thery committed
775
	{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
776
		for (unsigned int j = 0; j < NB_THREAD; ++j)
Sylvain Thery's avatar
Sylvain Thery committed
777
778
779
780
		{
			MarkSet ms = m_marksets[i][j];
			if (!ms.isClear())
			{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
781
				CGoGNout << "Orbit " << i << "  thread " << j << " : ";
Sylvain Thery's avatar
Sylvain Thery committed
782
				Mark m(1);
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
783
				for (unsigned i = 0; i < Mark::getNbMarks(); ++i)
Sylvain Thery's avatar
Sylvain Thery committed
784
785
				{
					if (ms.testMark(m))
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
786
						CGoGNout << m.getMarkVal() << ", ";
Sylvain Thery's avatar
Sylvain Thery committed
787
788
789
790
791
792
793
794
					m.setMarkVal(m.getMarkVal()<<1);
				}
				CGoGNout << CGoGNendl;
			}
		}
	}
}

795
796
void GenericMap::compact()
{
797
798
799
800
	// if MR compact the MR attrib container
	std::vector<unsigned int> oldnewMR;
	if (m_isMultiRes)
		m_mrattribs.compact(oldnewMR);
801
802

	// compacting the orbits attributes
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
//	for (unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
//	{
//		if ((orbit != DART) && (isOrbitEmbedded(orbit)))
//		{
//			m_attribs[orbit].compact(oldnew);
//
//			for (unsigned int i = m_attribs[DART].begin(); i != m_attribs[DART].end(); m_attribs[DART].next(i))
//			{
//				unsigned int& idx = m_embeddings[orbit]->operator [](i);
//				unsigned int jdx = oldnew[idx];
//				if ((jdx != 0xffffffff) && (jdx != idx))
//					idx = jdx;
//			}
//		}
//	}

	// compact embedding attribs
	std::vector< std::vector<unsigned int>* > oldnews;
	oldnews.resize(NB_ORBITS);
822
823
824
825
	for (unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
	{
		if ((orbit != DART) && (isOrbitEmbedded(orbit)))
		{
826
827
828
829
			oldnews[orbit] = new std::vector<unsigned int>;
			m_attribs[orbit].compact(*(oldnews[orbit]));
		}
	}
830

831
832
833
834
835
836
	// 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)))
837
			{
838
839
				unsigned int& idx = m_embeddings[orbit]->operator[](i);
				unsigned int jdx = oldnews[orbit]->operator[](idx);
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
840
				if ((jdx != 0xffffffff) && (jdx != idx))
841
842
843
844
					idx = jdx;
			}
		}
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
845

846
847
848
849
850
	// delete allocated vectors
	for (unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
		if ((orbit != DART) && (isOrbitEmbedded(orbit)))
			delete[] oldnews[orbit];

851
	//compacting the topo
852
	std::vector<unsigned int> oldnew;
853
854
	m_attribs[DART].compact(oldnew);

855
856
857
858
859
860
861
862
	// update MR indices to attribs[DART]
	if (m_isMultiRes)
	{
		unsigned int nbl = m_mrDarts.size();
		for (unsigned int i = m_mrattribs.begin(); i != m_mrattribs.end(); m_mrattribs.next(i))
		{
			for (unsigned int j=0; j<nbl; ++j)
			{
863
				unsigned int d_index = m_mrDarts[j]->operator[](i);
864
				if (d_index != oldnew[d_index])
865
					m_mrDarts[j]->operator[](i) = oldnew[d_index];
866
867
868
869
			}
		}
	}

Pierre Kraemer's avatar
Pierre Kraemer committed
870
	// update topo relations from real map
871
	compactTopoRelations(oldnewMR);
872
873
874
875

//	dumpAttributesAndMarkers();
}

Pierre Kraemer's avatar
Pierre Kraemer committed
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
/****************************************
 *           DARTS TRAVERSALS           *
 ****************************************/

bool GenericMap::foreach_dart(FunctorType& f, const FunctorSelect& good)
{
	for (Dart d = begin(); d != end(); next(d))
	{
		if(good(d))
			if (f(d))
				return true;
	}
	return false;
}

} // namespace CGoGN
Sylvain Thery's avatar
Sylvain Thery committed
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065

//
//bool GenericMap::saveMapXml(const std::string& filename, bool compress)
//{
//	xmlTextWriterPtr writer = xmlNewTextWriterFilename(filename.c_str(), 0);
//	xmlTextWriterStartDocument(writer,NULL,NULL,NULL);
//
//	// Entete
//	int rc = xmlTextWriterStartElement(writer, BAD_CAST "CGoGN_Map");
//	rc = xmlTextWriterWriteAttribute(writer,  BAD_CAST "type", BAD_CAST this->mapTypeName().c_str());
//	rc = xmlTextWriterWriteFormatAttribute(writer,  BAD_CAST "nb_max_orbits","%u",NB_ORBITS);
//
//	// save m_attribs
//	for (unsigned int i=0; i<NB_ORBITS; ++i)
//	{
//		m_attribs[i].saveXml(writer,i);
//	}
//
//	// save m_orbMarker
//	rc = xmlTextWriterStartElement(writer, BAD_CAST "Orbit_MarkerSet");
//	for (unsigned int i=0; i<NB_ORBITS; ++i)
//	{
//		int rc = xmlTextWriterStartElement(writer, BAD_CAST "MarkerSet");
//		rc = xmlTextWriterWriteFormatAttribute(writer,  BAD_CAST "orbit", "%u", i);
//		rc = xmlTextWriterWriteFormatAttribute(writer,  BAD_CAST "val", "%u", m_marksets[i][0].getMarkVal());
////		rc = xmlTextWriterWriteAttribute(writer,  BAD_CAST "bin", BAD_CAST m_marksets[i][0].getMarkerAsBinaryString().c_str());
////		m_marksets[i] ;
//		rc = xmlTextWriterEndElement(writer);
//	}
//	rc = xmlTextWriterEndElement(writer);
//
//	// save m_BoundaryMarkerVal
////	rc = xmlTextWriterStartElement(writer, BAD_CAST "Boundary_Marker");
////	rc = xmlTextWriterWriteFormatAttribute(writer,  BAD_CAST "val", "%u", m_BoundaryMarkerVal.getMarkVal());
////	rc = xmlTextWriterEndElement(writer);
//
//	rc = xmlTextWriterEndElement(writer); // of map
//
//	xmlTextWriterEndDocument(writer);
//	xmlFreeTextWriter(writer);
//
//	return true;
//}
//
//
//bool GenericMap::loadMapXml(const std::string& filename, bool compress)
//{
//
//	xmlDocPtr doc = xmlReadFile(filename.c_str(), NULL, 0);
//	xmlNodePtr map_node = xmlDocGetRootElement(doc);
//
//	// Check if it is a CGoGN_Map file
////	if (strcmp((char*)(map_node->name),(char*)"CGoGN_Map")!=0)
//	if (!chechXmlNode(map_node,"CGoGN_Map"))
//	{
//		CGoGNerr << "Wrong xml format: Root node != CGoGN_Map"<< CGoGNendl;
//		return false;
//	}
//
//	// check the file type
//	xmlChar *prop = xmlGetProp(map_node, BAD_CAST "type");
//	CGoGNout << "Loading "<< prop <<" xml file"<< CGoGNendl;
//
//	// check the nb max orbits
//	prop = xmlGetProp(map_node, BAD_CAST "nb_max_orbits");
//	unsigned int nbo = atoi((char*)prop);
//	if (nbo != NB_ORBITS)
//	{
//		CGoGNerr << "Wrong nb max orbits in xml map"<< CGoGNendl;
//		return false;
//	}
//
//	/***********************************************
//	*	 			 load attributs
//	************************************************/
//	for (xmlNode* cur_node = map_node->children; cur_node; cur_node = cur_node->next)
//	{
//		// for each attribute
////		if (strcmp((char*)(cur_node->name),(char*)"Attributes_Container")==0)
//		if (chechXmlNode(cur_node, "Attributes_Container"))
//		{
//			CGoGNout << "LOAD ATTRIBUT"<< CGoGNendl;
//			// get the orbit id
//			unsigned int id = AttributeContainer::getIdXmlNode(cur_node);
//			// and load container
//			unsigned int nba = m_attribs[id].getNbAttributes();
//
//
//			CGoGNout << "attribut "<<id<<" size="<< m_attribs[id].size()<< "  nbatt="<< nba<< CGoGNendl;
//
//			m_attribs[id].loadXml(cur_node);
//		}
//	}
//
//	/***********************************************
//	*   creation of the m_embeddings pointers table
//	************************************************/
////	// get attribute names of dart orbit
////	AttributeContainer& contDart = m_attribs[DART] ;
////	std::vector< std::string > tableNames;
////	contDart.getAttributesStrings(tableNames);
////
////	// find orbit frome name and store pointer in right place
////	for (unsigned int i = 0; i< tableNames.size(); ++i)
////	{
//////		CGoGNout << i <<" : "<< tableNames[i]<<CGoGNendl;
////
////		std::string& name = tableNames[i];
////		std::string is_an_emb = name.substr(0,4);
////		if (is_an_emb == "EMB_")
////		{
////			AttributeMultiVector<unsigned int>& amv = contDart.getDataVector<unsigned int>(i) ;
////
////			std::string orbitname = name.substr(4, name.size()-4);
////			std::istringstream iss(orbitname);
////			unsigned int orbit;
////			iss >> orbit;
////
////			m_embeddings[orbit] = &amv ;
////		}
////	}
////
////
//	update_m_emb_afterLoad();
//
//	/***********************************************
//	*	 load Orbit_MarkerSet & BoundaryMarker
//	************************************************/
//	xmlNode* cur_node = map_node->children;
//	bool read1=false;
//	bool read2=false;
//	while (!(read1 || read2) && cur_node)// scan nodes to find the one with right name
//	{
////		if (strcmp((char*)(cur_node->name),(char*)"Orbit_MarkerSet") == 0)
//		if (chechXmlNode(cur_node, "Orbit_MarkerSet"))
//		{
//			for (xmlNode* mark_node = cur_node->children; mark_node; mark_node = mark_node->next)
//			{
//				xmlChar* prop = xmlGetProp(mark_node, BAD_CAST "orbit");
//				unsigned int orb = atoi((char*)prop);
//				prop = xmlGetProp(mark_node, BAD_CAST "val");
//				unsigned int val = atoi((char*)prop);
//				m_marksets[orb][0].setMarkVal(val);
//			}
//			read1 =true;
//		}
//		else
//		{
////			if (strcmp((char*)(cur_node->name),(char*)"Boundary_Marker") == 0)
//			CGoGNout << "Orbit_MarkerSet"<<CGoGNendl;
////			if (chechXmlNode(cur_node, "Boundary_Marker"))
////			{
////				xmlChar* prop = xmlGetProp(cur_node, BAD_CAST "val");
////				unsigned int val = atoi((char*)prop);
////				m_BoundaryMarkerVal.setMarkVal(val);
////				read2 =true;
////			}
//		}
//		// next node
//		cur_node = cur_node->next;
//	}
//
//	if (!(read1 && read2))
//	{
//		CGoGNerr <<"Error reading Marker in xml node"<<CGoGNendl;
//		return false;
//	}
//
//	xmlFreeDoc(doc);
//
//	update_m_emb_afterLoad();
//
//	return true ;
//}