BendGeometry.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 __BENDGEOMETRY_HPP__
#define __BENDGEOMETRY_HPP__
 
#include "BaseManipulation.hpp"
 
namespace mimmo{
 
class BendGeometry: public BaseManipulation{
private:
    darray3E            m_origin;       
    dmatrix33E          m_system;       
    bool                m_local;        
    umatrix33E          m_degree;        
    dmat33Evec          m_coeffs;        
    dmpvector1D           m_filter;       
    dmpvecarr3E           m_displ;        
public:
    BendGeometry();
    BendGeometry(const bitpit::Config::Section & rootXML);
    ~BendGeometry();
 
    BendGeometry(const BendGeometry & other);
    BendGeometry & operator =(BendGeometry other);
    void buildPorts();
 
    darray3E    getOrigin();
    dmatrix33E  getRefSystem();
    umatrix33E     getDegree();
    dmat33Evec*    getCoeffs();
    dmpvecarr3E*   getDisplacements();
 
    void    setFilter(dmpvector1D *filter);
    void    setOrigin(darray3E origin);
    void    setRefSystem(dmatrix33E axes);
    void    setDegree(umatrix33E degree);
    void    setDegree(int i, int j, uint32_t degree);
    void    setCoeffs(dmat33Evec *coeffs);
    void    setCoeffs(int i, int j, dvector1D coeffs);
 
    void     execute();
    void     apply();
 
    //XML utilities from reading writing settings to file
    virtual void absorbSectionXML(const bitpit::Config::Section & slotXML, std::string name="");
    virtual void flushSectionXML(bitpit::Config::Section & slotXML, std::string name="");
 
protected:
    void swap(BendGeometry & x) noexcept;
    void    checkFilter();
    static darray3E   matmul(const darray3E & vec, const dmatrix33E & mat);
    static darray3E   matmul(const dmatrix33E & mat,const darray3E & vec);
 
private:
    darray3E    toLocalCoord(darray3E point);
    darray3E    toGlobalCoord(darray3E point);
 
};
 
REGISTER_PORT(M_FILTER, MC_SCALAR, MD_MPVECFLOAT_,__BENDGEOMETRY_HPP__)
REGISTER_PORT(M_POINT, MC_ARRAY3, MD_FLOAT,__BENDGEOMETRY_HPP__)
REGISTER_PORT(M_AXES, MC_ARR3ARR3, MD_FLOAT,__BENDGEOMETRY_HPP__)
REGISTER_PORT(M_BMATRIX, MC_ARR3ARR3, MD_INT,__BENDGEOMETRY_HPP__)
REGISTER_PORT(M_BCOEFFS, MC_SCALAR, MD_MATRIXCOEFF_,__BENDGEOMETRY_HPP__)
REGISTER_PORT(M_GEOM, MC_SCALAR, MD_MIMMO_,__BENDGEOMETRY_HPP__)
REGISTER_PORT(M_GDISPLS, MC_SCALAR, MD_MPVECARR3FLOAT_,__BENDGEOMETRY_HPP__)
 
 
REGISTER(BaseManipulation, BendGeometry, "mimmo.BendGeometry")
 
}
 
#endif /* __BENDGEOMETRY_HPP__ */