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