RefineGeometry.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 __REFINEGEOMETRY_HPP__
#define __REFINEGEOMETRY_HPP__
 
#include "BaseManipulation.hpp"
 
namespace mimmo{
 
enum class RefineType{
    TERNARY = 0, 
            REDGREEN = 1 
};
 
class RefineGeometry: public BaseManipulation{
 
protected:
    RefineType  m_type;     
    int         m_refinements; 
    int         m_steps;    
    std::vector<long> m_activecells; 
public:
    RefineGeometry();
    RefineGeometry(const bitpit::Config::Section & rootXML);
    virtual ~RefineGeometry();
 
    RefineGeometry(const RefineGeometry & other);
    RefineGeometry & operator=(RefineGeometry other);
 
    void buildPorts();
 
    RefineType  getRefineType();
 
    void    setRefineType(RefineType type);
    void    setRefineType(int type);
    void    setRefineSteps(int steps);
    void    setSmoothingSteps(int steps);
    void    clear();
    void    execute();
 
    virtual void absorbSectionXML(const bitpit::Config::Section & slotXML, std::string name="");
    virtual void flushSectionXML(bitpit::Config::Section & slotXML, std::string name="");
 
    void     plotOptionalResults();
protected:
    void swap(RefineGeometry & x) noexcept;
 
    void ternaryRefine(std::unordered_map<long,long> * mapping = nullptr, mimmo::MimmoSharedPointer<MimmoObject> coarsepatch = nullptr, mimmo::MimmoSharedPointer<MimmoObject> refinepatch = nullptr);
    std::vector<long> ternaryRefineCell(const long & cellId, const std::vector<bitpit::Vertex> & vertices, const std::array<double,3> & center);
    void redgreenRefine(std::unordered_map<long,long> * mapping = nullptr, mimmo::MimmoSharedPointer<MimmoObject> coarsepatch = nullptr, mimmo::MimmoSharedPointer<MimmoObject> refinepatch = nullptr);
    std::vector<long> redRefineCell(const long & cellId, const std::vector<long> & newVertexIds);
    std::vector<long> greenRefineCell(const long & cellId, const long newVertexId, int iface);
    void smoothing(std::set<long> * constrainedVertices = nullptr);
    bool checkTriangulation();
};
 
REGISTER_PORT(M_GEOM, MC_SCALAR, MD_MIMMO_, __REFINEGEOMETRY_HPP__)
REGISTER_PORT(M_GEOM2, MC_SCALAR, MD_MIMMO_, __REFINEGEOMETRY_HPP__)
 
REGISTER(BaseManipulation, RefineGeometry, "mimmo.RefineGeometry")
 
};
 
#endif /* __REFINEGEOMETRY_HPP__ */