u3d_svrg_math_bspl3_surface.hpp

// u3d_svrg_math_bspl3_surface.hpp
// $Id: u3d_svrg_math_bspl3_surface.hpp,v 1.8 2007/06/26 13:21:32 guser1 Exp $
#ifndef U3D_SVRG_MATH_BSPL3_SURFACE_
#define U3D_SVRG_MATH_BSPL3_SURFACE_

#include "u3d_svrg_math_abstract_surface.hpp"
#include "u3d_bspl3_family.hpp"


class U3dSVRGMathBspl3Surface : public U3dSVRGMathAbstractSurface
{
  friend class U3dSVRGMathBspl3Surf4d;
  
public:
  U3dSVRGMathBspl3Surface(U3dPointsAr &points, double g_step);
  
  U3dSVRGMathBspl3Surface(U3dSVRGMathBspl3Surface *surf);
  
  virtual ~U3dSVRGMathBspl3Surface();
  
  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);

protected:
  
  typedef std::vector<U3dSimpleBspl3Surface*> BsplLine;
  
  std::vector<U3dSimpleBspl3Curve*> m_bspls_x, m_bspls_y;
  
  std::vector<BsplLine> m_bspls_z;
  
  U3dBspl3AbstractMatrix *m_bm;
  
  double m_pq_accur;
  
  virtual bool approximate();
  
  static bool find_element(double p, int n, int &nnp, double &pp);
  
  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;
  }
  
  static void update_index (int &i, int m) {
    if (i < 2)
      i = 0;
    else if (i >= m + 1)
      i = m - 1;
    else
      i--;
  }
  
};

#endif /*U3D_SVRG_MATH_BSPL3_SURFACE_*/