parspline_service.hpp

// parspline_service.hpp
// $Id: parspline_service.hpp,v 1.4 2007/08/05 09:50:35 guser1 Exp $
//---------------------------------------------------------------------------------------------------------
#ifndef ParSplineOneDH
#define ParSplineOneDH
#include <math.h>
#include <string.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
//#include <mix/wenv.h>
//#include <mix/DoubleAr.H>
//#include <mix/IntAr.H>


#include "dvoev_ma.hpp"
#include "objectar.hpp"
#include "s_message.hpp"

//---------------------------------------------------------------------------------------------------------
struct AlgData
{
  double  **MN;
  double  *P;
  double  *abc;
};

struct AlgData3d
{
  double MN[3][3];
  double P[3];
  double abc[3];
};

struct Gauss3dData
{
  double MN[3][3];
  double P[3];
  int    index1[3];
  int    index2[3];
};

//---------------------------------------------------------------------------------------------------------
//bool Gauss (int N, double **LN, double *P, double *x);
bool Gauss (int N, double **LN, double *P, double *x);
bool Gauss (int N, AlgData *alg);
bool PointLineIntersection(double x01,double z01,double x11,double z11,double x02,double z02,
               double x12,double z12,double &xt,double &zt,bool flag1,bool flag2);

//int TwoParSplineIntersection(int N1,double_ar *xx1,double_ar *zz1,int N2,double_ar *xx2,double_ar *zz2,
//double_ar *xp,double_ar *zp);

//int TwoParSplineIntersection(int N1,double_ar *xx1,double_ar *zz1,int N2,double_ar *xx2,double_ar *zz2,
//           double_ar *xp,double_ar *zp,bool *flag1,bool *flag2);

int LineParSplineIntersection(double x0,double z0,double xk,double zk,int N1,double_ar *xx1,double_ar *zz1,
                  double_ar *xp,double_ar *zp,bool flag,int_ar *NumNode,bool *flagd);
//, AlgData3d *alg, Gauss3dData *gad);

bool LineParSplineIntersection1(double x0,double z0,double xk,double zk,int N1,double_ar *xx1,double_ar *zz1,
                double &xp,double &zp,double &beta,double &curvity,bool flag,int NumNode,
                double_ar *Dx,double_ar *Dz,double_ar *curv,bool *flagd);
//, AlgData3d *alg, Gauss3dData *gad);
//---------------------------------------------------------------------------------------------------------
class ParSplineOneD
{
  double_ar  *a1x,*a2x,*a3x,*a4x,*a1z,*a2z,*a3z,*a4z,*G,*PHIX,*PHIZ;
  double    dphin3;
  int       N,err,NI;
  double    sigma,amax,MinRc;
  char      *strerr;

  void      DateEditorParSplineOneD();
  void      CubParSplineOneD(double_ar *e,double_ar *a,double_ar *b,double_ar *c,double_ar *d,int ind);
  void      ApproxCubSplineOneD(double_ar *e,double_ar *a,double_ar *b,double_ar *c,double_ar *d,int ind);
  double    MinRadCurvityForFixPoint(int N2,double *x1,double *x2);
  void      InsertNode(int &N0,int k0,int *num,double xt,double zt);
  void      SigmaCalculation();
  void      AmaxCalculation();
  void      InitParameters(double_ar *a,double_ar *b);
  double    Angle(double x01,double x11,double z01,double z11,double x02,double x12,double z02,double z12);

 public:

  double_ar  *S;
  double_ar  *xx,*zz,*x,*z,*phinorm,*curv,*Dx,*Dz;
  ParSplineOneD();
  ParSplineOneD(int N,double *x,double *z,int N1,double *x1,double *z1,double rc);
  ~ParSplineOneD();

  // void NewParSpline(int N,double *x,double *z,int N1,double *x1,double *z1,double rc);
  void      GetPoints(int N,double *x,double *z,int flag);
  int       GetN();
  void      SetN(int n);
  void      GetXZ(double *x1,double *z1);

  void      GetABCD(double *a1x,double *b1x,double *c1x,double *d1x, double *a1z,double *b1z,double *c1z,
            double *d1z);

  double    GetSigma();
  double    GetAmax();
  double    GetMinRc();
  void      InterpolCurve(double g,int &n);
  bool      ParSplineFunction(double s,double &xt,double &zt);
  bool      DerParSplineFunction(double s,double &dx,double &dz);
  bool      Der2ParSplineFunction(double s,double &d2x,double &d2z);
  bool      CurvityParSplineFunction(double s,double &cr);
  bool      RadiusCurvityParSplineFunction(double s,double &rcr);
};
//---------------------------------------------------------------------------------------------------------


#endif