BaseFunct2D.h

// =======================================================================
// @(#)BaseFunct2D.h        1.5 02/08/00
//
// Class:      TBaseFunct2D
//
// Purpose:    represents the set of base functions for a finite element
//             in two dimensions
//
// Author:     Gunar Matthies
//
// History:    08.07.98 restart implementation
// //     :    10.06.2010 methods for vector basis functions (Sashikumaar Ganesan)
// =======================================================================

#ifndef __BASEFUNCT2D__
#define __BASEFUNCT2D__

#include <QuadFormula1D.h>
#include <QuadFormula2D.h>
#include <Constants.h>
#include <GridCell.h>

#include <Enumerations.h>

class TBaseFunct2D
{
  protected:
    int Dimension;

    BaseFunct2D BaseFunct;

    DoubleFunct2D *Functions[N_MultiIndices2D];

    bool changable;

    BF2DRefElements RefElement;

    int PolynomialDegree;

    int Accuracy;

    int N_BF2Change;

    int **BF2Change;    
 
    int BaseVectDim;
    
    bool SpaceDeptBasis;    
    
  public:
    TBaseFunct2D(int dimension, BaseFunct2D basefunct,
                 BF2DRefElements refelement,
                 DoubleFunct2D* functions, 
                 DoubleFunct2D* derivativesxi,
                 DoubleFunct2D* derivativeseta,
                 DoubleFunct2D* derivativesxixi,
                 DoubleFunct2D* derivativesxieta,
                 DoubleFunct2D* derivativesetaeta,
                 int polynomialdegree,
                 int accuracy,
                 int n_bf2change,
                 int **bf2change
                );

    TBaseFunct2D(int dimension, BaseFunct2D basefunct,
                 BF2DRefElements refelement,
                 DoubleFunct2D* functions, 
                 DoubleFunct2D* derivativesxi,
                 DoubleFunct2D* derivativeseta,
                 DoubleFunct2D* derivativesxixi,
                 DoubleFunct2D* derivativesxieta,
                 DoubleFunct2D* derivativesetaeta,
                 int polynomialdegree,
                 int accuracy,
                 int n_bf2change,
                 int **bf2change,
                 int baseVectDim
                );


    TBaseFunct2D(int dimension, BaseFunct2D basefunct,
                 BF2DRefElements refelement,
                 DoubleFunct2D* functions, 
                 DoubleFunct2D* derivativesxi,
                 DoubleFunct2D* derivativeseta,
                 DoubleFunct2D* derivativesxixi,
                 DoubleFunct2D* derivativesxieta,
                 DoubleFunct2D* derivativesetaeta,
                 int polynomialdegree,
                 int accuracy,
                 int n_bf2change,
                 int **bf2change,
                 bool spaceDeptBasis
                );



    TBaseFunct2D(int dimension);

    int GetDimension() const
    { return Dimension; }

    BaseFunct2D GetID() const
    { return BaseFunct; }

    void GetDerivatives(MultiIndex2D MultiIndex, double xi,
                        double eta, double *values)
    { Functions[MultiIndex](xi, eta, values); };

    void GetDerivatives(MultiIndex2D MultiIndex, 
                        TQuadFormula2D *formula, double **values);

    void GetValues(int N_Points, double *zeta, int i, double **Values);

    void GetDerivatives(MultiIndex2D MultiIndex, int N_Points,
                        double *zeta, int i, double **Values);

   void GetValues(int N_Points, double *zeta, int i, 
                  MultiIndex2D index, double **Values);

    void SetUnchangable()
      { changable = false; };

    void SetFunction(MultiIndex2D MultiIndex, DoubleFunct2D* function);

    void MakeRefElementData(QuadFormula1D QuadFormula);

    void MakeRefElementData(QuadFormula2D QuadFormula);

    TGridCell *GenerateRefElement();

    BF2DRefElements GetRefElement() const
      { return RefElement; };

    int GetPolynomialDegree() const
      { return PolynomialDegree; };

    int GetAccuracy() const
      { return Accuracy; };

    int GetN_BF2Change() const
      { return N_BF2Change; }

    int **GetBF2Change() const
      { return BF2Change; }

    void ChangeBF(TCollection *Coll, TBaseCell *Cell, double *Values);

    void ChangeBF(TCollection *Coll, TBaseCell *Cell, int N_Points, double **Values);
    
    int GetBaseVectDim() const
    { return BaseVectDim; }
    

};

#endif