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>