gsHBSplineBasis.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsHSplines/gsHTensorBasis.h>
#include <gsHSplines/gsHBSpline.h>
 
namespace gismo
{
template<short_t d, class T>
class gsHBSplineBasis : public gsHTensorBasis<d,T>
{
public:
    typedef gsHBSpline<d,T> GeometryType;
    
    typedef typename gsHTensorBasis<d,T>::CMatrix CMatrix;
 
    typedef typename gsHTensorBasis<d,T>::tensorBasis tensorBasis;
 
    typedef typename 
    util::conditional<d==1, gsConstantBasis<T>, gsHBSplineBasis<static_cast<short_t>(d-1),T>
                      >::type BoundaryBasisType;
 
    typedef memory::shared_ptr< gsHBSplineBasis > Ptr;
 
    typedef memory::unique_ptr< gsHBSplineBasis > uPtr;
 
public:
 
    gsHBSplineBasis() { }
 
    gsHBSplineBasis(gsBasis<T> const&  tbasis, bool manualLevels=false )
        : gsHTensorBasis<d,T>(tbasis, manualLevels) 
    {
        // initialize(); // is done in the base constructor
    }
    
    gsHBSplineBasis( gsTensorBSplineBasis<d,T> const&  tbasis,
                     std::vector<index_t> & boxes)
        : gsHTensorBasis<d,T>(tbasis, boxes) 
    {
        // initialize(); // is done in the base constructor
    }
    
    gsHBSplineBasis( gsTensorBSplineBasis<d,T> const&  tbasis,
                     gsMatrix<T> const & boxes)
        : gsHTensorBasis<d,T>(tbasis, boxes) 
    {
        // initialize(); // is done in the base constructor
    }
 
#ifdef __DOXYGEN__
    typename BoundaryBasisType::uPtr boundaryBasis(boxSide const & s);
#endif
    GISMO_UPTR_FUNCTION_DEF(BoundaryBasisType, boundaryBasis, boxSide const &)
    {
        return basisSlice(n1.direction(),n1.parameter());
    }
 
public:
    BoundaryBasisType * basisSlice(index_t dir_fixed,T par ) const;
 
    short_t domainDim() const { return d; }
    
    void eval_into(const gsMatrix<T> & u, gsMatrix<T>& result) const;
 
    void deriv_into(const gsMatrix<T> & u, gsMatrix<T>& result) const;
 
    void deriv2_into(const gsMatrix<T> & u, gsMatrix<T>& result) const;
    
    void evalSingle_into(index_t i, const gsMatrix<T> & u, gsMatrix<T>& result) const;
    
    void derivSingle_into(index_t i, const gsMatrix<T> & u, gsMatrix<T>& result) const;
    
    void deriv2Single_into(index_t i, const gsMatrix<T> & u, gsMatrix<T>& result) const;
 
    GISMO_CLONE_FUNCTION(gsHBSplineBasis)
 
    
    std::ostream &print(std::ostream &os) const;
    void transferbyLvl(std::vector<gsSparseMatrix<T> >& result);
 
    GISMO_MAKE_GEOMETRY_NEW
    
private:
    
    void initialize();
 
    gsSparseMatrix<T> coarsening(const std::vector<gsSortedVector<index_t> >& old, const std::vector<gsSortedVector<index_t> >& n, const gsSparseMatrix<T,RowMajor> & transfer) const;
    gsSparseMatrix<T> coarsening_direct( const std::vector<gsSortedVector<index_t> >& old, const std::vector<gsSortedVector<index_t> >& n, const std::vector<gsSparseMatrix<T,RowMajor> >& transfer) const;
 
    gsSparseMatrix<T> coarsening_direct2( const std::vector<gsSortedVector<index_t> >& old, const std::vector<gsSortedVector<index_t> >& n, const std::vector<gsSparseMatrix<T,RowMajor> >& transfer) const;
 
    using gsHTensorBasis<d,T>::m_bases;
    using gsHTensorBasis<d,T>::m_xmatrix;
    using gsHTensorBasis<d,T>::m_xmatrix_offset;
    using gsHTensorBasis<d,T>::m_deg;
    
}; // class gsHBSplineBasis
 
#ifdef GISMO_WITH_PYBIND11
 
  void pybind11_init_gsHBSplineBasis2(pybind11::module &m);
  void pybind11_init_gsHBSplineBasis3(pybind11::module &m);
  void pybind11_init_gsHBSplineBasis4(pybind11::module &m);
 
#endif // GISMO_WITH_PYBIND11
 
} // namespace gismo
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsHBSplineBasis.hpp)
#endif