u3d_bspl3_curve3d.hpp

// u3d_bspl3_curve3d.hpp
// $Id: u3d_bspl3_curve3d.hpp,v 1.11 2007/09/07 12:17:03 guser1 Exp $
#ifndef U3D_BSPL3_CURVE3D_
#define U3D_BSPL3_CURVE3D_

#include <u3d_bspl3_family.hpp>

class U3dBspl3Curve3d
{
public:
  U3dBspl3Curve3d(const U3dBspl3AbstractMatrix &bm, const U3dPointsAr &ps);
  
  ~U3dBspl3Curve3d();
  
  
  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;
  
  
  bool add(U3dPoint &point, int ind = -1);
  
  bool remove(int ind);
  
  bool move(int ind);
  
  
  bool isReady() const
  {return m_is_ready;}
  
  
  bool getLineSegments(double ds, std::vector<double> &segs);
  
  bool getLineDeviations(double p1, double p2, double p, U3dPoint &dev);
  
  bool getMaxLineDeviation(double p1, double p2, double &md);
  
  
  void info();
  
protected:
  
  struct Bspl3dCurve 
  {
    U3dSimpleBspl3Curve *x, *y, *z;
  }; 
      
  std::vector<Bspl3dCurve> m_bspls;
  
  U3dBspl3AbstractMatrix *m_bm;
  
  U3dPointsAr m_ps;
  
  bool m_is_ready;
  
  void reconstruct();
  
  void clear();
  
  bool check_add_bps(const U3dPointsAr &ps);
  
  bool find_element(double p, int &nn, double &pp) const;
  
  bool get_line_distance(double pp, int nn, U3dPoint kP, U3dPoint bP, double p1, U3dPoint &dist);
  
  bool construct_isect_line(double p1, double p2, U3dPoint &kP, U3dPoint &bP);
  
  void devs_update(double p, int n,  U3dPoint &kP, U3dPoint &bP, double p1,
           U3dPoint &dev, double &smx, double &smy, double &smz) {
    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();
  }
  
  void next_segment (double p1, double p2, std::vector<double> &segs);
};

#endif /* U3D_BSPL3_CURVE3D_ */