TriaIsoparametric.h

// =======================================================================
// @(#)TriaIsoparametric.h        1.5 04/13/00
//
// Class:      TTriaIsoparametric
//
// Purpose:    isoparametric reference transformations for triangle
//
// Author:     Gunar Matthies
//
// History:    29.04.98 start implementation
// 
// =======================================================================

#ifndef __TRIAISOPARAMETRIC__
#define __TRIAISOPARAMETRIC__

#include <Enumerations.h>
#include <RefTrans2D.h>

class TTriaIsoparametric : public TRefTrans2D
{
  protected:
    double x[3];

    double y[3];

    double xc0, xc1, xc2;

    double yc0, yc1, yc2;

    int N_AuxPoints;

    double XDistance[MaxN_BaseFunctions2D];

    double YDistance[MaxN_BaseFunctions2D];

    int ApproximationOrder;

    double FctValues[MaxN_QuadPoints_2D][MaxN_BaseFunctions2D];

    double XiDerValues[MaxN_QuadPoints_2D][MaxN_BaseFunctions2D];

    double EtaDerValues[MaxN_QuadPoints_2D][MaxN_BaseFunctions2D];

    static BaseFunct2D BaseFunctFromOrder[]; 

    static FEDesc2D FEDescFromOrder[]; 

    double DoubleAux[MaxN_BaseFunctions2D];

    int IntAux[MaxN_BaseFunctions2D];

    QuadFormula2D QuadFormula;
  
    double *XI, *ETA, *W;

    int N_QuadPoints;

  public:
    TTriaIsoparametric();

    void GetOrigFromRef(double eta, double xi, double &x, double &y);

    void GetOrigFromRef(int N_Points, double *eta, double *xi, 
                        double *x, double *y, double *absdetjk);

    void GetOrigFromRef(double *ref, double *orig);

    void GetRefFromOrig(double x, double y, double &eta, double &xi);

    void GetRefFromOrig(double *orig, double *ref);

    void GetOrigValues(BaseFunct2D BaseFunct,
                       int N_Points, double *xi, double *eta,
                       int N_Functs, QuadFormula2D formula);

    void GetOrigValues(int N_Sets, BaseFunct2D *BaseFunct,
                       int N_Points, double *xi, double *eta,
                       QuadFormula2D formula,
                       bool *Needs2ndDer);

    void GetOrigValues(double xi, double eta, int N_BaseFunct,
                double *uref, double *uxiref, double *uetaref,
                double *uorig, double *uxorig, double *uyorig);

    void GetOrigValues(int joint, double zeta, int N_BaseFunct,
                double *uref, double *uxiref, double *uetaref,
                double *uorig, double *uxorig, double *uyorig);

    void SetCell(TBaseCell * cell);

    void SetApproximationOrder(int order)
    {
      if(order <= 0)
        ApproximationOrder = 1;
       else
        ApproximationOrder = order;
    }

    void SetQuadFormula(QuadFormula2D formula)
    { QuadFormula = formula; }

    void GetOuterNormal(int j, double zeta,
                                double &n1, double &n2);

    void GetTangent(int j, double zeta,
                                double &t1, double &t2);

    double GetVolume();

    void GetOrigBoundFromRef( int joint, int N_Points, double *zeta, double *X, double *Y);

};

#endif