gsCurve.h

Go to the documentation of this file.

#pragma once

#include <gsCore/gsGeometry.h>

namespace gismo
{

template<class T>
class gsCurve : public gsGeometry<T>
{
    
public:
    typedef memory::shared_ptr< gsCurve > Ptr;

    typedef memory::unique_ptr< gsCurve > uPtr;
    
    typedef T Scalar_t;
public:


    gsCurve() : gsGeometry<T>() { }

    gsCurve(const gsBasis<T> & basis, gsMatrix<T> coefs ) :
    gsGeometry<T>(basis, give(coefs))
    {
        GISMO_ASSERT( this->m_coefs.size() >= 1,
        "Coefficient matrix cannot be empty.\n");
    }


//    GISMO_UPTR_FUNCTION_PURE(gsCurve, clone)
private: virtual gsCurve * clone_impl() const = 0;
public: inline uPtr clone() const { return uPtr(clone_impl()); }

    short_t domainDim() const { return 1; }
    
    short_t degree() const;

    void toMesh(gsMesh<T> & msh, int npoints = 100) const;

    void reverse()
    {
        this->m_coefs = this->m_coefs.colwise().reverse().eval();
        this->basis().reverse();
    }

    virtual bool isOn(gsMatrix<T> const &, T = 1e-3) const
    { GISMO_NO_IMPLEMENTATION }

}; // class gsCurve

} // namespace gismo


#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsCurve.hpp)
#endif