gsTHBSpline.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsCore/gsLinearAlgebra.h>
#include <gsHSplines/gsTHBSplineBasis.h>
#include <gsCore/gsForwardDeclarations.h>
 
namespace gismo
{
 
template<short_t d, class T>
class gsTHBSpline : public gsGeoTraits<d,T>::GeometryBase
{
public:
    typedef gsTHBSplineBasis<d,T> Basis;
 
    typedef typename Basis::tensorBasis tensorBasis;
 
    typedef typename gsGeoTraits<d,T>::GeometryBase Base;
 
    typedef memory::shared_ptr< gsTHBSpline > Ptr;
 
    typedef memory::unique_ptr< gsTHBSpline > uPtr;
 
    typedef typename
    util::conditional<d==1, gsConstantFunction<T>, gsTHBSpline<static_cast<short_t>(d-1),T>
                      >::type BoundaryGeometryType;
 
    typedef typename gsTHBSplineBasis<d,T>::BoundaryBasisType BoundaryBasisType;
 
public:
 
    gsTHBSpline() { }
 
    gsTHBSpline( const Basis & basis, const gsMatrix<T> & coefs ) :
    Base( basis, coefs ) 
    { }
 
    gsTHBSpline( const Basis& basis, gsMatrix<T> & coefs ) :
        Base( basis, coefs )
    { }
 
    explicit gsTHBSpline( const gsTensorBSpline<d,T> & tbsp )
    {
        this->m_basis = new Basis( tbsp.basis() );
        this->m_coefs = tbsp.coefs();
    }
 
    GISMO_CLONE_FUNCTION(gsTHBSpline)
 
    GISMO_BASIS_ACCESSORS
 
public:
 
 
// ************************************************
//  Other member functions
// ************************************************
public:
    void getBsplinePatches(gsMatrix<index_t>& b1, gsMatrix<index_t>& b2, gsVector<index_t>& level) const;
    //\param[out] bpatches list of B-spline patches associated with the boxes
    //void getBsplinePatches(gsMatrix<index_t>& b1, gsMatrix<index_t>& b2, gsVector<index_t>& level, std::vector< gsTensorBSpline<2> > & bpatches) const;
 
//    void convertToBSpline( gsTensorBSpline<d,T,gsCompactKnotVector<T> >& result );
 
    void convertToBSpline( gsTensorBSpline<d,T>& result );
 
    void increaseMultiplicity(index_t lvl, int dir, T knotValue, int mult = 1);
 
 
private:
 
    void getBsplinePatchGlobal(gsVector<index_t> b1, gsVector<index_t> b2, unsigned l, gsTensorBSpline<2> geo) const;
    void globalRefinement(int level)const;
 
    void initialize_cmatrix(int col, int c_level) const;
 
    void return_cp_1D(const gsMatrix<T> & mat, int direction, gsMatrix<T>& cp)const;
 
public:
 
    void slice(index_t dir_fixed,T par,BoundaryGeometryType & result) const;
 
}; // class gsTHBSpline
 
#ifdef GISMO_WITH_PYBIND11
 
  void pybind11_init_gsTHBSpline2(pybind11::module &m);
  void pybind11_init_gsTHBSpline3(pybind11::module &m);
  void pybind11_init_gsTHBSpline4(pybind11::module &m);
 
#endif // GISMO_WITH_PYBIND11
 
} // namespace gismo
 
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsTHBSpline.hpp)
#endif