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