TetraIsoparametric.h

// =======================================================================
// @(#)TetraIsoparametric.h        1.3 02/22/00
//
// Class:      TTetraIsoparametric
//
// Purpose:    Isoparametric reference transformations for Tetrahedron
//
// Author:     Gunar Matthies
//
// History:    2000/11/20 start implementation
// 
// =======================================================================

#ifndef __TetraIsoparametric__
#define __TetraIsoparametric__

#include <Enumerations.h>
#include <RefTrans3D.h>

class TTetraIsoparametric : public TRefTrans3D
{
  protected:
    double x0, x1, x2, x3;

    double y0, y1, y2, y3;

    double z0, z1, z2, z3;

    double xc0, xc1, xc2, xc3;

    double yc0, yc1, yc2, yc3;

    double zc0, zc1, zc2, zc3;

    int N_AuxPoints;

    double XDistance[MaxN_BaseFunctions3D];

    double YDistance[MaxN_BaseFunctions3D];

    double ZDistance[MaxN_BaseFunctions3D];

    int ApproximationOrder;

    double FctValues[MaxN_QuadPoints_3D][MaxN_BaseFunctions3D];

    double XiDerValues[MaxN_QuadPoints_3D][MaxN_BaseFunctions3D];

    double EtaDerValues[MaxN_QuadPoints_3D][MaxN_BaseFunctions3D];

    double ZetaDerValues[MaxN_QuadPoints_3D][MaxN_BaseFunctions3D];

    static BaseFunct3D BaseFunctFromOrder[]; 

    static FEDesc3D FEDescFromOrder[]; 

    double DoubleAux[MaxN_BaseFunctions3D];

    int IntAux[MaxN_BaseFunctions3D];

    QuadFormula3D QuadFormula;

    double *XI, *ETA, *ZETA, *W;

    int N_QuadPoints;

    double detjk;

    double rec_detjk;

  public:
    TTetraIsoparametric();

    void GetOrigBoundFromRef(int Joint, double xi, double eta,
                            double &X, double &Y, double &Z);
    
    void GetOrigBoundFromRef(int Joint, int N_Points, double *xi, double *eta,
                            double *X, double *Y, double *Z);


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

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

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

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

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

    void GetOrigValues(BaseFunct3D BaseFunct,
                       int N_Points, double *xi, double *eta, double *zeta,
                       int N_Functs, QuadFormula3D QuadFormula);

    void GetOrigValues(int N_Sets, BaseFunct3D *BaseFunct,
                       int N_Points, double *xi, double *eta, double *zeta,
                       QuadFormula3D QuadFormula,
                       bool *Needs2ndDer);

    void GetOrigValues(int JointNr, double p1, double p2,
                       int N_BaseFunct,
                       double *uref, double *uxiref, double *uetaref, double *uzetaref,
                       double *uorig, double *uxorig, double *uyorig, double *uzorig);

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

    void SetCell(TBaseCell * cell);

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

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

    void GetOuterNormal(int j, double s, double t,
                        double &n1, double &n2, double &n3);

    void GetTangentVectors(int j, double p1, double p2,
        double &t11, double &t12, double &t13,
        double &t21, double &t22, double &t23);

};

#endif