src/libengrid/meshquality.cpp

   1 //
   2 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
   3 // +                                                                      +
   4 // + This file is part of enGrid.                                         +
   5 // +                                                                      +
   6 // + Copyright 2008-2013 enGits GmbH                                      +
   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 "meshquality.h"
  24
  25 MeshQuality::MeshQuality()
  26 {
  27 }
  28
  29 void MeshQuality::computeNodesFromCells()
  30 {
  31   EG_VTKDCN(vtkDoubleArray, node_mesh_quality, m_Grid, "node_mesh_quality");
  32   EG_VTKDCC(vtkDoubleArray, cell_mesh_quality, m_Grid, "cell_mesh_quality");
  33   EG_FORALL_NODES(id_node, m_Grid) {
  34     double mq = 1;
  35     for (int i = 0; i < m_Part.n2cGSize(id_node); ++i) {
  36       vtkIdType id_cell = m_Part.n2cGG(id_node, i);
  37       mq = min(cell_mesh_quality->GetValue(id_cell), mq);
  38     }
  39     node_mesh_quality->SetValue(id_node, mq);
  40   }
  41 }