SystemCST2D.h

#ifndef __SYSTEMCST2D__
#define __SYSTEMCST2D__

#include <SquareMatrix2D.h>

class TSystemCST2D
{
  protected:

    int N_S, N_Active, N_DirichletDof;
    
    TFESpace2D *FeSpace[5];
      
    TFEFunction2D *FeFct[3];    

    // used for LPS streamline term
    TFEFunction2D *u1, *u2;
    
    int Disctype;
    
    int SOLVER;
    
    int Tensortype;
    
    int N_Matrices;
           
    CoeffFct2D *LinCoeffs[1];    

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

    
  public:
     TSystemCST2D(TFESpace2D *stress_fespace, TFEVectFunct2D *Stress, int tensortype, int disctype, int solver, TFESpace2D *Velocity_FeSpace, TFESpace2D* Pressure_FeSpace, TFESpace2D* Deformation_FeSpace = NULL, TFEVectFunct2D *Velocity = NULL);

    ~TSystemCST2D();

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

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

    
};

#endif