gsMPBESSpline.h

Go to the documentation of this file.

#pragma once

#include <gsCore/gsMultiPatch.h>

#include <gsMSplines/gsMappedSpline.h>
#include <gsMSplines/gsMappedSingleBasis.h>
#include <gsUnstructuredSplines/src/gsMPBESBasis.h>

namespace gismo
{
template<short_t d,class T>
class gsMPBESSpline : public gsMappedSpline<d,T>
{

private:
    typedef gsMappedSingleBasis<d,T> Basis;
    typedef gsMappedSpline<d,T> Base;
public:
    gsMPBESSpline() : Base() { }

    gsMPBESSpline( const gsMPBESBasis<d,T> & basis, const gsMatrix<T> & coefs );

    gsMPBESSpline( gsMultiPatch<T> const & mp,index_t incrSmoothness = -1,index_t minEVDistance = -1 );

    gsMPBESSpline( gsMultiPatch<T> const & mp,std::vector<patchCorner> C0List,index_t incrSmoothness = -1,index_t minEVDistance = -1 );

    ~gsMPBESSpline() {   } //destructor

    using Base::m_mbases;
    using Base::m_global;

    void setCornerC0(patchCorner const & pc);

    void smoothCornerEdge(const patchCorner& pc,const patchSide& ps,bool updateBasis=true);

    void smoothEverything();

    virtual void uniformRefine(index_t numKnots=1, index_t mul=1);

    void refine(const index_t patch, const gsMatrix<T> &boxes);

    void refineElements(const index_t patch, std::vector<index_t> const & boxes);

    void refineElements(std::vector<std::vector<index_t> > const & boxes);

    void uniformRefineAndSmooth(index_t numKnots=1);

    void refineAndSmooth(const index_t patch, const gsMatrix<T> &boxes);

    void refineElementsAndSmooth(const index_t patch, std::vector<index_t> const & boxes);

}; // class gsMPBESSpline

} // namespace gismo