QuadBilinear.h

// =======================================================================
// @(#)QuadBilinear.h        1.5 04/13/00
//
// Class:      TQuadBilinear
//
// Purpose:    Bilinear reference transformations for parallelogram
//
// Author:     Gunar Matthies
//
// History:    08.07.97 start implementation
// 
// =======================================================================

#ifndef __QUADBilinear__
#define __QUADBilinear__

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

class TQuadBilinear : public TRefTrans2D
{
  protected:
    double x0, x1, x2, x3;

    double y0, y1, y2, y3;

    double xc0, xc1, xc2, xc3;

    double yc0, yc1, yc2, yc3;

  public:
    TQuadBilinear();

    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 QuadFormula);

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

    void GetOrigValues(double xi, double eta, int N_BaseFunct,
                       double *uref, double *uxiref, double *uetaref,
                       double *uorig, double *uxorig, double *uyorig,
                       int _BaseVectDim = 1);
    
    void GetOrigValues(int joint, double zeta, int N_BaseFunct,
                       double *uref, double *uxiref, double *uetaref,
                       double *uorig, double *uxorig, double *uyorig,
                       int _BaseVectDim = 1);

    void SetCell(TBaseCell * cell);

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

    void GetTangent(int j, double zeta,
                                double &t1, double &t2);
    double GetVolume();
    
    void PiolaMapOrigFromRefNotAffine(int N_Functs, double *refD00, 
                                      double *origD00, double xi, double eta);
    void PiolaMapOrigFromRefNotAffine(int N_Functs, double *refD00, 
                                      double *refD10, double *refD01, 
                                      double *origD10, double *origD01,
                                      double xi, double eta);
};

#endif