BaseFunct3D.h

// =======================================================================
// %W% %G%
//
// Class:      TBaseFunct3D
//
// Purpose:    represents the set of base functions for a finite element
//             in three dimensions
//
// Author:     Gunar Matthies
//
// History:    start of implementation 19.11.99
// 
// =======================================================================

#ifndef __BASEFUNCT3D__
#define __BASEFUNCT3D__

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

#include <Enumerations.h>

class TBaseFunct3D
{
  protected:
    int Dimension;

    BaseFunct3D BaseFunct;

    DoubleFunct3D *Functions[N_MultiIndices3D];

    bool changable;

    BF3DRefElements RefElement;

    int PolynomialDegree;

    int Accuracy;

    int N_BF2Change;

    int ***BF2Change;
    
    int BaseVectDim;

  public:
    TBaseFunct3D(int dimension, BaseFunct3D basefunct,
                 BF3DRefElements refelement,
                 DoubleFunct3D* functions, 
                 DoubleFunct3D* derivativesXi,
                 DoubleFunct3D* derivativesEta,
                 DoubleFunct3D* derivativesZeta,
                 DoubleFunct3D* derivativesXiXi,
                 DoubleFunct3D* derivativesXiEta,
                 DoubleFunct3D* derivativesXiZeta,
                 DoubleFunct3D* derivativesEtaEta,
                 DoubleFunct3D* derivativesEtaZeta,
                 DoubleFunct3D* derivativesZetaZeta,
                 int polynomialdegree,
                 int accuracy,
                 int n_bf2change,
                 int ***bf2change,
                 int baseVectDim = 1
                );

    TBaseFunct3D(int dimension);

    int GetDimension() const
    { return Dimension; }
    
    void GetDerivatives(MultiIndex3D MultiIndex, double xi,
                        double eta, double zeta, double *values)
    { Functions[MultiIndex](xi, eta, zeta, values); };

    void GetDerivatives(MultiIndex3D MultiIndex, 
                        TQuadFormula3D *formula, double **values);

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

   void GetValues(int N_Points, double *t, double *s, int i, 
                  MultiIndex3D index, double **Values);

    void SetUnchangable()
      { changable = false; };

    void SetFunction(MultiIndex3D MultiIndex, DoubleFunct3D* function);

    void MakeRefElementData(QuadFormula2D QuadFormula);

    void MakeRefElementData(QuadFormula3D QuadFormula);

    TGridCell *GenerateRefElement();

    BF3DRefElements 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);

    BaseFunct3D GetID() const
    { return BaseFunct; }
    
    int GetBaseVectDim() const
    { return BaseVectDim; }
    
};

#endif