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

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

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

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

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

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

44

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

47
int GenericMap::m_nbInstances = 0;
48

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

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

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


56

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

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

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

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

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

Pierre Kraemer's avatar
Pierre Kraemer committed
81 82 83
		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
84

Pierre Kraemer's avatar
Pierre Kraemer committed
85 86 87
		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
88

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

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

		registerAttribute<MarkerBool>("MarkerBool");
Pierre Kraemer's avatar
Pierre Kraemer committed
96 97
	}

Sylvain Thery's avatar
Sylvain Thery committed
98 99


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

104
	s_instances->push_back(this);
105

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


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

Pierre Kraemer's avatar
Pierre Kraemer committed
125
	init();
Pierre Kraemer's avatar
Pierre Kraemer committed
126 127 128 129
}

GenericMap::~GenericMap()
{
130

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

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

141 142 143

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

149
		delete m_attributes_registry_map;
150
		m_attributes_registry_map = NULL;
Sylvain Thery's avatar
Sylvain Thery committed
151 152 153 154 155 156 157 158 159 160

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


161
	}
162

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

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

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

		for(unsigned int j = 0; j < NB_THREAD; ++j)
			m_markVectors_free[i][j].clear();
Pierre Kraemer's avatar
Pierre Kraemer committed
185 186
	}

187

188 189 190 191 192 193
	if (addBoundaryMarkers)
	{
		m_boundaryMarkers[0] = m_attribs[DART].addAttribute<MarkerBool>("BoundaryMark0") ;
		m_boundaryMarkers[1] = m_attribs[DART].addAttribute<MarkerBool>("BoundaryMark1") ;
	}

Pierre Kraemer's avatar
Pierre Kraemer committed
194

Sylvain Thery's avatar
Sylvain Thery committed
195
	 for(std::multimap<AttributeMultiVectorGen*, AttributeHandlerGen*>::iterator it = attributeHandlers.begin(); it != attributeHandlers.end(); ++it)
Pierre Kraemer's avatar
Pierre Kraemer committed
196 197 198 199 200 201 202 203
		(*it).second->setInvalid() ;
	attributeHandlers.clear() ;
}

void GenericMap::clear(bool removeAttrib)
{
	if (removeAttrib)
		init();
Sylvain Thery's avatar
Sylvain Thery committed
204 205 206 207 208
	else
	{
		for(unsigned int i = 0; i < NB_ORBITS; ++i)
			m_attribs[i].clear(false) ;
	}
untereiner's avatar
untereiner committed
209 210
}

Pierre Kraemer's avatar
Pierre Kraemer committed
211 212 213 214
/****************************************
 *        ATTRIBUTES MANAGEMENT         *
 ****************************************/

untereiner's avatar
untereiner committed
215 216 217 218 219 220 221 222 223 224 225 226
void GenericMap::swapEmbeddingContainers(unsigned int orbit1, unsigned int orbit2)
{
	assert(orbit1 != orbit2 || !"Cannot swap a container with itself") ;
	assert((orbit1 != DART && orbit2 != DART) || !"Cannot swap the darts container") ;

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

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


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

Sylvain Thery's avatar
Sylvain Thery committed
229 230 231 232 233 234 235 236 237 238
//	for (unsigned int i = 0; i < NB_THREAD; ++i)
//	{
//		for(std::vector<CellMarkerGen*>::iterator it = cellMarkers[i].begin(); it != cellMarkers[i].end(); ++it)
//		{
//			if((*it)->m_cell == orbit1)
//				(*it)->m_cell = orbit2 ;
//			else if((*it)->m_cell == orbit2)
//				(*it)->m_cell = orbit1 ;
//		}
//	}
untereiner's avatar
untereiner committed
239
}
Pierre Kraemer's avatar
Pierre Kraemer committed
240

Pierre Kraemer's avatar
Pierre Kraemer committed
241 242 243 244 245 246 247 248 249 250 251 252
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)
253 254 255 256 257
		{
			unsigned int id = cont.getAttributeIndex(*it);
			std::cout << "    " << *it << " ("<<id<<")"<<std::endl ;
		}

Pierre Kraemer's avatar
Pierre Kraemer committed
258 259 260 261 262 263 264
		std::cout << "-------------------------" << std::endl ;
	}
	std::cout << "m_embeddings: " << std::hex ;
	for (unsigned int i = 0; i < NB_ORBITS; ++i)
		std::cout << (long)(m_embeddings[i]) << " / " ;
	std::cout << std::endl << "-------------------------" << std::endl ;

