OBBox.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 __OBBox_HPP__
#define __OBBox_HPP__
 
#include "BaseManipulation.hpp"
 
namespace mimmo{
 
 
enum class OBBStrategy{
    OBB = 0 ,
    AABB = 1 ,
    MINVOL = 2 
};
 
 
 
class OBBox: public BaseManipulation {
 
protected:
    darray3E    m_origin;       
    darray3E    m_span;         
    dmatrix33E    m_axes;       
    std::unordered_map<MimmoSharedPointer<MimmoObject>, int> m_listgeo; 
    OBBStrategy m_strategy; 
    bool m_writeInfo; 
public:
    OBBox();
    OBBox(const bitpit::Config::Section & rootXML);
    virtual ~OBBox();
 
    //copy operators/constructors
    OBBox(const OBBox & other);
 
    void buildPorts();
 
    //clean structure;
    void         clearOBBox();
 
    //internal methods
    std::vector<MimmoSharedPointer<MimmoObject> >   getGeometries();
    darray3E                         getOrigin();
    darray3E                         getSpan();
    dmatrix33E                       getAxes();
    BITPIT_DEPRECATED(
    bool                             isForcedAABB());
    OBBStrategy                      getOBBStrategy();
 
    void        setGeometry(MimmoSharedPointer<MimmoObject> geo);
    void        setGeometries(std::vector<MimmoSharedPointer<MimmoObject> > listgeo);
    BITPIT_DEPRECATED(
    void        setForceAABB(bool flag));
    void        setOBBStrategy(OBBStrategy strategy);
    void        setOBBStrategyInt(int strategyflag);
    void        setWriteInfo(bool flag);
 
    //plotting wrappers
    void        plot(std::string directory, std::string filename, int counter, bool binary);
 
    //building method
    void execute();
 
    virtual void absorbSectionXML(const bitpit::Config::Section & slotXML, std::string name = "");
    virtual void flushSectionXML(bitpit::Config::Section & slotXML, std::string name= "");
 
 
protected:
    virtual void plotOptionalResults();
    void swap(OBBox & x) noexcept;
    dmatrix33E transpose(const dmatrix33E & mat);
    dmatrix33E inverse (const dmatrix33E & mat);
    void computeOBB(std::vector<MimmoSharedPointer<MimmoObject>> & vector_listgeo, darray3E& origin, darray3E & span, dmatrix33E & axes);
    void computeAABB(std::vector<MimmoSharedPointer<MimmoObject>> & vector_listgeo, darray3E& origin, darray3E & span, dmatrix33E & axes);
 
private:
    dmatrix33E      evaluatePointsCovarianceMatrix(std::vector<MimmoSharedPointer<MimmoObject> > list);
    dmatrix33E      evaluateElementsCovarianceMatrix(std::vector<MimmoSharedPointer<MimmoObject> > list);
    dmatrix33E      eigenVectors( dmatrix33E &, darray3E & eigenValues);
    void            adjustBasis( dmatrix33E &, darray3E & eigenValues);
    std::array<long,3> get3RepPoints(long cellID, mimmo::MimmoSharedPointer<MimmoObject> geo);
};
 
REGISTER_PORT(M_GEOM, MC_SCALAR, MD_MIMMO_,__OBBox_HPP__)
REGISTER_PORT(M_VECGEOM, MC_VECTOR, MD_MIMMO_,__OBBox_HPP__)
REGISTER_PORT(M_POINT, MC_ARRAY3, MD_FLOAT,__OBBox_HPP__)
REGISTER_PORT(M_AXES, MC_ARR3ARR3, MD_FLOAT,__OBBox_HPP__)
REGISTER_PORT(M_SPAN, MC_ARRAY3, MD_FLOAT,__OBBox_HPP__)
 
 
REGISTER(BaseManipulation, OBBox, "mimmo.OBBox")
 
};
 
#endif /* __LATTICE_HPP__ */