u3d_bspl3_surface3d.hpp

/* u3d_bspl3_surface3d.hpp */
/* $Id: u3d_bspl3_surface3d.hpp,v 1.6 2008/06/04 12:22:37 sad Exp $ */
#ifndef __u3d_bspl3_surface3d_hpp
#define __u3d_bspl3_surface3d_hpp

#include "u3d_math_abstract_surface.hpp"
#include <mth/u3d_bspl3_primitives.hpp>
#include "u3d_abstract_prgrid2d.hpp"

#include "u3d_prgrid2d.hpp"

class U3dBspl3Surface3d : public U3dMathAbstractSurface
{
  
  friend class U3dBspl3Surface4d;
  
public:
  
  U3dBspl3Surface3d(const U3dBspl3AbstractMatrix &bm, const U3dAbstractPRGrid2d &grd);
  
  U3dBspl3Surface3d(const U3dBspl3Surface3d &surf);
  
  virtual ~U3dBspl3Surface3d();
  
  
  virtual void approximate(const U3dBspl3AbstractMatrix &bm, const U3dAbstractPRGrid2d &grd);
    
  
  virtual bool getPoint(double p, double q, U3dPoint &point);
  
  virtual bool getXYZp(double p, double q, U3dPoint &drp);
  
  virtual bool getXYZq(double p, double q, U3dPoint &drq);
  
  virtual bool getXYZpp(double p, double q, U3dPoint &ddrpp);
  
  virtual bool getXYZqq(double p, double q, U3dPoint &ddrqq);
  
  virtual bool getXYZpq(double p, double q, U3dPoint &ddrpq);
  
  const U3dAbstractPRGrid2d* grtGrid() {return m_grid;}
  
  U3dBspl3MatrixType getType();
  
  virtual int getInitialPointsNum();
  
  virtual U3dPoint* getInitialPoint(unsigned int ind);
  
  virtual void getInitialPoints(U3dPointsAr &ps);

  virtual const U3dPointsAr& getInitialPoints() const;

  
protected:
  
  typedef std::vector<U3dSimpleBspl3Surface*> BsplLine;
  
  U3dAbstractPRGrid2d *m_grid;
  
  std::vector<U3dSimpleBspl3Curve*> m_bspls_x, m_bspls_y;
  
  std::vector<BsplLine> m_bspls_z;
  
  U3dBspl3AbstractMatrix *m_bm;
  
  double m_pq_accur;
  
  bool get_expanded_gp(int i, int j, U3dPoint &point);
  
  static void correct_vars(double &p, double acc) {
    if (p > 1. - acc)
      p = 1. - acc;
    if (p < acc)
      p = acc;
  }
    
  void clear();
  
};

#endif /* u3d_bspl3_surface3d.hpp */