u3d_bspl3_surface4d.hpp

/* u3d_bspl3_surface4d.hpp */
/* $Id: u3d_bspl3_surface4d.hpp,v 1.2 2008/03/25 11:46:58 ibadm Exp $ */
#ifndef __u3d_bspl3_surface4d_hpp
#define __u3d_bspl3_surface4d_hpp

#include "u3d_math_abstract_surface.hpp"
#include "u3d_bspl3_primitives.hpp"
#include "u3d_prgrid3d.hpp"

class  U3dBspl3Surface4d 
{
public:
  U3dBspl3Surface4d(const U3dBspl3AbstractMatrix &bm, const U3dPoints4dAr &points, double g_step);
  
  U3dBspl3Surface4d(const U3dBspl3Surface4d &obj);
  
  virtual ~U3dBspl3Surface4d();
  
  
  virtual void approximate(const U3dBspl3AbstractMatrix &bm, const U3dPoints4dAr &points, double g_step);
  
  
  virtual bool getPoint(double p, double q, double r, U3dPoint4d &point);
  
  void getAreaSize(U3dPoint4d &min, U3dPoint4d &max);
  
  bool isInterpolated()
  {return m_is_ip_done;}
  
  U3dBspl3MatrixType getType();
  
protected:
  
  typedef std::vector<U3dSimpleBspl3Surf3d*> BsplLine;
  
  typedef std::vector<BsplLine> BsplMatrix;
  
  std::vector<BsplMatrix> m_bspls_x4;
  
  std::vector<U3dSimpleBspl3Curve*> m_bspls_x1, m_bspls_x2, m_bspls_x3;
  
  U3dPRGrid3d *m_grid; 
  
  U3dBspl3AbstractMatrix *m_bm;
  
  bool m_is_ip_done;
    
  void clear();
  
  bool get_expanded_gp(int i, int j, int k, U3dPoint4d &point);
};


#endif /* u3d_bspl3_surface4d.hpp */