genericmap.cpp 29 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
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)
42
		m_attributes_registry_map = new std::map<std::string, RegisteredBaseAttribute*> ;
43

44
45

	m_nbInstances++;
Pierre Kraemer's avatar
Pierre Kraemer committed
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
	// register all known types
	registerAttribute<Dart>("Dart");
	registerAttribute<Mark>("Mark");

	registerAttribute<long>("long");
	registerAttribute<int>("int");
	registerAttribute<short>("short");
	registerAttribute<char>("char");

	registerAttribute<unsigned long>("unsigned long");
	registerAttribute<unsigned int>("unsigned int");
	registerAttribute<unsigned short>("unsigned short");
	registerAttribute<unsigned char>("unsigned char");

	registerAttribute<Geom::Vec2f>(Geom::Vec2f::CGoGNnameOfType());
	registerAttribute<Geom::Vec3f>(Geom::Vec3f::CGoGNnameOfType());
	registerAttribute<Geom::Vec4f>(Geom::Vec4f::CGoGNnameOfType());

	registerAttribute<Geom::Vec2d>(Geom::Vec2d::CGoGNnameOfType());
	registerAttribute<Geom::Vec3d>(Geom::Vec3d::CGoGNnameOfType());
	registerAttribute<Geom::Vec4d>(Geom::Vec4d::CGoGNnameOfType());

	registerAttribute<Geom::Matrix33f>(Geom::Matrix33f::CGoGNnameOfType());
	registerAttribute<Geom::Matrix44f>(Geom::Matrix44f::CGoGNnameOfType());

	registerAttribute<Geom::Matrix33d>(Geom::Matrix33d::CGoGNnameOfType());
	registerAttribute<Geom::Matrix44d>(Geom::Matrix44d::CGoGNnameOfType());

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

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

93
94
	// get & lock marker for boundary
	m_boundaryMarker =  m_marksets[DART][0].getNewMark();
95

Pierre Kraemer's avatar
Pierre Kraemer committed
96
97
98
99
}

GenericMap::~GenericMap()
{
100
101
102
	// release marker for boundary
	m_marksets[DART][0].releaseMark(m_boundaryMarker);

Pierre Kraemer's avatar
Pierre Kraemer committed
103
104
105
106
107
	for(unsigned int i = 0; i < NB_ORBITS; ++i)
	{
		if(isOrbitEmbedded(i))
			m_attribs[i].clear(true) ;
	}
108

109

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

114
115
116
117
118
	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() ;
119

120
121
122
123
		for(std::vector<CellMarkerGen*>::iterator it = cellMarkers[i].begin(); it != cellMarkers[i].end(); ++it)
			(*it)->setReleaseOnDestruct(false) ;
		cellMarkers[i].clear() ;
	}
124
125
126
127

	// clean type registry if necessary
	m_nbInstances--;
	if (m_nbInstances<=0)
128
	{
129
130
131
		for (std::map<std::string, RegisteredBaseAttribute*>::iterator it =  m_attributes_registry_map->begin(); it != m_attributes_registry_map->end(); ++it)
			delete it->second;

132
		delete m_attributes_registry_map;
133
134
		m_attributes_registry_map = NULL;
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
135
136
}

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

	if (m_isMultiRes)
158
		initMR() ;
159
160
}

Pierre Kraemer's avatar
Pierre Kraemer committed
161
162
163
164
/****************************************
 *           MULTIRES                   *
 ****************************************/

Pierre Kraemer's avatar
Pierre Kraemer committed
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
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))
	{
		std::cout << (*m_mrLevels)[i] << " / " ;
		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
185
186
	m_mrattribs.setRegistry(m_attributes_registry_map) ;

Pierre Kraemer's avatar
Pierre Kraemer committed
187
188
189
190
191
192
193
194
195
196
197
198
	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") ;
	addLevel() ;
	setCurrentLevel(0) ;
}

