DiscreteForm3D.h

// =======================================================================
// %W% %G%
// 
// Class:       TDiscreteForm3D
// Purpose:     assemble on one cell a couple of bilinear and linear forms
//
// Author:      Gunar Matthies (02.10.98)
//
// History:     start of implementation 02.10.98 (Gunar Matthies)
//              add new data 14.04.99 (Gunar Matthies)
//
// =======================================================================

#ifndef __DISCRETEFORM3D__
#define __DISCRETEFORM3D__

#include <Enumerations.h>
#include <Constants.h>

class TDiscreteForm3D
{
  protected:
    char *Name;

    char *Description;

    int N_Terms;

    int N_Spaces;

    bool *Needs2ndDerivatives;

    MultiIndex3D *Derivatives;

    int *FESpaceNumber;

    int N_Matrices;

    int N_Rhs;

    int *RowSpace;

    int *ColumnSpace;

    int *RhsSpace;

    CoeffFct3D *Coeffs;

    AssembleFct3D *Assemble;

    AssembleFctParam3D *AssembleParam;

    ManipulateFct3D *Manipulate;

    double ***AllOrigValues;

    double **OrigValues;

  public:
    TDiscreteForm3D(char *name, char *description, int n_terms,
                    MultiIndex3D *derivatives, int *fespacenumber,
                    int n_matrices, int n_rhs,
                    int *rowspace, int *columnspace, int *rhsspace,
                    AssembleFct3D *assemble, CoeffFct3D *coeffs,
                    ManipulateFct3D *manipulate);

    TDiscreteForm3D(char *name, char *description, int n_terms,
                    MultiIndex3D *derivatives, int *fespacenumber,
                    int n_matrices, int n_rhs,
                    int *rowspace, int *columnspace, int *rhsspace,
                    AssembleFctParam3D *assembleparam, CoeffFct3D *coeffs,
                    ManipulateFct3D *manipulate);

    ~TDiscreteForm3D();

    char *GetName()
    { return Name; }

    char *GetDescription()
    { return Description; }

    void GetLocalForms(int N_Points, double *weights, double *AbsDetjk,
                        double hK, double *X, double *Y, double *Z,
                        int *N_BaseFuncts, BaseFunct3D *BaseFuncts, 
                        double **Parameters, double **AuxArray,
                        TBaseCell *Cell,
                        int N_Matrices, int N_Rhs,
                        double ***LocMatrix, double **LocRhs);

    bool *GetNeeds2ndDerivatives()
    { return Needs2ndDerivatives; };

    CoeffFct3D *GetCoeffFct()
    { return Coeffs; }
};

/******************************************************************************/
//
// Routines for initializing discrete forms
//
/******************************************************************************/

void InitializeDiscreteForms(  
  TDiscreteForm3D *&DiscreteFormGalerkin,
  TDiscreteForm3D *&DiscreteFormSDFEM,
  TDiscreteForm3D *&DiscreteFormUpwind,
  TDiscreteForm3D *&DiscreteFormSGS,
  TDiscreteForm3D *&DiscreteFormVMSProjection,
  TDiscreteForm3D *&DiscreteFormNLGalerkin,
  TDiscreteForm3D *&DiscreteFormNLSDFEM,
  TDiscreteForm3D *&DiscreteFormNLUpwind,
  TDiscreteForm3D *&DiscreteFormNLSGS,
  TDiscreteForm3D *&DiscreteFormNLVMSProjection,
  TDiscreteForm3D *&DiscreteFormNLSDFEM_DivDiv,
  TDiscreteForm3D *&DiscreteFormPressSep,
  TDiscreteForm3D *&DiscreteFormPressSepAuxProb,
  TDiscreteForm3D *&DiscreteFormNSRFBRhs,
  CoeffFct3D *LinCoeffs, int NSTYPE);

void InitializeDiscreteFormsFriction(  
  TDiscreteForm3D *&DiscreteFormGalerkin,
  TDiscreteForm3D *&DiscreteFormSDFEM,
  TDiscreteForm3D *&DiscreteFormUpwind,
  TDiscreteForm3D *&DiscreteFormSGS,
  TDiscreteForm3D *&DiscreteFormNLGalerkin,
  TDiscreteForm3D *&DiscreteFormNLSDFEM,
  TDiscreteForm3D *&DiscreteFormNLUpwind,
  TDiscreteForm3D *&DiscreteFormNLSGS,
  TDiscreteForm3D *&DiscreteFormGalerkinFriction,
  TDiscreteForm3D *&DiscreteFormNLGalerkinFriction,
  CoeffFct3D *LinCoeffs, int NSTYPE);

