AssembleMat3D.h

#ifndef __ASSEMBLEMAT3D__
#define __ASSEMBLEMAT3D__

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

#ifdef _MPI

#endif

#ifdef _OMPONLY

#endif

class TAssembleMat3D: public TAssemble
{
  protected:  

  TFESpace3D **FeSpaces,  **FeRhs;
  TSquareMatrix3D **SqMatrices; 
  TMatrix3D **RecMatrices;
  TDiscreteForm3D *DiscreteForm;
  BoundCondFunct3D **BoundaryConditions;
  BoundValueFunct3D **BoundaryValues;
  TAuxParam3D *AuxParam;                     
  double *Param[MaxN_QuadPoints_3D], *AuxArray[MaxN_QuadPoints_3D];  

 
  public:
    TAssembleMat3D(int n_fespaces, TFESpace3D **fespaces,
                      int n_sqmatrices, TSquareMatrix3D **sqmatrices,
                      int n_matrices, TMatrix3D **matrices,
                      int n_rhs, double **rhs, TFESpace3D **ferhs,
                      TDiscreteForm3D *discreteform,
                      BoundCondFunct3D **boundarybonditions,
                      BoundValueFunct3D **boundaryvalues,
                      TAuxParam3D *parameters);
    
    ~TAssembleMat3D();
    
  void Init();
    
  void DeAllocate();

  void Reset();
  
  void Assemble3D(); 
  
  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, BaseFunct3D *BaseFuncts, RefTrans3D reftrans);
  
  void ModifyMatHang();

  void PrintAllMat();
  
};

#endif