199
//TODO rename to addBackLevel() ?
Pierre Kraemer's avatar
Pierre Kraemer committed
200
201
void GenericMap::addLevel()
{
Pierre Kraemer's avatar
Pierre Kraemer committed
202
	unsigned int newLevel = m_mrDarts.size() ;
Pierre Kraemer's avatar
Pierre Kraemer committed
203
	std::stringstream ss ;
Pierre Kraemer's avatar
Pierre Kraemer committed
204
	ss << "MRdart_"<< newLevel ;
205
	AttributeMultiVector<unsigned int>* newAttrib = m_mrattribs.addAttribute<unsigned int>(ss.str()) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
206

207
	m_mrDarts.push_back(newAttrib) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
208
	m_mrNbDarts.push_back(0) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
209
210

	if(m_mrDarts.size() > 1)
Pierre Kraemer's avatar
Pierre Kraemer committed
211
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
212
213
		AttributeMultiVector<unsigned int>* prevAttrib = m_mrDarts[newLevel - 1] ;	// copy the indices of
		m_mrattribs.copyAttribute(newAttrib->getIndex(), prevAttrib->getIndex()) ;	// previous level into new level
214

215
216
		for(unsigned int i = m_mrattribs.begin(); i != m_mrattribs.end(); m_mrattribs.next(i))
		{
217
			//++(*m_mrLevels)[i]; // increase the insertion level of the darts
218
219
220
			unsigned int oldi = (*prevAttrib)[i] ;	// get the index of the dart in previous level
			(*newAttrib)[i] = copyDartLine(oldi) ;	// copy the dart and affect it to the new level
		}
Pierre Kraemer's avatar
Pierre Kraemer committed
221
222
223
	}
}

224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
void GenericMap::addFrontLevel()
{
	unsigned int newLevel = m_mrDarts.size() ;
	std::stringstream ss ;
	ss << "MRdart_"<< newLevel ;
	AttributeMultiVector<unsigned int>* newAttrib = m_mrattribs.addAttribute<unsigned int>(ss.str()) ;

	m_mrDarts.insert(m_mrDarts.begin(), newAttrib) ;
	m_mrNbDarts.insert(m_mrNbDarts.begin(), 0) ;

	if(m_mrDarts.size() > 1)
	{
		AttributeMultiVector<unsigned int>* prevAttrib = m_mrDarts[1] ;	// copy the indices of
		m_mrattribs.copyAttribute(newAttrib->getIndex(), prevAttrib->getIndex()) ;	// previous level into new level

		for(unsigned int i = m_mrattribs.begin(); i != m_mrattribs.end(); m_mrattribs.next(i))
		{
			unsigned int oldi = (*prevAttrib)[i] ;	// get the index of the dart in previous level
			(*newAttrib)[i] = copyDartLine(oldi) ;	// copy the dart and affect it to the new level
		}
	}
}

247
//TODO rename to removeBackLevel() ?
Pierre Kraemer's avatar
Pierre Kraemer committed
248
249
void GenericMap::removeLevel()
{
Pierre Kraemer's avatar
Pierre Kraemer committed
250
251
	unsigned int maxL = getMaxLevel() ;
	if(maxL > 0)
Pierre Kraemer's avatar
Pierre Kraemer committed
252
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
253
254
		AttributeMultiVector<unsigned int>* maxMR = m_mrDarts[maxL] ;
		AttributeMultiVector<unsigned int>* prevMR = m_mrDarts[maxL - 1] ;
Pierre Kraemer's avatar
Pierre Kraemer committed
255
256
		for(unsigned int i = m_mrattribs.begin(); i != m_mrattribs.end(); m_mrattribs.next(i))
		{
Pierre Kraemer's avatar
Pierre Kraemer committed
257
258
259
260
261
262
263
264
265
266
267
			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
268
269
		}

Pierre Kraemer's avatar
Pierre Kraemer committed
270
		m_mrattribs.removeAttribute<unsigned int>(maxMR->getIndex()) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
271
272
273
274
275
		m_mrDarts.pop_back() ;
		m_mrNbDarts.pop_back() ;

		if(m_mrCurrentLevel == maxL)
			--m_mrCurrentLevel ;
Pierre Kraemer's avatar
Pierre Kraemer committed
276
277
278
	}
}

