genericmap.hpp 19.9 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
* Contact information: cgogn@unistra.fr                                        *
*                                                                              *
*******************************************************************************/

25
26
#include "Topology/generic/dartmarker.h"
#include "Topology/generic/traversorCell.h"
27
#include "Topology/generic/traversorFactory.h"
28

Pierre Kraemer's avatar
Pierre Kraemer committed
29
30
namespace CGoGN
{
Pierre Kraemer's avatar
Pierre Kraemer committed
31

32
33
34
35
36
37
38
39
40
41
42
43
44
45
/****************************************
 *           MULTIRES                   *
 ****************************************/

inline unsigned int GenericMap::getCurrentLevel()
{
	return m_mrCurrentLevel ;
}

inline void GenericMap::setCurrentLevel(unsigned int l)
{
	if(l < m_mrDarts.size())
		m_mrCurrentLevel = l ;
	else
Pierre Kraemer's avatar
Pierre Kraemer committed
46
		CGoGNout << "setCurrentLevel : try to access nonexistent resolution level" << CGoGNendl ;
47
48
}

Pierre Kraemer's avatar
Pierre Kraemer committed
49
50
inline void GenericMap::incCurrentLevel()
{
51
	if(m_mrCurrentLevel < m_mrDarts.size() - 1)
Pierre Kraemer's avatar
Pierre Kraemer committed
52
53
54
55
56
57
58
59
60
61
62
63
64
		++m_mrCurrentLevel ;
	else
		CGoGNout << "incCurrentLevel : already at maximum resolution level" << CGoGNendl ;
}

inline void GenericMap::decCurrentLevel()
{
	if(m_mrCurrentLevel > 0)
		--m_mrCurrentLevel ;
	else
		CGoGNout << "decCurrentLevel : already at minimum resolution level" << CGoGNendl ;
}

65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
inline void GenericMap::pushLevel()
{
	m_mrLevelStack.push_back(m_mrCurrentLevel) ;
}

inline void GenericMap::popLevel()
{
	m_mrCurrentLevel = m_mrLevelStack.back() ;
	m_mrLevelStack.pop_back() ;
}

inline unsigned int GenericMap::getMaxLevel()
{
	return m_mrDarts.size() - 1 ;
}

81

Pierre Kraemer's avatar
Pierre Kraemer committed
82
83
84
85
86
87
/****************************************
 *           DARTS MANAGEMENT           *
 ****************************************/

inline Dart GenericMap::newDart()
{
88
	unsigned int di = m_attribs[DART].insertLine();		// insert a new dart line
Pierre Kraemer's avatar
Pierre Kraemer committed
89
	for(unsigned int i = 0; i < NB_ORBITS; ++i)
Pierre Kraemer's avatar
Pierre Kraemer committed
90
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
91
92
		if (m_embeddings[i])								// set all its embeddings
			(*m_embeddings[i])[di] = EMBNULL ;				// to EMBNULL
Pierre Kraemer's avatar
Pierre Kraemer committed
93
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
94

Pierre Kraemer's avatar
Pierre Kraemer committed
95
96
	if (m_isMultiRes)
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
97
98
		unsigned int mrdi = m_mrattribs.insertLine() ;		// insert a new MRdart line
		(*m_mrLevels)[mrdi] = m_mrCurrentLevel ;			// set the introduction level of the dart
Pierre Kraemer's avatar
Pierre Kraemer committed
99
		m_mrNbDarts[m_mrCurrentLevel]++ ;
Pierre Kraemer's avatar
Pierre Kraemer committed
100

Pierre Kraemer's avatar
Pierre Kraemer committed
101
		for(unsigned int i = 0; i < m_mrCurrentLevel; ++i)	// for all previous levels
102
			(*m_mrDarts[i])[mrdi] = MRNULL ;					// this MRdart does not exist
Pierre Kraemer's avatar
Pierre Kraemer committed
103

untereiner's avatar
untereiner committed
104
105
		for(unsigned int i = m_mrCurrentLevel; i < m_mrDarts.size(); ++i)	// for all levels from current to max
  			(*m_mrDarts[i])[mrdi] = di ;									// make this MRdart point to the new dart line
106

Pierre Kraemer's avatar
Pierre Kraemer committed
107
108
		return Dart::create(mrdi) ;
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
109

Pierre Kraemer's avatar
Pierre Kraemer committed
110
	return Dart::create(di) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
111
112
113
114
}

inline void GenericMap::deleteDart(Dart d)
{
Pierre Kraemer's avatar
Pierre Kraemer committed
115
116
	if(m_isMultiRes)
	{
untereiner's avatar
untereiner committed
117
		unsigned int index = (*m_mrDarts[m_mrCurrentLevel])[d.index] ;
118

119
		if(getDartLevel(d) > m_mrCurrentLevel)
120
		{
121
122
123
124
125
126
127
128
			unsigned int di = (*m_mrDarts[m_mrCurrentLevel + 1])[d.index];
			// si le brin de niveau i pointe sur le meme brin que le niveau i-1
			if(di != index)
			{
				if(isDartValid(d))//index))
					deleteDartLine(index) ;
			}

untereiner's avatar
untereiner committed
129
130
			(*m_mrDarts[m_mrCurrentLevel])[d.index] = MRNULL ;
			return;
131
		}
untereiner's avatar
untereiner committed
132

133
134
135
		// a MRdart can only be deleted on its insertion level
		if(getDartLevel(d) == m_mrCurrentLevel)
		{
136
			if(isDartValid(d))//index))
137
				deleteDartLine(index) ;
untereiner's avatar
untereiner committed
138
139
140

			m_mrattribs.removeLine(d.index);
			m_mrNbDarts[m_mrCurrentLevel]--;
141
		}
untereiner's avatar
untereiner committed
142
143
		else
		{
144
			unsigned int di = (*m_mrDarts[m_mrCurrentLevel - 1])[d.index];
145
			// si le brin de niveau i pointe sur le meme brin que le niveau i-1
untereiner's avatar
untereiner committed
146
147
			if(di != index)
			{
148
				if(isDartValid(d))//index))
untereiner's avatar
untereiner committed
149
150
151
					deleteDartLine(index) ;
			}

152
153
154
155
156
			for(unsigned int i = m_mrCurrentLevel; i <= getMaxLevel(); ++i) // for all levels from current to max
			{
				(*m_mrDarts[i])[d.index] = di ; //copy the index from previous level
			}
		}
Pierre Kraemer's avatar
Pierre Kraemer committed
157
158
159
160
161
162
163
164
165
166
	}
	else
		deleteDartLine(dartIndex(d)) ;
}

inline void GenericMap::deleteDartLine(unsigned int index)
{
	m_attribs[DART].removeLine(index) ;	// free the dart line

	for (unsigned int t = 0; t < m_nbThreads; ++t)	// clear markers of
167
		(*m_markTables[DART][t])[index].clear() ;		// the removed dart
Pierre Kraemer's avatar
merge..    
Pierre Kraemer committed
168

Pierre Kraemer's avatar
Pierre Kraemer committed
169
170
	for(unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
171
		if (m_embeddings[orbit])									// for each embedded orbit
Pierre Kraemer's avatar
Pierre Kraemer committed
172
		{
Pierre Kraemer's avatar
Pierre Kraemer committed
173
			unsigned int emb = (*m_embeddings[orbit])[index] ;		// get the embedding of the dart
Pierre Kraemer's avatar
Pierre Kraemer committed
174
			if(emb != EMBNULL)
Pierre Kraemer's avatar
merge..    
Pierre Kraemer committed
175
			{
Pierre Kraemer's avatar
Pierre Kraemer committed
176
				if(m_attribs[orbit].unrefLine(emb))					// unref the pointed embedding line
Pierre Kraemer's avatar
merge..    
Pierre Kraemer committed
177
				{
Pierre Kraemer's avatar
Pierre Kraemer committed
178
					for (unsigned int t = 0; t < m_nbThreads; ++t)	// and clear its markers if it was
179
						(*m_markTables[orbit][t])[emb].clear() ;	// its last unref (and was thus freed)
Pierre Kraemer's avatar
merge..    
Pierre Kraemer committed
180
181
				}
			}
Pierre Kraemer's avatar
Pierre Kraemer committed
182
183
		}
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
184
}
Pierre Kraemer's avatar
Pierre Kraemer committed
185

Pierre Kraemer's avatar
Pierre Kraemer committed
186
inline unsigned int GenericMap::copyDartLine(unsigned int index)
Pierre Kraemer's avatar
Pierre Kraemer committed
187
{
188
189
	unsigned int newindex = m_attribs[DART].insertLine() ;	// create a new dart line
	m_attribs[DART].copyLine(newindex, index) ;				// copy the given dart line
Pierre Kraemer's avatar
Pierre Kraemer committed
190
191
	for(unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
192
193
194
195
196
197
		if (m_embeddings[orbit])
		{
			unsigned int emb = (*m_embeddings[orbit])[newindex] ;	// add a ref to the cells pointed
			if(emb != EMBNULL)										// by the new dart line
				m_attribs[orbit].refLine(emb) ;
		}
Pierre Kraemer's avatar
Pierre Kraemer committed
198
199
	}
	return newindex ;
Pierre Kraemer's avatar
Pierre Kraemer committed
200
201
}

Pierre Kraemer's avatar
Pierre Kraemer committed
202
203
inline void GenericMap::duplicateDart(Dart d)
{
204
	assert(getDartLevel(d) <= m_mrCurrentLevel || !"duplicateDart : called with a dart inserted after current level") ;
Pierre Kraemer's avatar
Pierre Kraemer committed
205

206
	if(getDartLevel(d) == m_mrCurrentLevel)	// no need to duplicate
Pierre Kraemer's avatar
Pierre Kraemer committed
207
208
209
210
		return ;								// a dart from its insertion level

	unsigned int oldindex = dartIndex(d) ;

211
212
213
214
215
	if(m_mrCurrentLevel > 0)
	{
		if((*m_mrDarts[m_mrCurrentLevel - 1])[d.index] != oldindex)	// no need to duplicate if the dart is already
			return ;												// duplicated with respect to previous level
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
216
217
218

	unsigned int newindex = copyDartLine(oldindex) ;

219
	for(unsigned int i = m_mrCurrentLevel; i <= getMaxLevel(); ++i) // for all levels from current to max
Pierre Kraemer's avatar
Pierre Kraemer committed
220
221
222
223
224
225
	{
		assert((*m_mrDarts[i])[d.index] == oldindex || !"duplicateDart : dart was already duplicated on a greater level") ;
		(*m_mrDarts[i])[d.index] = newindex ;						// make this MRdart points to the new dart line
	}
}

untereiner's avatar
untereiner committed
226
inline void GenericMap::duplicateDartAtOneLevel(Dart d, unsigned int level)
227
{
untereiner's avatar
untereiner committed
228
	(*m_mrDarts[level])[d.index] = copyDartLine(dartIndex(d)) ;
229
230
}

Sylvain Thery's avatar
Sylvain Thery committed
231
inline unsigned int GenericMap::dartIndex(Dart d) const
Pierre Kraemer's avatar
Pierre Kraemer committed
232
{
Pierre Kraemer's avatar
Pierre Kraemer committed
233
234
235
236
237
	if (m_isMultiRes)
		return (*m_mrDarts[m_mrCurrentLevel])[d.index] ;
	return d.index;
}

238
inline unsigned int GenericMap::getDartLevel(Dart d) const
Pierre Kraemer's avatar
Pierre Kraemer committed
239
240
241
242
{
	return (*m_mrLevels)[d.index] ;
}

243
244
inline void GenericMap::incDartLevel(Dart d) const
{
245
	++((*m_mrLevels)[d.index]) ;
246
247
248
}


Pierre Kraemer's avatar
Pierre Kraemer committed
249
250
251
252
253
254
255
256
257
258
259
inline unsigned int GenericMap::getNbInsertedDarts(unsigned int level)
{
	if(level < m_mrDarts.size())
		return m_mrNbDarts[level] ;
	else
		return 0 ;
}

inline unsigned int GenericMap::getNbDarts(unsigned int level)
{
	if(level < m_mrDarts.size())
Pierre Kraemer's avatar
Pierre Kraemer committed
260
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
261
262
263
264
		unsigned int nb = 0 ;
		for(unsigned int i = 0; i <= level; ++i)
			nb += m_mrNbDarts[i] ;
		return nb ;
Pierre Kraemer's avatar
Pierre Kraemer committed
265
	}
Pierre Kraemer's avatar
Pierre Kraemer committed
266
267
268
269
270
271
272
273
274
	else
		return 0 ;
}

inline unsigned int GenericMap::getNbDarts()
{
	if(m_isMultiRes)
		return getNbDarts(m_mrCurrentLevel) ;

275
	return m_attribs[DART].size() ;
Pierre Kraemer's avatar
Pierre Kraemer committed
276
277
}

Pierre Kraemer's avatar
Pierre Kraemer committed
278
279
280
281
282
inline bool GenericMap::isDartValid(Dart d)
{
	return !d.isNil() && m_attribs[DART].used(dartIndex(d)) ;
}

Pierre Kraemer's avatar
Pierre Kraemer committed
283
284
285
286
/****************************************
 *         EMBEDDING MANAGEMENT         *
 ****************************************/

287
288
template <unsigned int ORBIT>
inline bool GenericMap::isOrbitEmbedded() const
Pierre Kraemer's avatar
Pierre Kraemer committed
289
{
290
	return (ORBIT == DART) || (m_embeddings[ORBIT] != NULL) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
291
292
}

293
inline bool GenericMap::isOrbitEmbedded(unsigned int orbit) const
Pierre Kraemer's avatar
Pierre Kraemer committed
294
{
295
	return (orbit == DART) || (m_embeddings[orbit] != NULL) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
296
297
}

298
299
template <unsigned int ORBIT>
inline unsigned int GenericMap::getEmbedding(Dart d)
Pierre Kraemer's avatar
Pierre Kraemer committed
300
{
301
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
Pierre Kraemer's avatar
Pierre Kraemer committed
302

303
304
	unsigned int d_index = dartIndex(d);

305
	if (ORBIT == DART)
306
		return d_index;
Pierre Kraemer's avatar
Pierre Kraemer committed
307

308
	return (*m_embeddings[ORBIT])[d_index] ;
Pierre Kraemer's avatar
Pierre Kraemer committed
309
310
}

311
312
313
template <unsigned int ORBIT>
void GenericMap::setDartEmbedding(Dart d, unsigned int emb)
{
314
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332

	unsigned int old = getEmbedding<ORBIT>(d);

	if (old == emb)	// if same emb
		return;		// nothing to do

	if (old != EMBNULL)	// if different
	{
		if(m_attribs[ORBIT].unrefLine(old))	// then unref the old emb
		{
			for (unsigned int t = 0; t < m_nbThreads; ++t)	// clear the markers if it was the
				(*m_markTables[ORBIT][t])[old].clear();		// last unref of the line
		}
	}

	if (emb != EMBNULL)
		m_attribs[ORBIT].refLine(emb);	// ref the new emb

Pierre Kraemer's avatar
Pierre Kraemer committed
333
	(*m_embeddings[ORBIT])[dartIndex(d)] = emb ; // finally affect the embedding to the dart
334
335
}

336
337
338
339
template <unsigned int ORBIT>
void GenericMap::initDartEmbedding(Dart d, unsigned int emb)
{
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
340
341
	if(emb != EMBNULL)
		m_attribs[ORBIT].refLine(emb);	// ref the new emb
342
343
344
	(*m_embeddings[ORBIT])[dartIndex(d)] = emb ; // affect the embedding to the dart
}

345
346
template <unsigned int ORBIT>
inline void GenericMap::copyDartEmbedding(Dart dest, Dart src)
Pierre Kraemer's avatar
Pierre Kraemer committed
347
{
348
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
349
	setDartEmbedding<ORBIT>(dest, getEmbedding<ORBIT>(src));
Pierre Kraemer's avatar
Pierre Kraemer committed
350
351
}

352
353
template <unsigned int ORBIT>
inline unsigned int GenericMap::newCell()
Pierre Kraemer's avatar
Pierre Kraemer committed
354
{
355
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
356
	return m_attribs[ORBIT].insertLine();
Pierre Kraemer's avatar
Pierre Kraemer committed
357
358
}

359
template <unsigned int ORBIT>
360
inline void GenericMap::setOrbitEmbedding(Dart d, unsigned int em)
Pierre Kraemer's avatar
Pierre Kraemer committed
361
{
362
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
363
364
	FunctorSetEmb<GenericMap, ORBIT> fsetemb(*this, em);
	foreach_dart_of_orbit<ORBIT>(d, fsetemb);
Pierre Kraemer's avatar
Pierre Kraemer committed
365
366
}

367
template <unsigned int ORBIT>
368
369
370
371
372
373
374
375
inline void GenericMap::initOrbitEmbedding(Dart d, unsigned int em)
{
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
	FunctorInitEmb<GenericMap, ORBIT> fsetemb(*this, em);
	foreach_dart_of_orbit<ORBIT>(d, fsetemb);
}

template <unsigned int ORBIT>
376
inline unsigned int GenericMap::setOrbitEmbeddingOnNewCell(Dart d)
Pierre Kraemer's avatar
Pierre Kraemer committed
377
{
378
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
379
	unsigned int em = newCell<ORBIT>();
380
381
382
383
384
385
386
387
388
389
	setOrbitEmbedding<ORBIT>(d, em);
	return em;
}

template <unsigned int ORBIT>
inline unsigned int GenericMap::initOrbitEmbeddingNewCell(Dart d)
{
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
	unsigned int em = newCell<ORBIT>();
	initOrbitEmbedding<ORBIT>(d, em);
Pierre Kraemer's avatar
Pierre Kraemer committed
390
391
392
	return em;
}

393
394
template <unsigned int ORBIT>
inline void GenericMap::copyCell(Dart d, Dart e)
Pierre Kraemer's avatar
Pierre Kraemer committed
395
{
396
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
397
398
	unsigned int dE = getEmbedding<ORBIT>(d) ;
	unsigned int eE = getEmbedding<ORBIT>(e) ;
399
400
401
402
//	if(eE != EMBNULL)	// if the source is NULL, nothing to copy
//	{
//		if(dE == EMBNULL)	// if the dest is NULL, create a new cell
//			dE = setOrbitEmbeddingOnNewCell<ORBIT>(d) ;
403
		m_attribs[ORBIT].copyLine(dE, eE) ;	// copy the data
404
//	}
Pierre Kraemer's avatar
Pierre Kraemer committed
405
406
}

407
408
template <unsigned int ORBIT>
inline void GenericMap::copyCell(unsigned int i, unsigned int j)
Pierre Kraemer's avatar
Pierre Kraemer committed
409
{
410
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
411
	m_attribs[ORBIT].copyLine(i, j) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
412
413
}

414
415
template <unsigned int ORBIT>
inline void GenericMap::initCell(unsigned int i)
Pierre Kraemer's avatar
Pierre Kraemer committed
416
{
417
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded");
418
	m_attribs[ORBIT].initLine(i) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
419
420
}

Pierre Kraemer's avatar
Pierre Kraemer committed
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
/****************************************
 *     QUICK TRAVERSAL MANAGEMENT       *
 ****************************************/

template <unsigned int ORBIT>
inline void GenericMap::enableQuickTraversal()
{
	if(m_quickTraversal[ORBIT] == NULL)
	{
		if(!isOrbitEmbedded<ORBIT>())
			addEmbedding<ORBIT>() ;
		m_quickTraversal[ORBIT] = m_attribs[ORBIT].addAttribute<Dart>("quick_traversal") ;
	}
	updateQuickTraversal<ORBIT>() ;
}

template <unsigned int ORBIT>
438
inline void GenericMap::updateQuickTraversal(const FunctorSelect& good)
Pierre Kraemer's avatar
Pierre Kraemer committed
439
{
Pierre Kraemer's avatar
Pierre Kraemer committed
440
441
442
	assert(m_quickTraversal[ORBIT] != NULL || !"updateQuickTraversal on a disabled orbit") ;

	CellMarker<ORBIT> cm(*this) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
443
444
	for(Dart d = begin(); d != end(); next(d))
	{
445
		if(good(d) && !cm.isMarked(d))
Pierre Kraemer's avatar
Pierre Kraemer committed
446
447
		{
			cm.mark(d) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
448
			(*m_quickTraversal[ORBIT])[getEmbedding<ORBIT>(d)] = d ;
Pierre Kraemer's avatar
Pierre Kraemer committed
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
		}
	}
}

template <unsigned int ORBIT>
inline AttributeMultiVector<Dart>* GenericMap::getQuickTraversal()
{
	return m_quickTraversal[ORBIT] ;
}

template <unsigned int ORBIT>
inline void GenericMap::disableQuickTraversal()
{
	if(m_quickTraversal[ORBIT] != NULL)
	{
		m_attribs[ORBIT].removeAttribute<Dart>(m_quickTraversal[ORBIT]->getIndex()) ;
		m_quickTraversal[ORBIT] = NULL ;
	}
}

Pierre Kraemer's avatar
Pierre Kraemer committed
469
470
471
472
/****************************************
 *        ATTRIBUTES MANAGEMENT         *
 ****************************************/

473
474
template <unsigned int ORBIT>
inline AttributeContainer& GenericMap::getAttributeContainer()
Pierre Kraemer's avatar
Pierre Kraemer committed
475
{
476
	return m_attribs[ORBIT] ;
Pierre Kraemer's avatar
Pierre Kraemer committed
477
478
}

479
480
template <unsigned int ORBIT>
inline AttributeMultiVector<Mark>* GenericMap::getMarkVector(unsigned int thread)
Pierre Kraemer's avatar
Pierre Kraemer committed
481
{
482
	assert(isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit not embedded") ;
483
	return m_markTables[ORBIT][thread] ;
Pierre Kraemer's avatar
Pierre Kraemer committed
484
485
}

486
487
template <unsigned int ORBIT>
inline AttributeMultiVector<unsigned int>* GenericMap::getEmbeddingAttributeVector()
Pierre Kraemer's avatar
Pierre Kraemer committed
488
{
489
	return m_embeddings[ORBIT] ;
Pierre Kraemer's avatar
Pierre Kraemer committed
490
491
}

492
493
494
495
496
497
498
499
500
501
502
inline AttributeContainer& GenericMap::getMRAttributeContainer()
{
	return m_mrattribs ;
}

inline AttributeMultiVector<unsigned int>* GenericMap::getMRDartAttributeVector(unsigned int level)
{
	assert(level <= getMaxLevel() || !"Invalid parameter: level does not exist");
	return m_mrDarts[level] ;
}

untereiner's avatar
untereiner committed
503
504
505
506
507
inline AttributeMultiVector<unsigned int>* GenericMap::getMRLevelAttributeVector()
{
	return m_mrLevels ;
}

508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
template <typename R>
bool GenericMap::registerAttribute(const std::string &nameType)
{
	RegisteredBaseAttribute* ra = new RegisteredAttribute<R>;
	if (ra == NULL)
	{
		CGoGNerr << "Erreur enregistrement attribut" << CGoGNendl;
		return false;
	}

	ra->setTypeName(nameType);

	m_attributes_registry_map->insert(std::pair<std::string, RegisteredBaseAttribute*>(nameType,ra));
	return true;
}

524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
/****************************************
 *   EMBEDDING ATTRIBUTES MANAGEMENT    *
 ****************************************/

template <unsigned int ORBIT>
void GenericMap::addEmbedding()
{
	assert(!isOrbitEmbedded<ORBIT>() || !"Invalid parameter: orbit already embedded") ;

	std::ostringstream oss;
	oss << "EMB_" << ORBIT;

	AttributeContainer& dartCont = m_attribs[DART] ;
	AttributeMultiVector<unsigned int>* amv = dartCont.addAttribute<unsigned int>(oss.str()) ;
	m_embeddings[ORBIT] = amv ;

540
	FunctorInitEmb<GenericMap, ORBIT> fsetemb(*this);
541
542
	DartMarker dm(*this);
	for(Dart d = this->begin(); d != this->end(); this->next(d))
543
	{
544
545
546
547
548
549
550
		if(!dm.isMarked(d))
		{
			dm.markOrbit<ORBIT>(d);
			unsigned int em = newCell<ORBIT>();
			fsetemb.changeEmb(em);
			foreach_dart_of_orbit<ORBIT>(d, fsetemb);
		}
551
	}
552
553
}

Pierre Kraemer's avatar
Pierre Kraemer committed
554
555
556
557
/****************************************
 *           DARTS TRAVERSALS           *
 ****************************************/

558
inline Dart GenericMap::begin() const
Pierre Kraemer's avatar
Pierre Kraemer committed
559
{
560
561
562
563
564
565
566
567
568
	if (m_isMultiRes)
	{
		unsigned int d = m_mrattribs.begin() ;
		while (getDartLevel(d) > m_mrCurrentLevel)
			m_mrattribs.next(d) ;
		return Dart::create(d) ;
	}

	return Dart::create(m_attribs[DART].begin()) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
569
570
}

571
inline Dart GenericMap::end() const
Pierre Kraemer's avatar
Pierre Kraemer committed
572
{
573
574
575
	if (m_isMultiRes)
		return Dart::create(m_mrattribs.end()) ;

576
	return Dart::create(m_attribs[DART].end()) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
577
578
}

579
inline void GenericMap::next(Dart& d) const
Pierre Kraemer's avatar
Pierre Kraemer committed
580
{
581
582
	if (m_isMultiRes)
	{
Pierre Kraemer's avatar
Pierre Kraemer committed
583
584
		do
		{
585
			m_mrattribs.next(d.index) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
586
		} while (d.index != m_mrattribs.end() && getDartLevel(d) > m_mrCurrentLevel) ;
587
588
589
	}
	else
		m_attribs[DART].next(d.index) ;
Pierre Kraemer's avatar
Pierre Kraemer committed
590
591
}

592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
template <unsigned int ORBIT>
bool GenericMap::foreach_dart_of_orbit(Dart d, FunctorType& f, unsigned int thread)
{
	switch(ORBIT)
	{
		case DART:		return f(d);
		case VERTEX: 	return foreach_dart_of_vertex(d, f, thread);
		case EDGE: 		return foreach_dart_of_edge(d, f, thread);
		case FACE: 		return foreach_dart_of_face(d, f, thread);
		case VOLUME: 	return foreach_dart_of_volume(d, f, thread);
		case VERTEX1: 	return foreach_dart_of_vertex1(d, f, thread);
		case EDGE1: 	return foreach_dart_of_edge1(d, f, thread);
		case VERTEX2: 	return foreach_dart_of_vertex2(d, f, thread);
		case EDGE2:		return foreach_dart_of_edge2(d, f, thread);
		case FACE2:		return foreach_dart_of_face2(d, f, thread);
Pierre Kraemer's avatar
Pierre Kraemer committed
607
		default: 		assert(!"Cells of this dimension are not handled"); break;
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
	}
	return false;
}

template <unsigned int ORBIT>
bool GenericMap::foreach_orbit(FunctorType& fonct, const FunctorSelect& good, unsigned int thread)
{
	TraversorCell<GenericMap, ORBIT> trav(*this, good, true, thread);
	bool found = false;

	for (Dart d = trav.begin(); !found && d != trav.end(); d = trav.next())
	{
		if ((fonct)(d))
			found = true;
	}
	return found;
}

template <unsigned int ORBIT>
unsigned int GenericMap::getNbOrbits(const FunctorSelect& good)
{
	FunctorCount fcount;
	foreach_orbit<ORBIT>(fcount, good);
	return fcount.getNb();
}

634
template <typename MAP, unsigned int ORBIT, unsigned int INCIDENT>
635
636
unsigned int GenericMap::degree(Dart d)
{
637
	assert(ORBIT != INCIDENT || !"degree does not manage adjacency counting") ;
Pierre Kraemer's avatar
Pierre Kraemer committed
638
	Traversor<MAP>* t = TraversorFactory<MAP>::createIncident(*(reinterpret_cast<MAP*>(this)), d, dimension(), ORBIT, INCIDENT) ;
639
640
	FunctorCount fcount ;
	t->applyFunctor(fcount) ;
641
	delete t ;
642
643
644
	return fcount.getNb() ;
}

645
646
647
648
649
650
651
652
653
654
655
656
657
658
template <unsigned int ORBIT>
void GenericMap::boundaryMarkOrbit(Dart d)
{
	FunctorMark<GenericMap> fm(*this, m_boundaryMarker, m_markTables[DART][0]) ;
	foreach_dart_of_orbit<ORBIT>(d, fm, 0) ;
}

template <unsigned int ORBIT>
void GenericMap::boundaryUnmarkOrbit(Dart d)
{
	FunctorUnmark<GenericMap> fm(*this, m_boundaryMarker, m_markTables[DART][0]) ;
	foreach_dart_of_orbit<ORBIT>(d, fm, 0) ;
}

Sylvain Thery's avatar
Sylvain Thery committed
659
660
661
662
663
664
665
666
667
668
669
/****************************************
 *  TOPOLOGICAL ATTRIBUTES MANAGEMENT   *
 ****************************************/

inline AttributeMultiVector<Dart>* GenericMap::addRelation(const std::string& name)
{
	AttributeContainer& cont = m_attribs[DART] ;
	AttributeMultiVector<Dart>* amv = cont.addAttribute<Dart>(name) ;

	// set new relation to fix point for all the darts of the map
	for(unsigned int i = cont.begin(); i < cont.end(); cont.next(i))
Pierre Kraemer's avatar
Pierre Kraemer committed
670
		(*amv)[i] = i ;
Sylvain Thery's avatar
Sylvain Thery committed
671
672
673
674

	return amv ;
}

Sylvain Thery's avatar
Sylvain Thery committed
675
676
677
678
679
680
681
inline AttributeMultiVector<Dart>* GenericMap::getRelation(const std::string& name)
{
	AttributeContainer& cont = m_attribs[DART] ;
	AttributeMultiVector<Dart>* amv = cont.getDataVector<Dart>(cont.getAttributeIndex(name)) ;
	return amv ;
}

Sylvain Thery's avatar
Sylvain Thery committed
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701

/**************************
 *  BOUNDARY MANAGEMENT   *
 **************************/

inline void GenericMap::boundaryMark(Dart d)
{
	m_markTables[DART][0]->operator[](dartIndex(d)).setMark(m_boundaryMarker);
}

inline void GenericMap::boundaryUnmark(Dart d)
{
	m_markTables[DART][0]->operator[](dartIndex(d)).unsetMark(m_boundaryMarker);
}

inline bool GenericMap::isBoundaryMarked(Dart d) const
{
	return m_markTables[DART][0]->operator[](dartIndex(d)).testMark(m_boundaryMarker);
}

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