FE2DMapper.h

// =======================================================================
// @(#)FE2DMapper.h        1.2 04/12/99
//
// Class:       TFE2DMapper
// Purpose:     find out which of the given local degress of freedom
//              are equivalent to the same global degree of freedom
//
// Author:      Gunar Matthies  05.11.97
//
// =======================================================================

#ifndef __FE2DMAPPER__
#define __FE2DMAPPER__

#define HANGINGNODE -5
#define FIRSTMARK -10

#include <FEDesc2D.h>
#include "Vector.h"
#include <HangingNode.h>

class TFE2DMapper
{
  protected:
    char *Name;

    char *Description;

    int N_DOF0;

    int N_DOF1;

    int N_Pairs;

    int *Pairs;

    int N_NoOpposite;

    int *NoOpposite;

    int *Aux;

// =======================================================================
// local hanging nodes = coupling only this node on this edge
// =======================================================================
    int N_Hanging;

    int *Hanging;

    HNDesc *HangingTypes;

    int **Coupling;

// =======================================================================
// far hanging nodes = coupling also with node not on this edge
// =======================================================================
    int N_FarHanging;

    int *FarHanging;

    HNDesc *FarHangingTypes;

    int ****FarCoupling;

  public:
    TFE2DMapper(char *name, char *description, int n0, int n1, 
              int n_pairs, int *pairs, 
              int n_hanging, int *hanging,
              HNDesc *hangingtypes, int **coupling,
              int n_farhanging, int *farhanging,
              HNDesc *farhangingtypes, int ****farcoupling,
              int n_noopposite, int *noopposite,
              int n_nodes);

    // Methods
    char *GetName()
    { return Name; }

    char *GetDescription()
    { return Description; }

    int GetN_DOF0()
    { return N_DOF0; }

    int GetN_DOF1()
    { return N_DOF1; }

    void GetN_DOF(int &n0, int &n1)
    {
      n0 = N_DOF1;
      n1 = N_DOF1;
    }

    int GetN_Pairs()
    { return N_Pairs; }

    int *GetPairs()
    { return Pairs; }

    HNDesc *GetHangingTypes()
    { return HangingTypes; }

    int **GetCoupling()
    { return Coupling; }

    HNDesc *GetFarHangingTypes()
    { return FarHangingTypes; }

    int ****GetFarCoupling()
    { return FarCoupling; }

    int GetN_NoOpposite()
    { return N_NoOpposite; }

    int *GetNoOpposite()
    { return NoOpposite; }

    void Map(int *Global, int I_K0, int I_K1, 
             int *Indices0, int *Indices1,
             int LocEdge0, int LocEdge1,
             TFEDesc2D *Desc0, TFEDesc2D *Desc1,
             int &Counter,
             TVector<THangingNode *> *vect,
             TVector<int> *numbers);

    void MapBound(int *Global, int I_K, int *Indices, 
                  int &BoundCounter);

    void MapDOF(int *Global, int dof0, int dof1, int &Counter);

    ~TFE2DMapper()
    { delete Aux; }

};

#endif