279
280
281
282
283
void GenericMap::removeFrontLevel()
{
	std::cout << "TO DO" << std::endl;
}

Pierre Kraemer's avatar
Pierre Kraemer committed
284
285
286
287
/****************************************
 *        ATTRIBUTES MANAGEMENT         *
 ****************************************/

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
313
314
315
316
317
//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(std::vector<CellMarkerGen*>::iterator it = cellMarkers.begin(); it != cellMarkers.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
318

Pierre Kraemer's avatar
Pierre Kraemer committed
319
320
321
322
323
324
325
326
327
328
329
330
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)
331
332
333
334
335
		{
			unsigned int id = cont.getAttributeIndex(*it);
			std::cout << "    " << *it << " ("<<id<<")"<<std::endl ;
		}

Pierre Kraemer's avatar
Pierre Kraemer committed
336
337
338
339
340
341
342
343
344
345
346
347
348
		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 ;
}

349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
/****************************************
 *          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)
		{
364
365
366
367
368
			std::stringstream ss ;
			ss << "Mark_"<< th ;
			AttributeContainer& cellCont = m_attribs[i] ;
			AttributeMultiVector<Mark>* amvMark = cellCont.addAttribute<Mark>(ss.str()) ;
			m_markTables[i][th] = amvMark ;
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
		}
	}
}

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)
		{
387
388
389
390
391
			std::stringstream ss ;
			ss << "Mark_"<< th ;
			AttributeContainer& cellCont = m_attribs[i] ;
			cellCont.removeAttribute<Mark>(ss.str()) ;
			m_markTables[i][th] = NULL ;
392
393
394
		}
	}
}
Pierre Kraemer's avatar
Pierre Kraemer committed
395
396
397
398
399
400
401
402
403
404

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

bool GenericMap::saveMapBin(const std::string& filename)
{
	CGoGNostream fs(filename.c_str(), std::ios::out|std::ios::binary);
	if (!fs)
	{
405
		CGoGNerr << "Unable to open file for writing: " << filename << CGoGNendl;
Pierre Kraemer's avatar
Pierre Kraemer committed
406
407
408
409
410
		return false;
	}

	// Entete
	char* buff = new char[256];
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
411
	for (int i = 0; i < 256; ++i)
Pierre Kraemer's avatar
Pierre Kraemer committed
412
		buff[i] = char(255);
413
414
415
416
417
418
419
420
421
422
	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
423
424
	std::string mt = mapTypeName();
	const char* mtc = mt.c_str();
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
425
426
	memcpy(buff+32, mtc, mt.size()+1);
	unsigned int *buffi = reinterpret_cast<unsigned int*>(buff + 64);
Pierre Kraemer's avatar
Pierre Kraemer committed
427
428
429
430
	*buffi = NB_ORBITS;
	fs.write(reinterpret_cast<const char*>(buff), 256);
	delete buff;

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

435
	if (m_isMultiRes)
Sylvain Thery's avatar
Sylvain Thery committed
436
	{
437
438
		m_mrattribs.saveBin(fs, 00);

Sylvain Thery's avatar
Sylvain Thery committed
439
440
441
442
443
444
445
		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
446
447
448
449
450
	return true;
}

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

Pierre Kraemer's avatar
Pierre Kraemer committed
452
453
454
	CGoGNistream fs(filename.c_str(), std::ios::in|std::ios::binary);
	if (!fs)
	{
455
		CGoGNerr << "Unable to open file for loading" << CGoGNendl;
Pierre Kraemer's avatar
Pierre Kraemer committed
456
457
458
		return false;
	}

Sylvain Thery's avatar
Sylvain Thery committed
459
460
461
	GenericMap::clear(true);

	// read info
Pierre Kraemer's avatar
Pierre Kraemer committed
462
463
464
465
466
	char* buff = new char[256];
	fs.read(reinterpret_cast<char*>(buff), 256);

	std::string buff_str(buff);
	// Check file type
467
	if (m_isMultiRes)
Pierre Kraemer's avatar
Pierre Kraemer committed
468
	{
469
470
471
472
473
474
475
476
477
478
		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
479
	}
480
481
482
483
484
485
486
487
488
489
490
491
492
493
	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
494
495

	// Check map type
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
496
	buff_str = std::string(buff + 32);
Pierre Kraemer's avatar
Pierre Kraemer committed
497
498

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

	std::string fileType = buff_str;
Pierre Kraemer's avatar
Pierre Kraemer committed
501
502
503

	if (fileType != localType)
	{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
504
		CGoGNerr << "Not possible to load "<< fileType << " into " << localType << " object" << CGoGNendl;
Pierre Kraemer's avatar
Pierre Kraemer committed
505
506
507
508
		return false;
	}

	// Check max nb orbit
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
509
	unsigned int *ptr_nbo = reinterpret_cast<unsigned int*>(buff + 64);
Pierre Kraemer's avatar
Pierre Kraemer committed
510
511
512
	unsigned int nbo = *ptr_nbo;
	if (nbo != NB_ORBITS)
	{
513
		CGoGNerr << "Wrond max orbit number in file" << CGoGNendl;
Pierre Kraemer's avatar
Pierre Kraemer committed
514
515
516
517
		return  false;
	}

	// load attrib container
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
518
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
Pierre Kraemer's avatar
Pierre Kraemer committed
519
	{
520
		unsigned int id = AttributeContainer::loadBinId(fs);
Pierre Kraemer's avatar
Pierre Kraemer committed
521
522
523
		m_attribs[id].loadBin(fs);
	}

524
	if (m_isMultiRes)
Sylvain Thery's avatar
Sylvain Thery committed
525
526
	{
		AttributeContainer::loadBinId(fs); // not used but need to read to skip
527
528
		m_mrattribs.loadBin(fs);

Sylvain Thery's avatar
Sylvain Thery committed
529
530
531
532
533
534
535
536
		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));
	}


537
	// retrieve m_embeddings (from m_attribs)
Pierre Kraemer's avatar
Pierre Kraemer committed
538
	update_m_emb_afterLoad();
539

Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
540
	// 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
541
	update_topo_shortcuts();
Sylvain Thery's avatar
Sylvain Thery committed
542

Pierre Kraemer's avatar
Pierre Kraemer committed
543
544
545
	return true;
}

Sylvain Thery's avatar
Sylvain Thery committed
546
547
548
549
bool GenericMap::copyFrom(const GenericMap& map)
{
	if (mapTypeName() != map.mapTypeName())
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
550
		CGoGNerr << "try to copy from incompatible type map" << CGoGNendl;
Sylvain Thery's avatar
Sylvain Thery committed
551
552
553
554
555
556
557
558
559
560
561
562
563
564
		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
565
		unsigned int nb = map.m_mrNbDarts.size();
Sylvain Thery's avatar
Sylvain Thery committed
566
		m_mrNbDarts.resize(nb);
Pierre Kraemer's avatar
Pierre Kraemer committed
567
		for (unsigned int i = 0; i < nb; ++i)
Sylvain Thery's avatar
Sylvain Thery committed
568
569
570
571
572
573
574
575
576
577
578
579
			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
580
581
582
583
584
585
586
587
588
589
590
591
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
592
		std::string sub = listeNames[i].substr(0, listeNames[i].size() - 1);
Sylvain Thery's avatar
Sylvain Thery committed
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
		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
614
			std::string sub = listeNames[i].substr(0, 5);
Sylvain Thery's avatar
Sylvain Thery committed
615
616
617
618
			if (sub == "Mark_")
			{
				// get thread number
				unsigned int thread = listeNames[i][5]-'0';
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
619
				if (listeNames[i].size() > 6) 					// thread number is >9
Pierre Kraemer's avatar
Pierre Kraemer committed
620
					thread = 10*thread + (listeNames[i][6]-'0');
Sylvain Thery's avatar
Sylvain Thery committed
621
622
623
624

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

Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
625
				if ((orbit == DART) && (thread == 0))	// for Marker of dart of thread O keep the boundary marker
Sylvain Thery's avatar
Sylvain Thery committed
626
627
628
				{
					Mark m(m_boundaryMarker);
					m.invert();
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
629
					for (unsigned int i = cont.begin(); i != cont.end(); cont.next(i))
Sylvain Thery's avatar
Sylvain Thery committed
630
631
632
633
						amvMark->operator[](i).unsetMark(m);
				}
				else								// for others clear all
				{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
634
					for (unsigned int i = cont.begin(); i != cont.end(); cont.next(i))
Sylvain Thery's avatar
Sylvain Thery committed
635
636
637
638
639
						amvMark->operator[](i).clear();
				}
			}
		}
	}
640
641
642
643
644

	if (m_isMultiRes)
	{
		std::vector<std::string> names;
		m_mrattribs.getAttributesNames(names);
Pierre Kraemer's avatar
Pierre Kraemer committed
645
646
		m_mrDarts.resize(names.size() - 1);
		for (unsigned int i = 0; i < m_mrDarts.size(); ++i)
647
648
649
650
651
652
			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
653
			if (sub == "MRLevel")
654
				m_mrLevels = m_mrattribs.getDataVector<unsigned int>(i);
655

Pierre Kraemer's avatar
Pierre Kraemer committed
656
			if (sub == "MRdart_")
657
658
			{
				sub = names[i].substr(7);	// compute number following MT_Dart_
Pierre Kraemer's avatar
Pierre Kraemer committed
659
660
661
662
				unsigned int idx = 0;
				for (unsigned int j = 0; j < sub.length(); j++)
					idx = 10*idx + (sub[j]-'0');
				if (idx < names.size() - 1)
663
					m_mrDarts[idx] = m_mrattribs.getDataVector<unsigned int>(i);
664
				else
Pierre Kraemer's avatar
Pierre Kraemer committed
665
					CGoGNerr << "Warning problem updating MR_DARTS" << CGoGNendl;
666
667
668
			}
		}
		// check if all pointers are != NULL
Pierre Kraemer's avatar
Pierre Kraemer committed
669
		for (unsigned int i = 0; i < m_mrDarts.size(); ++i)
670
671
		{
			if (m_mrDarts[i] == NULL)
Pierre Kraemer's avatar
Pierre Kraemer committed
672
				CGoGNerr << "Warning problem MR_DARTS = NULL" << CGoGNendl;
673
674
		}
	}
Sylvain Thery's avatar
Sylvain Thery committed
675
676
677
678
}

void GenericMap::dumpAttributesAndMarkers()
{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
679
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
Sylvain Thery's avatar
Sylvain Thery committed
680
681
	{
		std::vector<std::string> names;
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
682
		names.reserve(32); 				//just to limit reallocation
Sylvain Thery's avatar
Sylvain Thery committed
683
684
		m_attribs[i].getAttributesNames(names);
		unsigned int nb = names.size();
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
685
		if (nb > 0)
Sylvain Thery's avatar
Sylvain Thery committed
686
		{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
687
			CGoGNout << "ORBIT "<< i << CGoGNendl;
Sylvain Thery's avatar
Sylvain Thery committed
688
689
690
			std::vector<std::string> types;
			types.reserve(nb);
			m_attribs[i].getAttributesTypes(types);
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
691
692
			for (unsigned int j = 0; j < nb; ++j)
				CGoGNout << "    " << j << " : " << types[j] << " " << names[j] << CGoGNendl;
Sylvain Thery's avatar
Sylvain Thery committed
693
694
695
		}
	}
	CGoGNout << "RESERVED MARKERS "<< CGoGNendl;
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
696
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
Sylvain Thery's avatar
Sylvain Thery committed
697
	{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
698
		for (unsigned int j = 0; j < NB_THREAD; ++j)
Sylvain Thery's avatar
Sylvain Thery committed
699
700
701
702
		{
			MarkSet ms = m_marksets[i][j];
			if (!ms.isClear())
			{
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
703
				CGoGNout << "Orbit " << i << "  thread " << j << " : ";
Sylvain Thery's avatar
Sylvain Thery committed
704
				Mark m(1);
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
705
				for (unsigned i = 0; i < Mark::getNbMarks(); ++i)
Sylvain Thery's avatar
Sylvain Thery committed
706
707
				{
					if (ms.testMark(m))
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
708
						CGoGNout << m.getMarkVal() << ", ";
Sylvain Thery's avatar
Sylvain Thery committed
709
710
711
712
713
714
715
716
					m.setMarkVal(m.getMarkVal()<<1);
				}
				CGoGNout << CGoGNendl;
			}
		}
	}
}

717
718
void GenericMap::compact()
{
719
720
721
722
	// if MR compact the MR attrib container
	std::vector<unsigned int> oldnewMR;
	if (m_isMultiRes)
		m_mrattribs.compact(oldnewMR);
723
724

	// compacting the orbits attributes
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
//	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);
744
745
746
747
	for (unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
	{
		if ((orbit != DART) && (isOrbitEmbedded(orbit)))
		{
748
749
750
751
			oldnews[orbit] = new std::vector<unsigned int>;
			m_attribs[orbit].compact(*(oldnews[orbit]));
		}
	}
752

753
754
755
756
757
758
	// 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)))
759
			{
760
761
				unsigned int& idx = m_embeddings[orbit]->operator[](i);
				unsigned int jdx = oldnews[orbit]->operator[](idx);
Pierre Kraemer's avatar
merges    
Pierre Kraemer committed
762
				if ((jdx != 0xffffffff) && (jdx != idx))
763
764
765
766
					idx = jdx;
			}
		}
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
767

768
769
770
771
772
	// delete allocated vectors
	for (unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
		if ((orbit != DART) && (isOrbitEmbedded(orbit)))
			delete[] oldnews[orbit];

773
	//compacting the topo
774
	std::vector<unsigned int> oldnew;
775
776
	m_attribs[DART].compact(oldnew);

777
778
779
780
781
782
783
784
	// 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)
			{
785
				unsigned int d_index = m_mrDarts[j]->operator[](i);
786
				if (d_index != oldnew[d_index])
787
					m_mrDarts[j]->operator[](i) = oldnew[d_index];
788
789
790
791
			}
		}
	}

Pierre Kraemer's avatar
Pierre Kraemer committed
792
	// update topo relations from real map
793
	compactTopoRelations(oldnewMR);
794
795
796
797

//	dumpAttributesAndMarkers();
}

Pierre Kraemer's avatar
Pierre Kraemer committed
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
/****************************************
 *           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;
}

813
814
void GenericMap::boundaryUnmarkAll()
{
815
	AttributeContainer& cont = getAttributeContainer<DART>() ;
816
817
818
819
	for (unsigned int i = cont.begin(); i != cont.end(); cont.next(i))
		m_markTables[DART][0]->operator[](i).unsetMark(m_boundaryMarker);
}

Pierre Kraemer's avatar
Pierre Kraemer committed
820
} // namespace CGoGN
Sylvain Thery's avatar
Sylvain Thery committed
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
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

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