updated to modern VTK
[engrid-github.git] / src / libengrid / facefinder.h
blob680092cc11b0f780c7303d1ed8ac01673d501eb4
1 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2 // + +
3 // + This file is part of enGrid. +
4 // + +
5 // + Copyright 2008-2014 enGits GmbH +
6 // + +
7 // + enGrid is free software: you can redistribute it and/or modify +
8 // + it under the terms of the GNU General Public License as published by +
9 // + the Free Software Foundation, either version 3 of the License, or +
10 // + (at your option) any later version. +
11 // + +
12 // + enGrid is distributed in the hope that it will be useful, +
13 // + but WITHOUT ANY WARRANTY; without even the implied warranty of +
14 // + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +
15 // + GNU General Public License for more details. +
16 // + +
17 // + You should have received a copy of the GNU General Public License +
18 // + along with enGrid. If not, see <http://www.gnu.org/licenses/>. +
19 // + +
20 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
21 #ifndef FACEFINDER_H
22 #define FACEFINDER_H
24 #include "octree.h"
25 #include "triangle.h"
26 #include "timer.h"
28 #include <QVector>
29 #include <QList>
31 class FaceFinder : public EgVtkObject
34 Octree m_Octree;
35 vtkUnstructuredGrid* m_Grid;
36 QVector<QList<vtkIdType> > m_Faces;
37 double m_MinSize;
38 int m_MaxFaces;
39 QVector<Triangle> m_Triangles;
40 QVector<vec3_t> m_Centres;
41 QVector<double> m_CritLengthOctreeCell;
42 QVector<double> m_CritLengthNode;
43 Timer m_Timer;
44 QVector<vec3_t> m_CollectedPoints;
45 int m_NumCollectedPoints;
46 int m_CollectedPointsIncrement;
47 bool m_UseImprovedFaceSearch;
50 private: // methods
52 void getPointsOfFace(vtkIdType id_face);
53 double calcCritLength(vtkIdType id_cell);
54 void calcCritLengthForAllNodes();
55 int refine();
58 public: // methods
60 FaceFinder();
62 void setGrid(vtkUnstructuredGrid *grid);
63 void setMaxNumFaces(int N) { m_MaxFaces = N; }
64 void getCloseFaces(vec3_t x, QVector<vtkIdType> &faces);
65 vtkIdType getClosestFace(vec3_t x, double &L);
69 #endif // FACEFINDER_H