2 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
4 // + This file is part of enGrid. +
6 // + Copyright 2008,2009 Oliver Gloth +
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>
27 SetBoundaryCode::SetBoundaryCode()
29 feature_angle
= 180.0;
31 setSurfaceIteration();
34 void SetBoundaryCode::pass1()
36 using namespace GeometryTools
;
37 double fa
= feature_angle
*M_PI
/180.0;
38 for (int i
= 0; i
< pair
.size(); ++i
) {
39 vec3_t n1
= cellNormal(grid
, pair
[i
].item1
);
40 vec3_t n2
= cellNormal(grid
, pair
[i
].item2
);
41 double cosa
= (n1
*n2
)/(n1
.abs()*n2
.abs());
42 if (fabs(acos(cosa
)) > fa
) {
43 pair
[i
].terminate
= true;
45 pair
[i
].terminate
= false;
50 void SetBoundaryCode::pass2()
52 EG_VTKDCC(vtkIntArray
, cell_code
, grid
, "cell_code");
54 foreach(cellId
, item
) {
55 cell_code
->SetValue(cellId
, boundary_code
);