Sylvain Thery's avatar
Sylvain Thery committed
265 266 267 268
//	std::cout << "m_markTables: " ;
//	for (unsigned int i = 0; i < NB_ORBITS; ++i)
//		std::cout << (long)(m_markTables[i][0]) << " / " ;
//	std::cout << std::endl << "-------------------------" << std::endl << std::dec ;
Pierre Kraemer's avatar
Pierre Kraemer committed
269 270
}

271 272 273 274 275 276 277
void GenericMap::printDartsTable()
{
	std::cout << "======================="<< std::endl ;

	//m_attribs[DART]
}

278 279 280 281
/****************************************
 *          THREAD MANAGEMENT           *
 ****************************************/

Sylvain Thery's avatar
Sylvain Thery committed
282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323
//void GenericMap::addThreadMarker(unsigned int nb)
//{
//	unsigned int th ;

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

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

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

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

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

Pierre Kraemer's avatar
Pierre Kraemer committed
329
void GenericMap::restore_shortcuts()
Sylvain Thery's avatar
Sylvain Thery committed
330
{
Pierre Kraemer's avatar
Pierre Kraemer committed
331 332
	// EMBEDDING

Sylvain Thery's avatar
Sylvain Thery committed
333 334 335 336 337 338 339 340 341 342
	// 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
343
		std::string sub = listeNames[i].substr(0, listeNames[i].size() - 1);
Sylvain Thery's avatar
Sylvain Thery committed
344 345 346 347 348 349 350 351
		if (sub == "EMB_")
		{
			unsigned int orb = listeNames[i][4]-'0'; // easy atoi computation for one char;
			AttributeMultiVector<unsigned int>* amv = cont.getDataVector<unsigned int>(i);
			m_embeddings[orb] = amv ;
		}
	}

Sylvain Thery's avatar
Sylvain Thery committed
352 353
	// MARKERS
	m_attribs[DART].getAttributesNames(listeNames);
Pierre Kraemer's avatar
Pierre Kraemer committed
354

Sylvain Thery's avatar
Sylvain Thery committed
355
	for (unsigned int i = 0;  i < listeNames.size(); ++i)
