ADISystem1D.h

// =======================================================================
// @(#)ADISystem1D.h        4.1 13.11.09
// 
// Class:       TADISystem1D
// Purpose:     class for  ADISystem1D

//
// Author:      Sashikumaar Ganesan (13.11.09)
//
// History:     start of implementation 13.11.09 (Sashikumaar Ganesan)
//
// =======================================================================

#ifndef __ADISYSTEM1D__
#define __ADISYSTEM1D__

#include <ADISystem.h>
#include <FESpace1D.h>
#include <FEFunction1D.h>

class TADISystem1D : public TADISystem
{
  protected:

    double X;

    double Y;

    double Z;   

    TFESpace1D *FESpace1D_Intl;

    TFEFunction1D *FEFunction_Intl;

    TSquareMatrix1D *M_Intl;

    TSquareMatrix1D *A_Intl;

    TSquareMatrix1D *S_Intl;

    TSquareMatrix1D *K_Intl;


     DoubleFunctND *Initial;

     DoubleFunct2D  *BDValue_LMin;
     DoubleFunct2D  *BDVal_LMax;

  private:
   int ConstructAllInfo();

  public:
    TADISystem1D(TFEFunction1D *FEFunction_intl, TSquareMatrix1D *M_internal, TSquareMatrix1D *A_internal,
                           TSquareMatrix1D *S_internal, TSquareMatrix1D *K_internal,
                           double x, double y, DoubleFunctND *initial, 
                           DoubleFunct2D *bDValue_LMin, DoubleFunct2D *bDVal_LMax,
                           double *Sol,  double *OldSol, double *b, 
                           double *Defect, double *Intlposl, DoubleFunctND *growthAndNuc);

    TADISystem1D(TFEFunction1D *FEFunction_intl, TSquareMatrix1D *M_internal, TSquareMatrix1D *A_internal,
                 double x, double y, double z, double *Sol, double *OldSol, double *b, double *Defect,
                 DoubleFunctND *growthfunct);
  
    //void Interpolate(double *Sol, DoubleFunct3D *Exact);
    void SetDirichletBc(BoundCond cond_Lmin, BoundCond cond_Lmax, double BDValue_Lmin, double BDValue_Lmax);

    void Interpolate(double *Sol, DoubleFunctND *Exact);
    
    void SolveAllLevels(double *MatValues_orig, double *SolFromX, double *AggrRhs, double &C, double C_Sat, 
                        double &T, CoeffFctND *BilinearCoeffs, 
                        double tau, BoundCond cond_Lmin, BoundCond cond_Lmax,  DoubleFunct3D *Exact );

    void AssembleInitRhs(double C, double T, CoeffFctND *BilinearCoeffs, 
                         BoundCond cond_Lmin, BoundCond cond_Lmax);
 
    void AssembleARhs_FD(double Conv, CoeffFctND *Bilinear);
    
    void AssembleARhs(double Conv, CoeffFctND *Bilinear);
    
    void AssembleARhs_SUPG(double Conv, CoeffFctND *Bilinear);

    void AssembleARhs_DG(double Conv, double Bnuc, CoeffFctND *Bilinear, BoundCond cond_Lmin, BoundCond cond_Lmax);
    
    double GetWeightedF();

   int GetN_InternalLevels()
    {
     return N_V;
    }

    double GetQ3Max(double *currsol);

    ~TADISystem1D();

};

#endif