MeshChecker.hpp

/*---------------------------------------------------------------------------*\
 *
 *  mimmo
 *
 *  Copyright (C) 2015-2021 OPTIMAD engineering Srl
 *
 *  -------------------------------------------------------------------------
 *  License
 *  This file is part of mimmo.
 *
 *  mimmo is free software: you can redistribute it and/or modify it
 *  under the terms of the GNU Lesser General Public License v3 (LGPL)
 *  as published by the Free Software Foundation.
 *
 *  mimmo is distributed in the hope that it will be useful, but WITHOUT
 *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
 *  FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
 *  License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with mimmo. If not, see <http://www.gnu.org/licenses/>.
 *
 \ *---------------------------------------------------------------------------*/
 
#ifndef __MESH_CHECKER_HPP__
#define __MESH_CHECKER_HPP__
 
#include <BaseManipulation.hpp>
 
 
namespace mimmo{
 
class MeshChecker: public BaseManipulation {
 
public:
 
    enum CMeshOutput{               NOTRUN = -1,                GOOD = 0,       FACEVALIDITY = 1, VOLUMECHANGERATIO = 2,   BOUNDARYSKEWNESS = 3,            SKEWNESS = 4,      MINIMUMVOLUME = 5,      MAXIMUMVOLUME = 6
    };
 
 
    MeshChecker();
    MeshChecker(const bitpit::Config::Section & rootXML);
    virtual ~MeshChecker();
 
    MeshChecker(const MeshChecker & other);
    MeshChecker & operator=(MeshChecker other);
 
    void setGeometry(MimmoSharedPointer<MimmoObject> obj);
    void setMinimumVolumeTolerance(double tol);
    void setMaximumVolumeTolerance(double tol);
    void setMaximumSkewnessTolerance(double tol);
    void setMaximumBoundarySkewnessTolerance(double tol);
    void setMinimumFaceValidityTolerance(double tol);
    void setMinimumVolumeChangeTolerance(double tol);
    void setPrintResumeFile(bool flag);
 
    bool isGood();
    CMeshOutput getQualityStatus();
    int  getQualityStatusInt();
    void execute();
    void buildPorts();
 
    void plotOptionalResults();
    virtual void absorbSectionXML(const bitpit::Config::Section & slotXML, std::string name="");
    virtual void flushSectionXML(bitpit::Config::Section & slotXML, std::string name="");
 
 
protected:
 
    void swap(MeshChecker & x) noexcept;
    void setDefault();
    CMeshOutput  checkMeshQuality();
    bool checkVolume();
    bool checkSkewness();
    bool checkFaceValidity();
    void initializeVolumes();
    void clear();
    void printResumeFile();
 
protected:
 
    // check values
    double  m_minVolume;            
    double  m_maxVolume;            
    double  m_maxSkewness;          
    double  m_maxSkewnessBoundary;  
    double  m_minFaceValidity;      
    double  m_maxFaceValidity;      
    double  m_minVolumeChange;      
    // tolerance values
    double  m_minVolumeTol;             
    double  m_maxVolumeTol;             
    double  m_maxSkewnessTol;           
    double  m_maxSkewnessBoundaryTol;   
    double  m_minFaceValidityTol;       
    double  m_minVolumeChangeTol;       
    bool        m_isGood;       
    CMeshOutput m_qualityStatus;
    //Temporary aux variables
    MimmoPiercedVector<double>  m_volumes; 
    std::unique_ptr<MimmoObject>    m_volume; 
    std::unique_ptr<MimmoObject>    m_skewness; 
    std::unique_ptr<MimmoObject>    m_facevalidity; 
    std::unique_ptr<MimmoObject>    m_volumechange; 
    bool m_printResumeFile; 
};
 
REGISTER_PORT(M_GEOM, MC_SCALAR, MD_MIMMO_,__MESH_CHECKER_HPP__)
REGISTER_PORT(M_VALUEB, MC_SCALAR, MD_BOOL,__MESH_CHECKER_HPP__)
REGISTER_PORT(M_VALUEI, MC_SCALAR, MD_INT,__MESH_CHECKER_HPP__)
 
 
REGISTER(BaseManipulation, MeshChecker, "mimmo.MeshChecker")
}
 
#endif