gsMonomialBasis.h

 
#pragma once
 
#include <iostream>
#include <gsCore/gsBasis.h>
 
 
namespace gismo
{
 
template<class T>
class gsMonomialBasis : public gsBasis<T>
{
public:
    typedef memory::shared_ptr< gsMonomialBasis > Ptr;
 
    typedef memory::unique_ptr< gsMonomialBasis > uPtr;
 
    static const int Dim = 1;
 
    gsMonomialBasis()  : m_p(0)
    { }
 
    explicit gsMonomialBasis(int degree);
    
    // Pure virtual member functions required by the base class
    virtual short_t domainDim() const;
 
    GISMO_CLONE_FUNCTION(gsMonomialBasis)
 
    virtual memory::unique_ptr<gsGeometry<T> > makeGeometry( gsMatrix<T> coefs ) const;
    virtual std::ostream &print(std::ostream &os) const;
    virtual index_t size() const;
 
 
    // Virtual member functions of the base class which need to be redefined
 
    virtual short_t degree() const { return m_p; }
 
    virtual short_t setDegree() const { return m_p; } // TODO: really?
 
    inline short_t deg() const { return m_p; }
 
    void active_into(const gsMatrix<T> & u, gsMatrix<index_t>& result) const
    {
        result.resize(m_p+1,u.cols());
        result.colwise() = 
            gsVector<index_t>::LinSpaced(m_p+1, 0, m_p); //.template cast<unsigned>();
    }
    
    virtual void eval_into(const gsMatrix<T> & u, gsMatrix<T>& result) const;
 
    virtual void evalSingle_into(index_t i, const gsMatrix<T> & u, gsMatrix<T>& result) const;
 
    virtual void deriv_into(const gsMatrix<T> & u, gsMatrix<T>& result ) const;
 
    virtual void derivSingle_into(index_t i, const gsMatrix<T> & u, gsMatrix<T>& result ) const;
 
    virtual void deriv2_into(const gsMatrix<T> & u, gsMatrix<T>& result ) const;
 
    virtual void deriv2Single_into(index_t i, const gsMatrix<T> & u, gsMatrix<T>& result ) const;
 
    virtual void evalAllDers_into(const gsMatrix<T> & u, int n, gsMatrix<T>& result) const;
    virtual void evalAllDersSingle_into(index_t i, const gsMatrix<T> & u, int n, gsMatrix<T>& result) const;
 
    virtual void evalDerSingle_into(index_t i, const gsMatrix<T> & u, int n, gsMatrix<T>& result) const;
 
    virtual gsMatrix<T> support() const
    {
        gsMatrix<real_t> support(1,2);
        support<<0, 1;
        return support;
    }
 
    virtual gsMatrix<T> support(const index_t & ) const
    {
        gsMatrix<real_t> support(1,2);
        support<<0, 1;
        return support;
    }
 
    short_t degree(short_t i) const { return m_p; }
 
    virtual void evalFunc_into(const gsMatrix<T> & u,
                               const gsMatrix<T> & coefs,
                               gsMatrix<T>& result) const;
 
 
 
// Data members
private:
    short_t m_p;
 
};  // class gsMonomialBasis
 
}   // namespce gismo
 
 
#include <gsPolynomial/gsMonomialBasis.hpp>