Assemble2D.h

// =======================================================================
// @(#)Assemble2D.h        1.5 04/13/00
// 
// Purpose:     assemble matrix and right-hand side
//
// Author:      Gunar Matthies (10.08.98)
//
// History:     start of implementation 10.08.98 (Gunar Matthies)
//
// =======================================================================

#ifndef __ASSEMBLE2D__
#define __ASSEMBLE2D__

#include <AllClasses.h>
#include <Constants.h>
#include <FEDatabase2D.h>
#include <LocalAssembling2D.h>

#ifdef __3D__
  #include <Aux2D3D.h>
#endif

void Assemble2D(int n_fespaces, TFESpace2D **fespaces,
                int n_sqmatrices, TSquareMatrix2D **sqmatrices,
                int n_matrices, TMatrix2D **matrices,
                int n_rhs, double **rhs, TFESpace2D **ferhs,
                BoundCondFunct2D **BoundaryConditions,
                BoundValueFunct2D **BoundaryValues,
                LocalAssembling2D& la
#ifdef __3D__
                , TAux2D3D *Aux2D3D
#endif
, int AssemblePhaseID = -1 
               );





void Assemble2D(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 **BoundaryConditions,
                BoundValueFunct2D **BoundaryValues,
                TAuxParam2D *parameters
#ifdef __3D__
                , TAux2D3D *Aux2D3D
#endif
, int AssemblePhaseID = -1 
               );

void Assemble2D_FCT(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 **BoundaryConditions,
                    BoundValueFunct2D **BoundaryValues, TAuxParam2D *Parameters,
                    double factor
#ifdef __3D__
                    , TAux2D3D *Aux2D3D
#endif
                   );



void Assemble2DSlipBC(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 **BoundaryConditions,
                      BoundValueFunct2D **BoundaryValues,
                      TAuxParam2D *parameters,
                      TFEFunction2D *u1, TFEFunction2D *u2);

void Assemble2D_neigh(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 **BoundaryConditions,
                BoundValueFunct2D **BoundaryValues,
                TAuxParam2D *Parameters
#ifdef __3D__
                , TAux2D3D *Aux2D3D
#endif
                );

void Assemble2D_DG(CoeffFct2D *Coeff, int n_fespaces, TFESpace2D **fespaces,
           int n_sqmatrices, TSquareMatrix2D **sqmatrices,
           int n_matrices, TMatrix2D **matrices,
           int n_rhs, double **rhs, TFESpace2D **ferhs,
           BoundCondFunct2D **BoundaryConditions,
           BoundValueFunct2D **BoundaryValues,
           TAuxParam2D *Parameters);

void Assemble2D_CIP(CoeffFct2D *Coeff,int n_fespaces, TFESpace2D **fespaces,
            int n_sqmatrices, TSquareMatrix2D **sqmatrices,
            int n_matrices, TMatrix2D **matrices,
            int n_rhs, double **rhs, TFESpace2D **ferhs,
            BoundCondFunct2D **BoundaryConditions,
            BoundValueFunct2D **BoundaryValues,
            TAuxParam2D *Parameters);

void Assemble2D_VectFE(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 **BoundaryConditions,
           BoundValueFunct2D **BoundaryValues,
           TAuxParam2D *Parameters=NULL
           );

#ifdef __MORTAR__
void Assemble(TMatrix2D *matrix);

  #ifdef __ADD_LINK_SDFEM__
  void AddLinkSDFEM(TFESpace2D *Space2D, TSquareMatrix2D *Matrix,
  #endif
#endif // __MORTAR__

void Assemble2D(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 **BoundaryConditions,
                BoundValueFunct2D **BoundaryValues,
                TAuxParam2D *Parameters,
        TAuxParam2D *ParametersBound, 
        TypeBoundSwitchFunct2D *TypeBoundSwitcher,
        int *CounterBoundaryParam
#ifdef __3D__
                , TAux2D3D *Aux2D3D
#endif
                );
#ifdef __2D__
void ComputeInterfaceConditionDirichlet(TSquareMatrix2D  **sqmatrices,
                    TMatrix2D **matrices,
                    TFEVectFunct2D *u_other,
                    double *rhs,
                    int *interface_dof_Array,
                    int *interface_dof_other_Array,
                    int interface_dof, double density, 
                    double density_other);


void ComputeInterfaceConditionStress(TSquareMatrix2D  **sqmatrices,
                     TMatrix2D **matrices,
                     TFEVectFunct2D *u_other,
                     double *rhs,
                     int *interface_dof_Array,
                     int *interface_dof_other_Array,
                     int interface_dof, double density, 
                     double density_other, 
                     double *rhs_for_stress_other);
     
#endif // __2D__

#endif // __ASSEMBLE2D__