gsCPPInterface.h

Go to the documentation of this file.

#pragma once

#include <gsCore/gsAffineFunction.h>
#include <gsCore/gsBoundary.h>
#include <gsIO/gsOptionList.h>

namespace gismo {

template <class T>
class gsCPPInterface : public gsFunction<T>
{
public:

    typedef memory::shared_ptr< gsCPPInterface > Ptr;

    typedef memory::unique_ptr< gsCPPInterface > uPtr;

    gsCPPInterface(const gsMultiPatch<T>   & mp,
                     const gsMultiBasis<T>   & mb,
                     const boundaryInterface & bi,
                     const gsOptionList      & opt = defaultOptions() );

    static uPtr make (const gsMultiPatch<T>   & mp,
                     const gsMultiBasis<T>   & mb,
                     const boundaryInterface & bi,
                     const gsOptionList      & opt = defaultOptions() )
    { return uPtr(new gsCPPInterface(mp,mb,bi,opt)); }
    
private:
    const gsGeometry<T>* m_slaveGeom, *m_masterGeom; 
    typename gsGeometry<T>::Ptr m_masterBdr;  

    const gsBasis<T> * m_slaveBasis ;        
    const gsBasis<T> * m_masterBasis;        

    boundaryInterface m_boundaryInterface;   

    T m_Tolerance;                            
    
    std::vector<index_t> m_freeDirs;
    index_t m_fixedParam;
    index_t m_fixedDir;

    std::vector< std::vector<T> > m_breakpoints;  

public:
    static gsOptionList defaultOptions()
    {
        gsOptionList defaultOpts;
        defaultOpts.addInt ( "Tied",      "Is it a tied contact?", 0);
        defaultOpts.addReal( "Tolerance", "Tolerance to be used for closest point search.", 1e-5);
        return defaultOpts;
    }


    void updateBdr();


    virtual void eval_into(const gsMatrix<T> & u, gsMatrix<T> & result) const;

    virtual short_t domainDim() const { return m_slaveGeom->domainDim(); }

    typename gsDomainIterator<T>::uPtr makeDomainIterator() const;


    const std::vector< std::vector<T> > & breakPoints() const { return m_breakpoints; }

    virtual std::ostream & print(std::ostream& os) const;

private:
    void constructInterfaceBox();

    void constructBreaks();
}; // End gsCPPInterface


template <class T>
inline std::ostream & operator<<(std::ostream & os, const gsCPPInterface<T> & cppIf)
{ return cppIf.print(os); }

} // End namespace gismo

#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsCPPInterface.hpp)
#endif