u_bspl3_curve1d.hpp

/* u_bspl3_curve1d.hpp */
/* $Id$ */
#ifndef __u_bspl3_curve1d_hpp
#define __u_bspl3_curve1d_hpp

#include <mth/u_curve1d.hpp>
#include <mth/u3d_bspl3_matrices.hpp>
#include <mth/u3d_bspl3_primitives.hpp>

class UBspl3Curve1d
{
public:
  
  UBspl3Curve1d();
  
  UBspl3Curve1d(const U3dBspl3AbstractMatrix &bm);
  
  UBspl3Curve1d(const U3dBspl3AbstractMatrix &bm, 
        const std::vector<double> &ps);
  
  UBspl3Curve1d(const UBspl3Curve1d &obj);
  
  virtual ~UBspl3Curve1d();
  
  
  void approximate(const U3dBspl3AbstractMatrix &bm, 
           const std::vector<double> &ps);
  
  void approximate(const std::vector<double> &ps);
    
  double operator()(double p) const;
  
  double getFirstDerivative(double p) const;

  double getSecondDerivative(double p) const;

  
  bool ready() const
  {return m_is_ready;}
  
  
  U3dBspl3MatrixType getBMType() const;
  
  
  virtual void addNode(double x, int id = -1);
  
  virtual bool removeNode(int id);
  
  
  void info() const;
  
  
  static void update_index (int &i, int m, int bm) {
    if (i <= bm)
      i = 0;
    else if (i >= m + bm - 1)
      i = m - 1;
    else
      i -= bm;
  }

  static int get_bp_num(U3dBspl3MatrixType type); 
  
  static void find_element(double p, int n, int &nn, double &pp);
  
protected:
        
  std::vector<U3dSimpleBspl3Curve> m_bspls;
  
  U3dBspl3AbstractMatrix *m_bm;
  
  std::vector<double> m_ps;
  
  bool m_is_ready;
  
  void clear();
  
  void construct_segment(int id, int n, int bpnum, U3dSimpleBspl3Curve &curve);
};


#endif /* u_bspl3_curve1d.hpp */