AssembleMat2D.h

#ifndef __ASSEMBLEMAT2D__
#define __ASSEMBLEMAT2D__

#include <AllClasses.h>
#include <Assemble.h>
#include <Enumerations.h>

#ifdef _MPI

#endif

#ifdef _OMPONLY

#endif

class TAssembleMat2D: public TAssemble
{
  protected:  

  TFESpace2D **FeSpaces,  **FeRhs;
  TSquareMatrix2D **SqMatrices; 
  TMatrix2D **RecMatrices;
  TDiscreteForm2D *DiscreteForm;
  BoundCondFunct2D **BoundaryConditions;
  BoundValueFunct2D **BoundaryValues;
  TAuxParam2D *AuxParam;                     
  double *Param[MaxN_QuadPoints_2D], *AuxArray[MaxN_QuadPoints_2D];  

 
  public:
    TAssembleMat2D(int n_fespaces, TFESpace2D **fespaces,
                      int n_sqmatrices, TSquareMatrix2D **sqmatrices,
                      int n_matrices, TMatrix2D **matrices,
                      int n_rhs, double **rhs, TFESpace2D **ferhs,
                      TDiscreteForm2D *discreteform,
                      BoundCondFunct2D **boundarybonditions,
                      BoundValueFunct2D **boundaryvalues,
                      TAuxParam2D *parameters);
    
    ~TAssembleMat2D();
    
  void Init();
    
  void DeAllocate();

  void Reset();
  
  void Assemble2D(); 
  
  void AssembleNavierSlip();
  
private:  
  void AddLocalSqMatToGlobal(int i, TBaseCell *cell, int *N_BaseFunct);
 
  void AddLocalRecMatToGlobal(int i, TBaseCell *cell, int *N_BaseFunct);

  void AddLocalRhsToGlobal(int i, TBaseCell *cell, int *N_BaseFunct, BaseFunct2D *BaseFuncts, RefTrans2D reftrans);
  
  void ModifyMatHang();

  void PrintAllMat();
  
};

#endif