gsCurve.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <gsUtils/gsPointGrid.h>
#include <gsUtils/gsMesh/gsMesh.h>
 
namespace gismo
{
 
template<class T>
short_t gsCurve<T>::degree() const
{
    return this->basis().degree(0);
}
 
template<class T>
void gsCurve<T>::toMesh(gsMesh<T> & msh, int npoints) const
{
    const gsMatrix<T> param        = this->parameterRange();
    const gsMatrix<T> paramSamples = gsPointGrid<T>(param(0,0), param(0,1), npoints);
    //const bool        zzero        = ( this->geoDim()==2 );
 
    gsMatrix<T> cp; // Curve point
 
    this->eval_into(paramSamples.col(0), cp);
 
    typename gsMesh<T>::VertexHandle v1,
    //v0 = msh.addVertex( cp(0,0), cp(1,0), zzero ? (T)(0) : cp(2,0) );
    v0 = msh.addVertex( cp );
 
    for ( int i = 1; i<npoints; ++i)
    {
        this->eval_into(paramSamples.col(i), cp);
        //v1 = msh.addVertex( cp(0,0), cp(1,0), zzero ? (T)(0) : cp(2,0) );
        v1 = msh.addVertex( cp );
        msh.addEdge(v0, v1);
        v0 = v1;
    }
}
 
} // namespace gismo