src/libengrid/deletevolumegrid.cpp

   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 #include "deletevolumegrid.h"
  22
  23 void DeleteVolumeGrid::operate()
  24 {
  25   EG_VTKSP(vtkUnstructuredGrid, sgrid);
  26   QVector<vtkIdType> scells, snodes;
  27   QVector<int>       _snodes;
  28   getAllSurfaceCells(scells, m_Grid);
  29   getNodesFromCells(scells, snodes, m_Grid);
  30   createNodeMapping(snodes, _snodes, m_Grid);
  31   allocateGrid(sgrid, scells.size(), snodes.size());
  32   {
  33     vtkIdType newId = 0;
  34     foreach (vtkIdType nodeId, snodes) {
  35       vec3_t x;
  36       m_Grid->GetPoint(nodeId, x.data());
  37       sgrid ->GetPoints()->SetPoint(newId, x.data());
  38       copyNodeData(m_Grid, nodeId, sgrid, newId);
  39       ++newId;
  40     }
  41   }
  42   foreach (vtkIdType cellId, scells) {
  43     vtkIdType *pts, Npts;
  44     m_Grid->GetCellPoints(cellId, Npts, pts);
  45     for (int i = 0; i < Npts; ++i) {
  46       pts[i] = _snodes[pts[i]];
  47     }
  48     vtkIdType cellType = m_Grid->GetCellType(cellId);
  49     vtkIdType newId = sgrid->InsertNextCell(cellType, Npts, pts);
  50     copyCellData(m_Grid, cellId, sgrid, newId);
  51   }
  52   makeCopy(sgrid, m_Grid);
  53 }
  54