CDSystemTimeDG.h

#ifndef __CDSYSTEMTIMEDG__
#define __CDSYSTEMTIMEDG__

#include <SquareMatrix2D.h>
#include <string.h>

class TCDSystemTimeDG
{
  protected:
    int N_Eqn;   
    
    int N_U, N_UActive;   
    
    double *Sol_dG, *rhs_dG;
    
    int N_RowBlockMat, N_ColBlockMat;  
    
    int Sys_N_Entries, *Sys_rowptr, *Sys_colindex;    

    TSquareMatrix2D *Mat_M, *Mat_A, *Sys_Mat;

    TSquareMatrix2D *Mat_M_Qp1, *Mat_A_Qp1; // dG(1)
    double *Rhs_Qp1;
    bool CONSERVATIVEALE, QpMatricsAdded;
     
    TSquareStructure2D *Sys_structure;
    
    TFESpace2D *RhsSpace;
    
    CoeffFct2D *BilinearCoeffs;
    
  public:
     TCDSystemTimeDG(int n_RowBlockMat, int n_ColBlockMat, TSquareMatrix2D *mat_M, TSquareMatrix2D *mat_A);

    ~TCDSystemTimeDG();
    
    void AddBilinear(CoeffFct2D *bilinearCoeffs)
          {BilinearCoeffs =  bilinearCoeffs;}

    void AddQp1Matrices(TSquareMatrix2D *mat_M_Qp1,TSquareMatrix2D *mat_A_Qp1, double *rhs_Qp1)
    {
     Mat_M_Qp1 = mat_M_Qp1; 
     Mat_A_Qp1 = mat_A_Qp1;
     Rhs_Qp1 = rhs_Qp1;
     QpMatricsAdded = TRUE;
    }
          
    virtual void AssembleSysMat(double *Mu_old, double *Rhs);
    
    virtual void AssembleALESysMat_Qp1(double *Mu_old, double *Rhs);
    
    void AssembleRhs(int N_Rhs, double *T, double *Rhs);

    void ResetSysMat()
    {  memset(Sys_Mat->GetEntries(), 0, Sys_N_Entries*SizeOfDouble);  }
    
    void SovedGSystem();
 
    virtual void SoveTimedG(double *Sol);   
    
    void SetALEForm(bool conservative)
    {  CONSERVATIVEALE = conservative; }
    
};

#endif