gsNurbsBasis.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <gsNurbs/gsNurbs.h>
#include <gsNurbs/gsTensorNurbs.h>
#include <gsIO/gsXml.h>
#include <gsIO/gsXmlGenericUtils.hpp>
 
namespace gismo
{
 
template <class T>
gsNurbsBasis<T>::~gsNurbsBasis() { }
 
template <class T>
typename gsNurbsBasis<T>::gsGeoPtr
gsNurbsBasis<T>::makeGeometry( gsMatrix<T> coefs ) const
{ return gsGeoPtr(new GeometryType(*this, give(coefs))); }
 
template <class T>
typename gsNurbsBasis<T>::gsBasisPtr
gsNurbsBasis<T>::create(std::vector<KnotVectorType> cKV, gsMatrix<T> weights)
{
    const index_t dd = cKV.size();
    switch (dd)
    {
    case 1:
        return gsBasisPtr(new gsNurbsBasis<T>(give(cKV.front()), give(weights)));
        break;
    case 2:
        return gsBasisPtr(new gsTensorNurbsBasis<2,T>(give(cKV), give(weights)));
        break;
    case 3:
        return gsBasisPtr(new gsTensorNurbsBasis<3,T>(give(cKV), give(weights)));
        break;
    case 4:
        return gsBasisPtr(new gsTensorNurbsBasis<4,T>(give(cKV), give(weights)));
        break;
    }
    GISMO_ERROR("Dimension should be between 1 and 4.");
}
 
namespace internal {
template<class T>
class gsXml< gsNurbsBasis<T> >
{
private:
    gsXml() { }
public:
    GSXML_COMMON_FUNCTIONS(gsNurbsBasis<T>);
    static std::string tag () { return "Basis"; }
    static std::string type () { return "NurbsBasis"; }
 
    static gsNurbsBasis<T> * get (gsXmlNode * node)
    {
        return getRationalBasisFromXml<gsNurbsBasis<T> >(node);
    }
 
    static gsXmlNode * put (const gsNurbsBasis<T> & obj,
                            gsXmlTree & data )
    {
        return putRationalBasisToXml(obj,data);
    }
};
 
}
 
} // namespace gismo