FE3DMapper.h

// =======================================================================
// %W% %G%
//
// Class:       TFE3DMapper
// 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
//
// History:     Added edge details by Sashi (08.03.12)
// =======================================================================

#ifndef __FE3DMAPPER__
#define __FE3DMAPPER__

#define HANGINGNODE -5
#define FIRSTMARK -10

#include <FEDesc3D.h>
#include "Vector.h"

class TFE3DMapper
{
  protected:
    char *Name;

    char *Description;

    int N_DOF0;

    int N_DOF1;

    int N_Pairs;

    int **Pairs;

    int N_NoOpposite;

    int **NoOpposite;

    int *Aux;

  public:
    TFE3DMapper(char *name, char *description, int n0, int n1, 
                int n_pairs, int **pairs, 
                int n_noopposite, int **noopposite,
                int n_nodes);

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

    char *GetDescription() const
    { return Description; }

    int GetN_DOF0() const
    { return N_DOF0; }

    int GetN_DOF1() const
    { return N_DOF1; }

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

    int GetN_Pairs() const
    { return N_Pairs; }

    int **GetPairs() const
    { return Pairs; }

    int *GetPairs(int i) const
    { return Pairs[i]; }

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

    void MapBound(int *Global, int I_K, int *Indices, 
                  int &BoundCounter,
                  TVector<THangingNode *> *vect,
                  TVector<int> *numbers);

    void MapDOF(int *Global, int dof0, int dof1, int &Counter,
                TVector<THangingNode *> *vect,
                TVector<int> *numbers);

    void MapBoundEdge(int N_EdgeDOF, int *Global, int I_K, int *Indices, 
                      int &BoundCounter,
                      TVector<THangingNode *> *vect,
                      TVector<int> *numbers);   
    
    void MapBoundVert(int *Global, int I_K, int Index, 
                     int &BoundCounter,
                     TVector<THangingNode *> *vect,
                     TVector<int> *numbers);  
    
    
    ~TFE3DMapper()
    { delete Aux; }

};

#endif