2 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
4 // + This file is part of enGrid. +
6 // + Copyright 2008-2010 enGits GmbH +
8 // + enGrid is free software: you can redistribute it and/or modify +
9 // + it under the terms of the GNU General Public License as published by +
10 // + the Free Software Foundation, either version 3 of the License, or +
11 // + (at your option) any later version. +
13 // + enGrid is distributed in the hope that it will be useful, +
14 // + but WITHOUT ANY WARRANTY; without even the implied warranty of +
15 // + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +
16 // + GNU General Public License for more details. +
18 // + You should have received a copy of the GNU General Public License +
19 // + along with enGrid. If not, see <http://www.gnu.org/licenses/>. +
21 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
23 #include "setboundarycode.h"
24 #include <vtkIntArray.h>
25 #include <vtkCellData.h>
26 #include "guimainwindow.h"
28 SetBoundaryCode::SetBoundaryCode()
30 feature_angle
= 180.0;
32 setSurfaceIteration();
36 void SetBoundaryCode::pass1()
38 if(!(SelectAllVisible
|| OnlyPickedCell
|| OnlyPickedCellAndNeighbours
)) {
39 using namespace GeometryTools
;
40 double fa
= feature_angle
*M_PI
/180.0;
43 GuiMainWindow::pointer()->getDisplayBoundaryCodes(DBC
);
44 DBC
.insert(boundary_code
);
47 EG_VTKDCC(vtkIntArray
, cell_code
, m_Grid
, "cell_code");
49 for (int i
= 0; i
< pair
.size(); ++i
) {
50 int bc1
= cell_code
->GetValue(pair
[i
].item1
);
51 int bc2
= cell_code
->GetValue(pair
[i
].item2
);
53 vec3_t n1
= cellNormal(m_Grid
, pair
[i
].item1
);
54 vec3_t n2
= cellNormal(m_Grid
, pair
[i
].item2
);
55 double cosa
= (n1
*n2
)/(n1
.abs()*n2
.abs());
56 if (fabs(acos(cosa
)) > fa
) {
57 pair
[i
].terminate
= true;
59 pair
[i
].terminate
= false;
62 if(DBC
.contains(bc1
) && DBC
.contains(bc2
)){
63 vec3_t n1
= cellNormal(m_Grid
, pair
[i
].item1
);
64 vec3_t n2
= cellNormal(m_Grid
, pair
[i
].item2
);
65 double cosa
= (n1
*n2
)/(n1
.abs()*n2
.abs());
66 if (fabs(acos(cosa
)) > fa
) {
67 pair
[i
].terminate
= true;
69 pair
[i
].terminate
= false;
72 pair
[i
].terminate
= true;
79 void SetBoundaryCode::pass2()
81 EG_VTKDCC(vtkIntArray
, cell_code
, m_Grid
, "cell_code");
83 if(SelectAllVisible
) {
85 GuiMainWindow::pointer()->getDisplayBoundaryCodes(DBC
);
86 DBC
.insert(boundary_code
);
87 foreach(cellId
, cells
) {
88 int bc
= cell_code
->GetValue(cellId
);
90 cell_code
->SetValue(cellId
, boundary_code
);
93 } else if (OnlyPickedCell
) {
94 cout
<<"this->getStart()="<<this->getStart()<<endl
;
95 cell_code
->SetValue(this->getStart(), boundary_code
);
96 } else if (OnlyPickedCellAndNeighbours
) {
97 cell_code
->SetValue(this->getStart(), boundary_code
);
98 foreach (cellId
, c2c
[this->getStart()]) {
99 cell_code
->SetValue(cellId
, boundary_code
);
102 foreach (cellId
, item
) {
103 cell_code
->SetValue(cellId
, boundary_code
);