u3d_bspl3_curve3d.hpp

// u3d_bspl3_curve3d.hpp
// $Id: u3d_bspl3_curve3d.hpp,v 1.4 2009/06/30 13:06:22 urij Exp $
#ifndef U3D_BSPL3_CURVE3D_
#define U3D_BSPL3_CURVE3D_

#include <mth/u3d_bspl3_primitives.hpp>
#include <mth/u3d_point.hpp>

class U3dBspl3Curve3d
{
public:
  U3dBspl3Curve3d(const U3dBspl3AbstractMatrix &bm);
  
  U3dBspl3Curve3d(const U3dBspl3AbstractMatrix &bm, const U3dPointsAr &ps);
  
  U3dBspl3Curve3d(const U3dBspl3Curve3d &obj);
  
  ~U3dBspl3Curve3d();
  
  
  void approximate(const U3dBspl3AbstractMatrix &bm, const U3dPointsAr &ps);
  
  void approximate(const U3dPointsAr &ps);
  
  
  bool getPoint(double p, U3dPoint &point) const;
  
  bool getXYZp(double p, U3dPoint &point) const;
  
  bool getXYZpp(double p, U3dPoint &point) const;
  
  bool getCurvature(double p, double &c);
  
  int getControlPointsNum() const;
  
  bool getControlPoint(int ind, U3dPoint &point) const;
  
  U3dBspl3MatrixType getBMType() const;
  
  
  void add(const U3dPoint &point, int id = -1);
  
  bool remove(int id);
  
  
  bool isReady() const
  {return m_is_ready;}
  
  
  bool getLineSegments(double ds, std::vector<double> &segs) const;
  
  bool getLineDeviations(double p1, double p2, double p, U3dPoint &dev) const;
  
  bool getMaxLineDeviation(double p1, double p2, double &md) const;
  
  
  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:
  
  struct Bspl3dCurve 
  {
    Bspl3dCurve() : x(0x0), y(0x0), z(0x0) {};
    
    void clear() {
      delete x; delete y; delete z;
    }
    
    U3dSimpleBspl3Curve *x, *y, *z;
  }; 
      
  std::vector<Bspl3dCurve> m_bspls;
  
  U3dBspl3AbstractMatrix *m_bm;
  
  U3dPointsAr m_ps;
  
  bool m_is_ready;
  
  void clear();
  
  void construct_segment(int id, int n, int bpnum, Bspl3dCurve &curve);
  
  bool get_line_distance(double pp, int nn, U3dPoint kP, U3dPoint bP, double p1, U3dPoint &dist) const;
  
  bool construct_isect_line(double p1, double p2, U3dPoint &kP, U3dPoint &bP) const;
  
  void devs_update(double p, int n,  U3dPoint &kP, U3dPoint &bP, double p1,
           U3dPoint &dev, double &smx, double &smy, double &smz) const {
    get_line_distance(p, n, kP, bP, p1, dev);
    if (dev.getX() > smx)
      smx = dev.getX();
    if (dev.getY() > smy)
      smy = dev.getY();
    if (dev.getZ() > smz)
      smz = dev.getZ();
  }

};

#endif /* U3D_BSPL3_CURVE3D_ */