1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
|
/*****************************************************************************
* $CAMITK_LICENCE_BEGIN$
*
* CamiTK - Computer Assisted Medical Intervention ToolKit
* (c) 2001-2025 Univ. Grenoble Alpes, CNRS, Grenoble INP - UGA, TIMC, 38000 Grenoble, France
*
* Visit http://camitk.imag.fr for more information
*
* This file is part of CamiTK.
*
* CamiTK is free software: you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License version 3
* only, as published by the Free Software Foundation.
*
* CamiTK 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 version 3 for more details.
*
* You should have received a copy of the GNU Lesser General Public License
* version 3 along with CamiTK. If not, see <http://www.gnu.org/licenses/>.
*
* $CAMITK_LICENCE_END$
****************************************************************************/
#ifndef PMLTRANSFORM_H
#define PMLTRANSFORM_H
#include "PhysicalModel.h"
#include <map>
#include <vector>
class CellProperties;
/**
* @ingroup group_cepmodeling_libraries_pml
*
* @brief
* Class facet to old and compare facet
*
**/
class Facet {
public:
/// create a facet using size nodes and their indexes
Facet(unsigned int size, unsigned int id[]);
/// destructor
virtual ~Facet();
/// if it is the same (equivalent) facet, increment used (return true if equivalence)
bool testEquivalence(unsigned int size, unsigned int id[]);
/// return the corresponding PML cell
Cell* getCell(PhysicalModel*) const;
/// print on stdout
void debug();
/// get the number of time it is being used
unsigned int getUsed() const;
private:
/// is this atom index present in this facet (no check on the order)
bool isIn(unsigned int) const;
/// the facet atom indexes
unsigned int* id;
/// nr of atoms composing the facet (3 = triangle, 4 = quad)
unsigned int size;
/// nr of times the facet is used
unsigned int used;
};
/** PML Transform is composed by static methods
* It performs transformations on pml object and do a lot of useful things
*/
class PMLTransform {
public :
//-- elem to neighborhhod methods
/// get the iterator on the correct atom index in the neighMap
/// if non existant create it
static std::map<unsigned int, Cell*>::iterator getIterator(unsigned int index);
/// generate the neighborhoods
static StructuralComponent* generateNeighborhood(StructuralComponent* sc);
/// check if equivalent of already existing facet
static void equivalent(int size, unsigned int id[]);
/// generate the outside surface
static MultiComponent* generateExternalSurface(StructuralComponent* sc);
private :
// -------------------- Neigborhood Map ------------------------
// associative map of all the neighboors for a given index of an atom
static std::map<unsigned int, Cell*> neighMap;
// -------------------- All border facets ------------------------
/// storing all the facets
static std::vector <Facet*> allFacets;
};
#endif
|
