/home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/core/spline/spline_nubs.h

Go to the documentation of this file.

#ifndef OPENTISSUE_CORE_SPLINE_SPLINE_NUBS_H
#define OPENTISSUE_CORE_SPLINE_SPLINE_NUBS_H
//
// OpenTissue Template Library
// - A generic toolbox for physics-based modeling and simulation.
// Copyright (C) 2008 Department of Computer Science, University of Copenhagen.
//
// OTTL is licensed under zlib: http://opensource.org/licenses/zlib-license.php
//
#include <OpenTissue/configuration.h>

#include <stdexcept>

namespace OpenTissue
{
  namespace spline
  {
    template<typename U, typename P>
    class NUBSpline
    {
    public:

      typedef U              knot_container;
      typedef P              point_container;

    protected:

      U      m_U;            
      P      m_P;            
      size_t m_k;            

    public:

      NUBSpline()
        : m_k(0u)
      {}

      NUBSpline(
          size_t const order
        , U const & knots
        , P const & controls
        )
      {

        if(order<= 0u)
          throw std::invalid_argument("order must be greater than zero");

        if(knots.size() <= 0u)
          throw std::invalid_argument("knots was empty");

        if(controls.size() <= 0u)
          throw std::invalid_argument("controls was empty");

        if( (knots.size() - controls.size())  != order)
          throw std::invalid_argument("Incorrect specification, n = m - k, was not fulfilled! ");

        this->m_U     = knots;
        this->m_P     = controls;        
        this->m_k     = order;
      }

      virtual ~NUBSpline() {      }

    public:

      U       & get_knot_container()       {   return this->m_U;   }
      U const & get_knot_container() const {   return this->m_U;   }

      P       & get_control_container()       {        return this->m_P;      }
      P const & get_control_container() const {        return this->m_P;      }

      size_t const get_dimension()   const 
      {
        if(this->m_P.size() == 0u)
          return 0u;
        return this->m_P[0].size(); 
      }

      size_t const & get_order() const { return this->m_k;   }

    };

  } // namespace spline
} // namespace OpenTissue

//OPENTISSUE_CORE_SPLINE_SPLINE_NUBS_H
#endif