src/libengrid/updatedesiredmeshdensity.h

   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 #ifndef UPDATEDESIREDMESHDENSITY_H
  24 #define UPDATEDESIREDMESHDENSITY_H
  25
  26 #include "surfaceoperation.h"
  27
  28 #include "vertexmeshdensity.h"
  29 #include "edgelengthsourcemanager.h"
  30
  31
  32 /// Update desired mesh density, i.e. the field used for surface meshing
  33
  34 class UpdateDesiredMeshDensity : public SurfaceOperation
  35 {
  36
  37 private: // data types
  38
  39   struct point_t {
  40     vec3_t x;
  41     vec3_t n;
  42     double L;
  43     QList<int> idx;
  44     vtkIdType id_face;
  45   };
  46
  47 private: //attributes
  48
  49   QSet<int>                   m_BCs;
  50   double                      m_GrowthFactor;
  51   QVector <VertexMeshDensity> m_VMDvector; ///< the mesh density rules
  52   double                      m_MaxEdgeLength;
  53   double                      m_MinEdgeLength;
  54   double                      m_NodesPerQuarterCircle;
  55   double                      m_FeatureResolution2D;
  56   double                      m_FeatureResolution3D;
  57   double                      m_FeatureThresholdAngle;
  58   QVector<double>             m_FeatureSize;
  59   QVector<bool>               m_Fixed;
  60   EdgeLengthSourceManager     m_ELSManager;
  61   bool                        m_OnlySurfaceCells;
  62
  63 protected: // methods
  64
  65   void   computeFeature(const QList<point_t> points, QVector<double> &cl_pre, double res);
  66   void   computeFeature2D(QVector<double> &cl_pre);
  67   void   computeFeature3D(QVector<double> &cl_pre);
  68   double computeSearchDistance(vtkIdType id_face);
  69   void   computeExistingLengths();
  70
  71
  72 public: //methods
  73
  74   UpdateDesiredMeshDensity();
  75   virtual void operate();
  76   void setVertexMeshDensityVector(QVector <VertexMeshDensity> const & vmd) { m_VMDvector = vmd; }
  77   void setMaxEdgeLength(double l) { m_MaxEdgeLength = l; }
  78   void setMinEdgeLength(double l) { m_MinEdgeLength = l; }
  79   void setNodesPerQuarterCircle(double N) { m_NodesPerQuarterCircle = N; }
  80   void setCellGrowthFactor(double cgf) { m_GrowthFactor = cgf; }
  81   void setVolumeCellsOn() { m_OnlySurfaceCells = false; }
  82   void setVolumeCellsOff() { m_OnlySurfaceCells = true; }
  83   void setFeatureResolution2D(double n) { m_FeatureResolution2D = n; }
  84   void setFeatureResolution3D(double n) { m_FeatureResolution3D = n; }
  85   void setFeatureThresholdAngle(double a) { m_FeatureThresholdAngle = a; }
  86
  87 };
  88
  89 #endif