parspline.hpp

// parspline.hpp
// $Id: parspline.hpp,v 1.6 2008/08/18 08:19:18 sasha Exp $
//-------------------------------------------------------------------------
#ifndef ParSplineH
#define ParSplineH
//=========================================================================
#include "parspline_service.hpp"
//-------------------------------------------------------------------------

const int S_DEBUG_MODE = 1;

class ParSpline
{
 public:

  /* Constructor for empty object */
  ParSpline();

  ParSpline(double *xx, double *zz, double *ddx, double *ddz, double *ccurv,bool *flag, int nn);
  
  virtual ~ParSpline();
  
  void Interpolate(double *xx, double *zz, int nn, bool *flag, double rc, double sp,
           double *xx1 = NULL, double *zz1 = NULL, int nn1 = 0, bool noInter = false);
  
  void       GetXZ(double *xx, double *zz);

  void       GetSPXZ(double_ar *sx, double_ar *sz);

  bool       GetKnot(int kn, double &x, double &z);

  void       GetDXDZ(double *ddx, double *ddz);

  void       GetCurv(double *ccurv);

  void       GetSPFlags(bool *flag);

  int        GetN();

  bool       GetPFlag(int index);

  bool       GetPSpX(double z, double &x);

  bool       GetPSpXs(double z, double *&xs, int &nx);

  bool       GetPSpXs(double z, double *&xs, int &nx, double_ar *xx1, double_ar *zz1);

  bool       GetPSpZ(double x, double &z);

  bool       GetPSpZs(double x, double *&zs, int &nz);

  bool       GetPSpZs(double x, double *&zs, int &nz, double_ar *xx1, double_ar *zz1);

  bool       PSpInt(double x1, double z1, double x2, double z2, double &ix,
            double &iz, int *nn=NULL);

  bool       PSpInt1(double x1, double z1, double x2, double z2, double &ix,
             double &iz, int *nn=NULL);

  bool       PSpInts(double x1, double z1, double x2, double z2, double *ixs,
             double *izs, int &in, int *inns=NULL);

  bool       PSpI(double x1, double z1, double x2, double z2, int nn,
             double &ix, double &iz, double &angle, double &curvity);

  bool       PSpI1(double x1, double z1, double x2, double z2, int nn,
              double &ix, double &iz, double &angle);

  bool       PSpIs(double x1, double z1, double x2, double z2, int nn,
              double *ixs, double *izs, int &in, double &angle);

  bool       PSpCutMiddle(double x1, double z1, double x2, double z2,
              double &ix, double &iz);
  bool       isPSpClean();

  void       SetXZMinMax(double xmn, double xmx, double zmn, double zmx);

  void       GetXZMinMax(double &xmn, double &xmx, double &zmn, double &zmx);

  bool       CheckSplineStatus();

protected:
 
  virtual int IntervalSmoothParSpline(int N1,double *x1,double *z1,bool *flag,int N2, double *xf,double *zf,
                      double rc, double sp, double_ar *xx1, double_ar *zz1, double_ar *dx, double_ar *dz,
                      double_ar *curvity,int_ar *NumNodeBreak);
  
private:
  
  bool       GetActiveKnotNs(double x1, double z1, int &n1, int &n2);

  bool       GetKnots(int n1, int n2, double_ar *xx, double_ar *zz);

  bool       PrepareKnots(double x1, double z1, double x2, double z2, double_ar *xx, double_ar *zz, bool fp=true);

  bool       CheckPPos(double x1, double z1, double x2, double z2);

  bool       CheckPPosX(double x1, double x2);

  bool       CheckPPosZ(double z1, double z2);
    
  double_ar    x, z, dx, dz, curv;
  double_ar    px, pz;
  double_ar    sx, sz;

  //AlgData3d   ald;
  //Gauss3dData gad;

  double_ar    resXs;
  double_ar    resZs;
  int_ar       NumNode;
  bool*        s_p;
  bool*        sPoints;

  double      xMin, xMax, zMin, zMax;

  bool        clean;

  int         n;
  bool        ready;

  ParSplineOneD *pSp;
};



//=========================================================================
#endif