u3d_bspl3_primitives.hpp

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

#include "u3d_bspl3_matrices.hpp"
#include <vector>

//------------------------------------------------------------------------------------------
class U3dSimpleBspl3Curve
{
public:
  
  U3dSimpleBspl3Curve();
  
  U3dSimpleBspl3Curve(const U3dBspl3AbstractMatrix *bm, const double (&cp)[4]);
  
  virtual ~U3dSimpleBspl3Curve(){};
  

  void construct(const U3dBspl3AbstractMatrix *bm, const double (&cp)[4]);
  

  double get(double p) const;
  
  double getSp(double p) const;
  
  double getSpp(double p) const;
  
  void getCoeffs(double (&c)[4]) const {
    c[0] = m_rv[0];
    c[1] = m_rv[1];
    c[2] = m_rv[2];
    c[3] = m_rv[3];
  }
    
  double getLineDeviation(double k, double b, double p);
  
protected:
  
  const U3dBspl3AbstractMatrix *m_bm;
  
  double m_rv[4];
};

typedef std::vector<U3dSimpleBspl3Curve*> U3dSimpleBspl3CurveAr;

//------------------------------------------------------------------------------------------
class U3dSimpleBspl3Surface
{
public:
  
  U3dSimpleBspl3Surface(const U3dBspl3AbstractMatrix *bm, const double (&cp)[4][4]);
  
  virtual ~U3dSimpleBspl3Surface(){};
  
  double get(double p, double q) const;
  
  double getSp(double p, double q) const;
  
  double getSq(double p, double q) const;
  
  double getSpp(double p, double q) const;
  
  double getSqq(double p, double q) const;
  
  double getSpq(double p, double q) const;
  
protected:
  
  const U3dBspl3AbstractMatrix *m_bm;
  
  double m_rm[4][4];
  
  double calc(double *pp, double *qq) const;
  
private:
  
  void cp2rm(const U3dBspl3AbstractMatrix *bm, const double (&cp)[4][4]);
};


//------------------------------------------------------------------------------------------
class U3dSimpleBspl3Surf3d
{
public:
  
  U3dSimpleBspl3Surf3d(const U3dBspl3AbstractMatrix *bm, const double (&cp)[4][4][4]);
  
  virtual ~U3dSimpleBspl3Surf3d(){};
  
  double get(double p, double q, double r) const;
  
protected:
  
  double m_rt[4][4][4];
  
  const U3dBspl3AbstractMatrix *m_bm;
  
  void make_tensor(const U3dBspl3AbstractMatrix *bm, const double (&cp)[4][4][4]);

  double calc(double *pp, double *qq, double *rr) const;
};


#endif /* u3d_bspl3_primitives.hpp */