BdNonUniformSpline.h

// =======================================================================
// @(#)BdNonUniformSpline.h        1.2 07/16/99
//
// Class:       TBdNonUniformSpline
// Superclass:  TBoundComp
// Purpose:     spline function as a component of a boundary part
//
// Author:      Volker Behns  18.06.97
//
// =======================================================================

#ifndef __BDNONUNIFORMSPLINE__
#define __BDNONUNIFORMSPLINE__

#include <BoundComp2D.h>

/*    determined by 10 parameters;
    every subspline is presented as X(t), Y(t), t=0..1,
    X(t) = Params[0]*phi1(t) + Params[2]*phi2(t) +
           Params[4]*phi3(t) + Params[6]*phi4(t)
    Y(t) = Params[1]*phi1(t) + Params[3]*phi2(t) +
           Params[5]*phi3(t) + Params[7]*phi4(t); 
    phi1(0) = 1, phi2(0) = 0, phi3(0) = 0, phi4(0) = 0;
    phi1(1) = 0, phi2(1) = 1, phi3(1) = 0, phi4(1) = 0;
    d_phi1(0) = 0, d_phi2(0) = 0, d_phi3(0) = 1, d_phi4(0) = 0;
    d_phi1(1) = 0, d_phi2(1) = 0, d_phi3(1) = 0, d_phi4(1) = 1; 
 Params[8] = T[i], T=0..1 through the whole spline;
 Params[9] is not used */

class TBdNonUniformSpline : public TBoundComp2D
{
  protected:
    int N_Splines;
    double *Params;
    double *Param9;

  public:
    // Constuctor
    TBdNonUniformSpline (int id, int N_Splines);
    // Destructor
    ~TBdNonUniformSpline () {delete Params; delete Param9;};

    // Methods
    void SetParams (double *params);
    int GetN_Splines();
    //double *GetBoundPoints(QuadFormula *formula, int &N_NewPt);

    virtual int GetXYofT(double T, double &X, double &Y);

    virtual int GetTofXY(double X, double Y, double &T);

    virtual int ReadIn(std::ifstream &dat);

    virtual int GetN_InitVerts()
    { return 4; }
    virtual int GenInitVerts(double *&points, int I_points,
                             int *&edges, int I_edges)
    { return -1; }
    
    double GetLocalXofT(int ISpline, double LocalT);
    
    double GetLocalYofT(int ISpline, double LocalT);
    
    void GenerateParams1(double *x, double *y, double dx0, double dy0, double dx1, double dy1);

    void Solver_3diag(double *a, double *b, double *c, double *rhs, double *sol);

    double AdXofT(int ISpline, double LocalT);
    double AdYofT(int ISpline, double LocalT);
    double GetLocalDXofT(int ISpline, double LocalT);
    double GetLocalDYofT(int ISpline, double LocalT);
};

#endif