src/libengrid/laplacesmoother.h

   1 //
   2 // ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
   3 // +                                                                      +
   4 // + This file is part of enGrid.                                         +
   5 // +                                                                      +
   6 // + Copyright 2008-2012 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 LAPLACESMOOTHER_H
  24 #define LAPLACESMOOTHER_H
  25
  26 #include "surfaceoperation.h"
  27 #include "surfaceprojection.h"
  28
  29 class LaplaceSmoother : public SurfaceOperation
  30 {
  31
  32 private:
  33
  34   QSet<int> m_BCs;
  35   int       m_NumberOfIterations;
  36   bool      m_UseProjection;
  37   bool      m_UseNormalCorrection;
  38   double    m_UnderRelaxation;
  39   bool      m_Success;
  40   int       m_ProjectionIterations;
  41   double    m_FeatureMagic;
  42   double    m_Limit;
  43
  44   QVector<QVector<int> > m_NodeToBc;
  45
  46   bool      m_CorrectCurvature;
  47   bool      m_NoCheck;
  48
  49   QSet<vtkIdType> m_AllowedCellTypes;
  50   QVector<bool> m_Fixed;
  51
  52
  53 private: // methods
  54
  55   bool setNewPosition(vtkIdType id_node, vec3_t x_new);
  56   bool moveNode(vtkIdType id_node, vec3_t &Dx);
  57   void featureCorrection(vtkIdType id_node, SurfaceProjection *proj, vec3_t &x_new);
  58
  59
  60 public:
  61
  62   LaplaceSmoother(); ///< default constructor
  63   virtual void operate(); ///< Run operation
  64   void setNumberOfIterations(int N) { m_NumberOfIterations = N;} ///< Set number of iterations
  65   void setProjectionOn() { m_UseProjection = true; }
  66   void setProjectionOff() { m_UseProjection = false; }
  67   void setNormalCorrectionOn() { m_UseNormalCorrection = true; }
  68   void setNormalCorrectionOff() { m_UseNormalCorrection = false; }
  69   bool succeeded() { return m_Success; }
  70   void fixNodes(const QVector<bool> &fixnodes);
  71
  72 public:
  73
  74   void setCorrectCurvature(bool b) { m_CorrectCurvature = b; }
  75   bool getCorrectCurvature() { return m_CorrectCurvature; }
  76   void setNoCheck(bool b) { m_NoCheck = b; }
  77   bool getNoCheck() { return m_NoCheck; }
  78   void setProjectionIterations(int n) { m_ProjectionIterations = n; }
  79
  80 };
  81
  82 #endif