SystemCD2D.h

#ifndef __SYSTEMCD2D__
#define __SYSTEMCD2D__

#include <SquareMatrix2D.h>

class TSystemCD2D
{
  protected:
    
    TFESpace2D *FeSpace;

    int N_DOF, N_Active;
    
    int Disctype;
    
    int SOLVER;
       
    int N_Matrices;

    TSquareStructure2D *sqstructure;

    TSquareMatrix2D *sqmatrixA, *SQMATRICES[3];;
    TSquareMatrix **sqmatrices;
    
    BoundCondFunct2D *BoundaryConditions[1];

    BoundValueFunct2D *BoundaryValues[1];
        
    TDiscreteForm2D *DiscreteFormARhs;
 
    
  public:
     TSystemCD2D(TFESpace2D *fespace, int disctype, int solver);

    ~TSystemCD2D();
    
    void Init(CoeffFct2D *BilinearCoeffs, BoundCondFunct2D *BoundCond, BoundValueFunct2D *BoundValue);
 
    void Assemble(TAuxParam2D *aux, double *sol, double *rhs);

    void  Solve(double *sol, double *rhs);
    
};

#endif