src/libengrid/meshquality.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 "meshquality.h"
  22
  23 MeshQuality::MeshQuality()
  24 {
  25   m_Name = "unknown mesh quality";
  26 }
  27
  28 void MeshQuality::computeNodesFromCells()
  29 {
  30   EG_VTKDCN(vtkDoubleArray, node_mesh_quality, m_Grid, "node_mesh_quality");
  31   EG_VTKDCC(vtkDoubleArray, cell_mesh_quality, m_Grid, "cell_mesh_quality");
  32   EG_FORALL_NODES(id_node, m_Grid) {
  33     double mq = 1;
  34     for (int i = 0; i < m_Part.n2cGSize(id_node); ++i) {
  35       vtkIdType id_cell = m_Part.n2cGG(id_node, i);
  36       mq = min(cell_mesh_quality->GetValue(id_cell), mq);
  37     }
  38     node_mesh_quality->SetValue(id_node, mq);
  39   }
  40 }
  41
  42 void MeshQuality::printCellInfo(int indent)
  43 {
  44   double min_quality = 1e99;
  45   double max_quality = 0;
  46   double mean_quality = 0;
  47   EG_VTKDCC(vtkDoubleArray, cell_mesh_quality, m_Grid, "cell_mesh_quality");
  48   EG_FORALL_CELLS(id_cell, m_Grid) {
  49     double q = cell_mesh_quality->GetValue(id_cell);
  50     min_quality = min(q, min_quality);
  51     max_quality = max(q, max_quality);
  52     mean_quality += q;
  53   }
  54   mean_quality /= m_Grid->GetNumberOfCells();
  55   QString indent_str;
  56   indent_str.fill(' ', indent);
  57   cout << qPrintable(indent_str) << qPrintable(name()) << endl;
  58   cout << qPrintable(indent_str) << "min  : " << min_quality << endl;
  59   cout << qPrintable(indent_str) << "mean : " << mean_quality << endl;
  60   cout << qPrintable(indent_str) << "max  : " << max_quality << endl;
  61 }