SystemCST2D_Giesekus.h

 
#ifndef __SYSTEMCST2DGiesekus__
#define __SYSTEMCST2DGiesekus__

#include <SquareMatrix2D.h>

class TSystemCST2D_Giesekus
{
  protected:

    int N_S, N_Active, N_DirichletDof;
    
    TFESpace2D *FeSpace;
    
    TFESpace2D *FeSpaces_All[2];
    
    TFEFunction2D *FeFct[3];    
    
    TFEFunction2D *FeFct_All[5]; 
    
    int Disctype;
    
    int SOLVER;
    
    int N_Matrices;
           
    CoeffFct2D *LinCoeffs[1];    

    TSquareStructure2D *sqstructure;
  
    TSquareMatrix2D *SqmatrixS11, *SqmatrixS12, *SqmatrixS21, *SqmatrixS22, *SqmatrixS23, *SqmatrixS32, *SqmatrixS33, *SQMATRICES[8];
    TSquareMatrix **sqmatrices;
  
    TAuxParam2D *CSTNSEaux, *CSTaux_error;
     
    BoundCondFunct2D *BoundaryConditions[3];
  
    BoundValueFunct2D *BoundaryValues[3];
        
    TDiscreteForm2D *DiscreteFormARhs;

    
  public:
     TSystemCST2D_Giesekus(TFESpace2D *stress_fespace, TFEVectFunct2D *Stress, int disctype, int solver);

    ~TSystemCST2D_Giesekus();

    void Init(CoeffFct2D *lincoeffs, BoundCondFunct2D *BoundCond, BoundValueFunct2D *S1BoundValue, BoundValueFunct2D *S2BoundValue, BoundValueFunct2D *S3BoundValue, TAuxParam2D *aux, TAuxParam2D *auxerror);
 
    void Assemble(double *sol, double *rhs);
    
    void GetResidual(double *sol, double *rhs, double *res);

    void  Solve(double *sol, double *rhs);
    
    void MeasureErrors(DoubleFunct2D *ExactS1, DoubleFunct2D *ExactS2, DoubleFunct2D *ExactS3, double *s_error);

};

#endif