u3d_svrg_math_bs_surface.hpp

// u3d_svrg_math_bs_surface.hpp
// $Id: u3d_svrg_math_bs_surface.hpp,v 1.3 2007/06/26 13:21:32 guser1 Exp $
#ifndef U3D_SVRG_MATH_BS_SURFACE_
#define U3D_SVRG_MATH_BS_SURFACE_

#include "u3d_svrg_math_abstract_surface.hpp"

class U3dSVRGMathBSSurface : public U3dSVRGMathAbstractSurface
{
public:
  U3dSVRGMathBSSurface(U3dPointsAr &points, double g_step, int ord_p = 4, int ord_q = 4);
  
  U3dSVRGMathBSSurface(U3dSVRGMathBSSurface *surf);
  
  virtual ~U3dSVRGMathBSSurface();
    
  virtual bool getPoint(double p, double q, U3dPoint &point);
  
  //virtual bool getCurvature(double p, double q, double &c);
  
  //virtual bool getMinCurvature(double p1, double q1, double p2, double q2, 
  //               double &c);
  
  int getOrderP()
  {return m_ord_p;}
  
  int getOrderQ()
  {return m_ord_q;}
  
  void setOrderP(int ord)
  {m_ord_p = ord; approximate();}
  
  void setOrderQ(int ord)
  {m_ord_q = ord; approximate();}
  
  void movePoint(int ind, double x, double y, double z);

  void addPoint(U3dPoint *point);
  
  
protected:
  int m_ord_p, m_ord_q;
  
  double *m_knots_p;
  
  double *m_knots_q;
  
  bool approximate();
  
  void make_knots(int dir, double *&knots);
  
  double getBSpl(double t, double *x, int i, int k);    
  
private:
  double *m_Np, *m_Nq; 
};

#endif /*U3D_SVRG_MATH_BS_SURFACE_*/