SystemTNSE3D.h

#ifndef __SYSTEMTNSE3D__
#define __SYSTEMTNSE3D__

#include <SystemNSE3D.h>

class TSystemTNSE3D : public TSystemNSE3D
{
  protected:
    
    TFESpace3D **Projection_Spaces;
    
    TSquareMatrix3D **SqmatrixM11, **SqmatrixM12, **SqmatrixM13, **SqmatrixM21, **SqmatrixM22, **SqmatrixM23,
                    **SqmatrixM31, **SqmatrixM32, **SqmatrixM33;

#ifdef __PRIVATE__  

    TSquareStructure3D *sqstructureL;

    TStructure3D *structure_G, *structure_tilde_G; 
    
    TMatrix3D *matrix_tilde_G11, *matrix_tilde_G22, *matrix_tilde_G33, *matrix_G11, *matrix_G22, *matrix_G33;
    TMatrix3D **Matrices_tilde_G11, **Matrices_tilde_G22, **Matrices_tilde_G33, **Matrices_G11, **Matrices_G22, **Matrices_G33;  

    TSquareMatrix3D *sqmatrixL, **MatricesL;
#endif
     
    // Assembling rhs*/
    TAssembleMat3D *RhsOnlyAssemble;
    
    double *B;   
   
    double *defect;   
    
    double gamma;      

    TDiscreteForm3D *DiscreteFormRhs; 
    
    TAuxParam3D *NSE_Rhsaux;
    
    bool SystMatAssembled;
    
    double olderror_l_2_l_2u;
    
     
    
  public:
     TSystemTNSE3D(int N_levels, TFESpace3D **velocity_fespace, TFESpace3D **presssure_fespace, TFEVectFunct3D **velocity, 
                     TFEFunction3D **pressure, double **sol, double **rhs, int disctype, int nsetype, int solver
#ifdef __PRIVATE__  
                                   ,TFESpace3D **Projection_space
#endif    
                      );
     
    ~TSystemTNSE3D();

    void Init(CoeffFct3D *lincoeffs, BoundCondFunct3D *BoundCond, BoundValueFunct3D *U1BoundValue, 
              BoundValueFunct3D *U2BoundValue, BoundValueFunct3D *U3BoundValue);
    
    void Assemble();
        
    void AssembleRhs(); 
    
    void AssembleSystMat(double scale, double *oldrhs, double *rhs, double *sol);

    void AssembleSystMatNonLinear();
    
    void RestoreMassMat();
    
    void RestoreMassMatNonLinear();
      
    void AssembleNonLinear();

    void  Solve(double *sol);
    
    void GetResidual(double *sol, double &impuls_residual, double &residual);
 
    void MeasureTNSEErrors(DoubleFunct3D *ExactU1, DoubleFunct3D *ExactU2, DoubleFunct3D *ExactU3, DoubleFunct3D *ExactP, double *AllErrors);

    void printall_matrix();
    
    void GetMesh();
    
    void All_levels_check();
    
    double BlockMatVect(TSquareMatrix *A);
    
    double BlockMatVect(TMatrix *A);
    
};

#endif