Sylvain Thery's avatar
Sylvain Thery committed
356
	{
Sylvain Thery's avatar
Sylvain Thery committed
357 358 359 360 361 362
		if (listeNames[i] == "BoundaryMark0")
			m_boundaryMarkers[0] = cont.getDataVector<MarkerBool>(i);

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

Sylvain Thery's avatar
Sylvain Thery committed
364
	// QUICK TRAVERSAL
Pierre Kraemer's avatar
Pierre Kraemer committed
365

Sylvain Thery's avatar
Sylvain Thery committed
366 367 368
	for(unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
	{
		AttributeContainer& cont = m_attribs[orbit];
Pierre Kraemer's avatar
Pierre Kraemer committed
369 370 371 372 373 374 375
		m_quickTraversal[orbit] = cont.getDataVector<Dart>("quick_traversal") ;
		for(unsigned int j = 0; j < NB_ORBITS; ++j)
		{
			std::stringstream ss;
			ss << "quickLocalIncidentTraversal_" << j;
			m_quickLocalIncidentTraversal[orbit][j] = cont.getDataVector< NoTypeNameAttribute<std::vector<Dart> > >(ss.str()) ;
			std::stringstream ss2;
Pierre Kraemer's avatar
Pierre Kraemer committed
376
			ss2 << "quickLocalAdjacentTraversal_" << j;
Pierre Kraemer's avatar
Pierre Kraemer committed
377 378
			m_quickLocalAdjacentTraversal[orbit][j] = cont.getDataVector< NoTypeNameAttribute<std::vector<Dart> > >(ss2.str()) ;
		}
379
	}
380

381 382 383 384
	// set Attribute handlers invalid
	for(std::multimap<AttributeMultiVectorGen*, AttributeHandlerGen*>::iterator it = attributeHandlers.begin(); it != attributeHandlers.end(); ++it)
		(*it).second->setInvalid() ;
	attributeHandlers.clear() ;
Sylvain Thery's avatar
Sylvain Thery committed
385 386 387 388
}

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

Sylvain Thery's avatar
Sylvain Thery committed
427 428

void GenericMap::compact(bool topoOnly)
429
{
Sylvain Thery's avatar
Sylvain Thery committed
430 431 432 433 434 435
	compactTopo();

	if (topoOnly)
		return;

	std::vector<unsigned int> oldnew;
Sylvain Thery's avatar
Sylvain Thery committed
436

437 438 439 440
	for (unsigned int orbit = 0; orbit < NB_ORBITS; ++orbit)
	{
		if ((orbit != DART) && (isOrbitEmbedded(orbit)))
		{
Sylvain Thery's avatar
Sylvain Thery committed
441 442
			m_attribs[orbit].compact(oldnew);
			for (unsigned int i = m_attribs[DART].begin(); i != m_attribs[DART].end(); m_attribs[DART].next(i))
443
			{
444
				unsigned int& idx = m_embeddings[orbit]->operator[](i);
Sylvain Thery's avatar
Sylvain Thery committed
445 446
				unsigned int jdx = oldnew[idx];
				if (jdx != 0xffffffff)
447 448 449 450 451 452
					idx = jdx;
			}
		}
	}
}

Sylvain Thery's avatar
Sylvain Thery committed
453

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

Sylvain Thery's avatar
Sylvain Thery committed
463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543

void GenericMap::moveData(GenericMap &mapf)
{
	GenericMap::init(false);

	for (unsigned int i=0; i<NB_ORBITS; ++i)
	{
		this->m_attribs[i].swap(mapf.m_attribs[i]);
		this->m_embeddings[i] = mapf.m_embeddings[i];
		this->m_quickTraversal[i] = mapf.m_quickTraversal[i];
		mapf.m_embeddings[i] = NULL ;
		mapf.m_quickTraversal[i] = NULL;

		for (unsigned int j=0; j<NB_ORBITS; ++j)
		{
			this->m_quickLocalIncidentTraversal[i][j] = mapf.m_quickLocalIncidentTraversal[i][j];
			this->m_quickLocalAdjacentTraversal[i][j] = mapf.m_quickLocalAdjacentTraversal[i][j];
			mapf.m_quickLocalIncidentTraversal[i][j] = NULL ;
			mapf.m_quickLocalAdjacentTraversal[i][j] = NULL ;
		}

		for (unsigned int j=0; j<NB_THREAD; ++j)
			this->m_markVectors_free[i][j].swap(mapf.m_markVectors_free[i][j]);
	}

	this->m_boundaryMarkers[0] = mapf.m_boundaryMarkers[0];
	this->m_boundaryMarkers[1] = mapf.m_boundaryMarkers[1];

	this->m_nextMarkerId = mapf.m_nextMarkerId;

	for(auto it = mapf.attributeHandlers.begin(); it != mapf.attributeHandlers.end(); ++it)
	   (*it).second->setInvalid() ;
	mapf.attributeHandlers.clear() ;

}


void GenericMap::garbageMarkVectors()
{
	unsigned int maxId=0;

	for (unsigned int orbit=0; orbit<NB_ORBITS;++orbit)
	{
		std::vector<std::string> attNames;
		m_attribs[orbit].getAttributesNames(attNames);
		for (auto sit=attNames.begin(); sit!=attNames.end();++sit)
		{
			if (sit->substr(0,7) == "marker_")
			{
				std::string num = sit->substr(sit->length()-3,3);
				unsigned int id = 100*(num[0]-'0')+10*(num[1]-'0')+(num[2]-'0');
				if (id > maxId)
					maxId = id;
				AttributeMultiVector<MarkerBool>* amv = m_attribs[orbit].getDataVector<MarkerBool>(*sit);
				amv->allFalse();
				m_markVectors_free[orbit][0].push_back(amv);
			}
		}
	}
	m_nextMarkerId = maxId+1;
}

void GenericMap::removeMarkVectors()
{
	for (unsigned int orbit=0; orbit<NB_ORBITS;++orbit)
	{
		std::vector<std::string> attNames;
		m_attribs[orbit].getAttributesNames(attNames);
		for (auto sit=attNames.begin(); sit!=attNames.end();++sit)
		{
			if (sit->substr(0,7) == "marker_")
			{
				m_attribs[orbit].removeAttribute<MarkerBool>(*sit);
			}
		}
	}
	m_nextMarkerId = 0;
}



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