RefTrans2D.h

// =======================================================================
// @(#)RefTrans2D.h        1.4 04/13/00
//
// Class:      TRefTrans2D
//
// Purpose:    reference transformations for 2D geometric objects
//
// Author:     Gunar Matthies
//
// History:    08.07.97 start implementation
// 
// =======================================================================

#ifndef __REFTRANS2D__
#define __REFTRANS2D__

#include <Constants.h>
#include <Enumerations.h>
#include <BaseCell.h>

class TRefTrans2D
{
  protected:
    TBaseCell *Cell;

  public:
    TRefTrans2D() {};

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

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

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

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

    void GetOrigValues(TBaseCell *cell);

    virtual void SetCell(TBaseCell *cell)
    {  Cell = cell; }

    static RefTrans2D FindRefTrans2D
        (int N_LocalUsedElements, FE2D *LocalUsedElements);

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

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

    double GetVolume();
    
    // piola map, needed for vector values basis functions such as 
    // Raviart-Thomas (RT) or Brezzi-Douglas-Marini (BDM)
    virtual void PiolaMapOrigFromRef(int N_Functs, double *refD00, double *origD00 ) 
    { cout << " Piola Map not defined for this element " << endl; };
    
    // piola map for derivatives
    virtual void PiolaMapOrigFromRef(int N_Functs, double *refD10, 
                                     double *refD01, double *origD10, 
                                     double *origD01)
    { cout << " Piola Map not defined for this element " << endl; };
};

#endif