void InitializeDiscreteForms(  
  TDiscreteForm3D *&DiscreteFormGalerkin,
  TDiscreteForm3D *&DiscreteFormUpwind,
  TDiscreteForm3D *&DiscreteFormUpwindNC,
  TDiscreteForm3D *&DiscreteFormSmagorinsky,
  TDiscreteForm3D *&DiscreteFormClassicalLES,
  TDiscreteForm3D *&DiscreteFormGL00Convolution, 
  TDiscreteForm3D *&DiscreteFormGL00AuxProblem, 
  TDiscreteForm3D *&DiscreteFormVMS_Projection, 
  TDiscreteForm3D *&DiscreteFormVMS_SUPG, 
  TDiscreteForm3D *&DiscreteFormLerayAlpha,
  TDiscreteForm3D *&DiscreteFormNLGalerkin,
  TDiscreteForm3D *&DiscreteFormNLUpwind, 
  TDiscreteForm3D *&DiscreteFormNLUpwindNC, 
  TDiscreteForm3D *&DiscreteFormNLSmagorinsky,
  TDiscreteForm3D *&DiscreteFormNLClassicalLES, 
  TDiscreteForm3D *&DiscreteFormNLGL00Convolution,
  TDiscreteForm3D *&DiscreteFormNLGL00AuxProblem,
  TDiscreteForm3D *&DiscreteFormNLVMS_Projection, 
  TDiscreteForm3D *&DiscreteFormNLVMS_ProjectionExpl, 
  TDiscreteForm3D *&DiscreteFormNSRFBExplRhs,
  TDiscreteForm3D *&DiscreteFormNLVMS_SUPG,
  TDiscreteForm3D *&DiscreteFormNLLerayAlpha,
  TDiscreteForm3D *&DiscreteFormRHS,
  TDiscreteForm3D *&DiscreteFormRHSClassicalLES,
  TDiscreteForm3D *&DiscreteFormRHSLES,
  TDiscreteForm3D *&DiscreteFormRHSSUPG,
  TDiscreteForm3D *&DiscreteFormMatrixGL00AuxProblem,
  TDiscreteForm3D *&DiscreteFormGL00AuxProblemRHS,
  TDiscreteForm3D *&DiscreteFormMatrixAuxProblemU,
  TDiscreteForm3D *&DiscreteFormRHSAuxProblemU,
  TDiscreteForm3D *&DiscreteFormRHSNewton,
  TDiscreteForm3D *&DiscreteFormRHSNewtonNL,            
  TDiscreteForm3D *&DiscreteFormC, 
  TDiscreteForm3D *&DiscreteFormJ,
  CoeffFct3D *LinCoeffs, int NSTYPE);

void InitializeDiscreteFormsVMS(  
  TDiscreteForm3D *&DiscreteFormGalerkin,
  TDiscreteForm3D *&DiscreteFormUpwind,
  TDiscreteForm3D *&DiscreteFormSmagorinsky,
  TDiscreteForm3D *&DiscreteFormNLGalerkin,
  TDiscreteForm3D *&DiscreteFormNLUpwind, 
  TDiscreteForm3D *&DiscreteFormNLSmagorinsky,
  TDiscreteForm3D *&DiscreteFormRHS,
  TDiscreteForm3D *&DiscreteForm_ho_RHS,
  CoeffFct3D *LinCoeffs, int NSTYPE);
  
void InitializeDiscreteFormsOS_ST(
  TDiscreteForm3D *&DiscreteFormGalerkin,
  TDiscreteForm3D *&DiscreteFormStiffRhsOS_ST,
  CoeffFct3D *LinCoeffs, int NSTYPE);

void InitializeDiscreteForms
     ( TDiscreteForm3D *&DiscreteFormGalerkin,
       TDiscreteForm3D *&DiscreteFormNLGalerkin,
       CoeffFct3D *LinCoeffs, int NSTYPE);
       
void InitializeDiscreteForms
     ( TDiscreteForm3D *&DiscreteFormGalerkin,
       TDiscreteForm3D *&DiscreteFormNLGalerkin,
       TDiscreteForm3D *&DiscreteFormRHS,
       CoeffFct3D *LinCoeffs, int NSTYPE);
       
void InitializeDiscreteForms (TDiscreteForm3D *&DiscreteForm, CoeffFct3D *LinCoeff);

void InitializeDiscreteFormsScalar(TDiscreteForm3D *&DiscreteFormMRhs_Galerkin, TDiscreteForm3D *&DiscreteFormARhs_Galerkin, 
                                   TDiscreteForm3D *&DiscreteFormRhs, CoeffFct3D *LinCoeff);

void InitializeDiscreteFormGrid(TDiscreteForm3D *&DiscreteFormGrid,
                CoeffFct3D *GridCoeffs);
#endif