NSE2DErrorEstimator.h

// =======================================================================
// @(#)NS2DErrorEstimator.h        1.3 11/15/99
// 
// Class:       TNS2DErrorEstimator
//
// Purpose:     define error estimator
//
// Author:      Volker John
//
// History:     29.01.98 start implementation
//
// =======================================================================
#ifdef __2D__

#ifndef __NS_ERROR_ESTIMATOR__
#define __NS_ERROR_ESTIMATOR__

#include <FEFunction2D.h>
#include <FEVectFunct2D.h>

#define ns_gradient_indicator                      0
#define ns_residual_estimator_h1                   1
#define ns_residual_estimator_l2                   2
#define ns_residual_estimator_energy_quasi_robust  3
#define ns_gradient_recovery                       4
#define ns_implicit_estimator_neumann              5

class TNS2DErrorEstimator
{
  protected:
    TCollection *Collection_U;
    TCollection *Collection_P;
    
    TFESpace2D *FESpace2D_U;
    TFESpace2D *FESpace2D_P;

    TFEVectFunct2D *U;
    TFEFunction2D *P;

    int FELocalEstimator;

    int ErrorControl;
 
    int NavierStokes;

public:
    TNS2DErrorEstimator(int fe_local_estimator, 
                      TFEVectFunct2D *u, 
                      TFEFunction2D *p, 
                      int error_control,
                      int navierstokes);

    TCollection *GetCollection_U()
      { return Collection_U; };
    TCollection *GetCollection_P()
      { return Collection_P; };

    TFEVectFunct2D *GetU()
      { return U;};

    TFEFunction2D *GetP()
      { return P;};

    int GetFELocalEstimator()
      { return FELocalEstimator;};

    int GetErrorControl()
      { return ErrorControl;}

    int GetNavierStokes()
      { return NavierStokes;}

    void GetErrorEstimate(int N_Derivatives,
                          MultiIndex2D *NeededDerivatives,
                          int N_DerivativesP,
                          MultiIndex2D *NeededDerivativesP,
                          CoeffFct2D *Coeff, 
                          BoundCondFunct2D **BoundaryConds,
                          BoundValueFunct2D **BoundaryValues,
                          TAuxParam2D *Aux,
                          int n_fespaces,
                          TFESpace2D **fespaces,
                          double *eta_K, 
                          double *eta_max,
                          double *estimated_global_error);
      
    void  EstimateCellError(TFESpace2D **fespaces,
                            TBaseCell *cell,
                            int N_Points, 
                            double *X, 
                            double *Y, 
                            double *AbsDetjk, 
                            double *weights,
                            double **Derivatives, 
                            double **AuxArray, 
                            BoundCondFunct2D **BoundaryConds,
                            BoundValueFunct2D **BoundaryValues,
                            int N_Points1D, 
                            double *zeta,
                            double *X1D[4], 
                            double *Y1D[4], 
                            double *weights1D,
                            double *xyval_ref1D[4],
                            double *xderiv_ref1D[4],
                            double *yderiv_ref1D[4],
                            int *GlobalNumbers, 
                            int *BeginIndex,
                            int *DOF,
                            double *Values,                                           
                            int *GlobalNumbersP, 
                            int *BeginIndexP,
                            int *DOFP,
                            double *ValuesP,
                            double *local_error);

};

#endif

#endif // #ifdef __2D__