FEVectFunct2D.h

// =======================================================================
// @(#)FEVectFunct2D.h        1.2 07/20/99
// 
// Class:       TFEVectFunct2D
// Purpose:     a function from a finite element space in 2D
//
// Author:      Gunar Matthies (17.01.98)
//
// History:     start of implementation 17.01.98 (Gunar Matthies)
//
//              start of reimplementation 06.08.1998 (GM)
//
// =======================================================================


#ifndef __FEVECTFUNCT2D__
#define __FEVECTFUNCT2D__


#include <FEFunction2D.h>

class TFEVectFunct2D : public TFEFunction2D
{
  protected:
    int N_Components;

  public:
    TFEVectFunct2D(TFESpace2D *fespace2D, char *name, char *description,
                  double *values, int length, int n_components);

    int GetN_Components()
    { return N_Components; }

    TFEFunction2D *GetComponent(int i)
    {
      // the name of the component will include the index i
      std::ostringstream os;
      os.seekp(std::ios::beg);
      os << Name << i <<ends;
      return new TFEFunction2D(FESpace2D, (char*)os.str().c_str(), Description,
                               Values+i*Length, Length);
    }

    void GridToData();

    void DataToGrid();

    void GetDeformationTensorErrors( 
        DoubleFunct2D *Exact, DoubleFunct2D *Exact1,
        int N_Derivatives,
        MultiIndex2D *NeededDerivatives,
        int N_Errors, ErrorMethod2D *ErrorMeth, 
        CoeffFct2D *Coeff, TAuxParam2D *Aux,
        int n_fespaces, TFESpace2D **fespaces,
        double *errors);

    double GetL2NormDivergence();

   void WriteSol(double t);

   void ReadSol(char *BaseName);
   
   void FindVectGradient(double x, double y, double *val1, double *val2);

   void Interpolate(TFEVectFunct2D *OldVectFunct);
   
    TFEVectFunct2D& operator*=(double alpha);
   
    TFEVectFunct2D & operator+=(const TFEVectFunct2D & rhs);
   
    TFEVectFunct2D & operator=(const TFEVectFunct2D & rhs);
};

#endif

// #ifdef __2D__
// # endif // #ifdef __2D__