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