SystemDarcy2D.h

#ifndef __SYSTEMDARCY2D__
#define __SYSTEMDARCY2D__

#include <SquareMatrix2D.h>
#include <LocalAssembling2D.h>

class TSystemDarcy2D
{
  protected:
    
    TFESpace2D *fe_spaces[2];
    
    int N_Matrices;

    // ( A  B1' )   ( 0 2 )
    // ( B2 C   )   ( 3 1 )
    
    TSquareMatrix2D *sq_matrices[2];
    
    TMatrix2D *rect_matrices[2];
    
    BoundCondFunct2D *BoundaryConditions[2];

    BoundValueFunct2D *BoundaryValues[2];
    
  public:
     TSystemDarcy2D(TFESpace2D **fespaces);

    ~TSystemDarcy2D();
    
    void Init(BoundCondFunct2D **BoundCond, BoundValueFunct2D **BoundValue);
 
    void Assemble(LocalAssembling2D& la, double *sol, double *rhs);

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

#endif // __SYSTEMMATDARCY2D__