FESpace2D.h

// =======================================================================
// @(#)FESpace2D.h        1.5 06/13/00
// 
// Class:       TFESpace2D
// Purpose:     class for all 2D finite element spaces
//
// Author:      Gunar Matthies (03.11.97)
//
// History:     start of implementation 03.11.97 (Gunar Matthies)
//
//              split FESpace into TFESpacexD
//              15.04.1998 (Volker Behns)
//
//              start of reimplementation
//              30.07.98 (Gunar Matthies)
//
// =======================================================================

#ifndef __FESPACE2D__
#define __FESPACE2D__

#include <FESpace.h>
#include <FE2D.h>

class THangingNode;

class TFESpace2D : public TFESpace
{
  protected:
    int N_ActiveDegrees;

    int N_SlaveDegrees;

    int HangingBound;

    FE2D *UsedElements; 

    FE2D *ElementForShape;

    THangingNode **HangingNodeArray;

    FE2D *AllElements;
    
   BoundCondFunct2D *BoundCondition;

#ifdef __MORTAR__

    TCollection *MortarColl;
#endif

   int N_IntlPts;

   double *X_Intl;

   double *Y_Intl;

   int *IntlPtIndexOfPts;


  public:
    TFESpace2D(TCollection *coll, char *name, char *description);

    // if no mortar is used, the last argument can be set to NULL
    TFESpace2D(TCollection *coll, char *name, char *description, 
               BoundCondFunct2D *BoundaryCondition, int k,
               TCollection *mortarcoll);

    TFESpace2D(TCollection *coll, char *name, char *description, 
               BoundCondFunct2D *BoundaryCondition, SpaceType type,
               int k, TCollection *mortarcoll);

    TFESpace2D(TCollection *coll, char *name, char *description,
               BoundCondFunct2D *BoundaryCondition,
               FE2D *fes, TCollection *mortarcoll);

    ~TFESpace2D();

    void FindUsedElements();

    void ConstructSpace(BoundCondFunct2D *BoundaryCondition);

    int GetN_ActiveDegrees() const
    { return N_ActiveDegrees; }

    int GetN_SlaveDegrees() const
    { return N_SlaveDegrees; }

    int GetHangingBound() const
    { return HangingBound; }

    int GetN_Hanging() const
    { return N_SlaveDegrees; }

    FE2D *GetUsedElements() const
    { return UsedElements; }

    THangingNode **GetHangingNodes() const
    { return HangingNodeArray; }

    FE2D GetFE2D(int i, TBaseCell *cell);

    void GetDOFPosition(int dof, double &x, double &y);

    void GetDOFPosition(double *x, double *y);

    void SetIntlPtIndexOfPts(int *intlPtIndexOfPts)
     { IntlPtIndexOfPts = intlPtIndexOfPts; }

    int *GetIntlPtIndexOfPts() const
     { return IntlPtIndexOfPts; }
     
     FE2D *GetAllElements() const
     { return AllElements; }
     
    BoundCondFunct2D *GetBoundCondition() const
    { return BoundCondition; }
    
    friend  bool operator== (const TFESpace2D &lhs, const TFESpace2D &rhs);
    friend  bool operator!= (const TFESpace2D &lhs, const TFESpace2D &rhs);
};

#ifdef __MORTAR__
  double GetLambda(double , double , TVertex *, double , double );
#endif // __MORTAR__

#endif