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