mesh_local_smoothing_process.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics 
//
//  License:         BSD License 
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Pooyan Dadvand
//                    
//
               
#if !defined(KRATOS_MESH_LOCAL_SMOOTHING_PROCESS_H_INCLUDED )
#define  KRATOS_MESH_LOCAL_SMOOTHING_PROCESS_H_INCLUDED
 
 
 
// System includes
#include <string>
#include <iostream> 
 
 
// External includes 
 
 
// Project includes
#include "includes/define.h"
#include "processes/process.h"
#include "includes/model_part.h"
 
 
 
namespace Kratos
{
 
  
 
  class KRATOS_API(KRATOS_CORE) MeshLocalSmoothingProcess : public Process
    {
    public:
      
      KRATOS_CLASS_POINTER_DEFINITION(MeshLocalSmoothingProcess);
 
      typedef Node NodeType;
 
      typedef GlobalPointersVector< Node > NeighboursVectorType;
 
      typedef std::vector<Point > PointsVectorType;
  
 
 
      MeshLocalSmoothingProcess(
          ModelPart &rModelPart,
          double AptQuality = 0.5,
          std::size_t MaxIterationsNumber = 10,
          const Flags &rBoundaryFlag = BOUNDARY);
 
      ~MeshLocalSmoothingProcess() override;
      
 
      
      
      
      void Execute() override;
      
      
      
      
      
 
      std::string Info() const override;
      
      void PrintInfo(std::ostream& rOStream) const override;
 
      void PrintData(std::ostream& rOStream) const override;
      
            
      
            
      
    protected:
        
        
        
        
        
        
 
        virtual void FindOptimumPositionsAndWeights(NodeType& rNode, PointsVectorType& rOptimumPoints, Vector& rWeights);
 
        
        
        
        
        
      
            
      
    private:
 
        ModelPart& mrModelPart;
 
        std::size_t mMaxIterationsNumber;
 
        double mAptQuality;
 
        std::size_t mNumberOfLowQualityElements;
 
        double mMeshMinQuality;
 
        double mMeshQualityNorm;
 
        const Flags &mrBoundaryFlag;
        
        
        
        
        
        
 
        void SelectLowQualityElementNodes();
 
        void PerformSmoothing();
 
        void InterpolateNodeOptimumPosition(PointsVectorType const& rOptimumPoints, Vector const& rWeights, Point& OptimumPosition);
 
        void MoveNodeIfImprovesMinimumQuality(NodeType& rNode, Point const& OptimumPosition);
 
 
        
        
        
        
      
      MeshLocalSmoothingProcess& operator=(MeshLocalSmoothingProcess const& rOther);
 
      MeshLocalSmoothingProcess(MeshLocalSmoothingProcess const& rOther);
 
        
        
    }; // Class MeshLocalSmoothingProcess 
 
  
  
  
        
 
  inline std::istream& operator >> (std::istream& rIStream, 
                    MeshLocalSmoothingProcess& rThis);
 
  inline std::ostream& operator << (std::ostream& rOStream, 
                    const MeshLocalSmoothingProcess& rThis)
    {
      rThis.PrintInfo(rOStream);
      rOStream << std::endl;
      rThis.PrintData(rOStream);
 
      return rOStream;
    }
 
 
}  // namespace Kratos.
 
#endif // KRATOS_MESH_LOCAL_SMOOTHING_PROCESS_H_INCLUDED  defined