setboundarycode.cpp

   1 //
   2 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
   3 // +                                                                      +
   4 // + This file is part of enGrid.                                         +
   5 // +                                                                      +
   6 // + Copyright 2008,2009 Oliver Gloth                                     +
   7 // +                                                                      +
   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.                                  +
  12 // +                                                                      +
  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.                         +
  17 // +                                                                      +
  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/>.       +
  20 // +                                                                      +
  21 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  22 //
  23 #include "setboundarycode.h"
  24 #include <vtkIntArray.h>
  25 #include <vtkCellData.h>
  26
  27 SetBoundaryCode::SetBoundaryCode()
  28 {
  29   feature_angle = 180.0;
  30   boundary_code = 1;
  31   setSurfaceIteration();
  32 };
  33
  34 void SetBoundaryCode::pass1()
  35 {
  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;
  44     } else {
  45       pair[i].terminate = false;
  46     };
  47   };
  48 };
  49
  50 void SetBoundaryCode::pass2()
  51 {
  52   EG_VTKDCC(vtkIntArray, cell_code, grid, "cell_code");
  53   vtkIdType cellId;
  54   foreach(cellId, item) {
  55     cell_code->SetValue(cellId, boundary_code);
  56   };
  57 };
  58