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