Department of Scientific Computing    Institute for Numerical Simulation    University of Bonn

Documentation

Download

Programming References

Bug Reports / Suggestions

FAQ

Authors

           

---

LevelAdaptive.hpp

//
//  Main classes of non adaptive sparse grid code
//

#ifndef _LEVELADAPTIVE_
# define _LEVELADAPTIVE_

#include "UniformData.hpp"


template<int DIM>
struct LevelAdaptiveData;

template<int DIM>
class BlendingData ;

template<int DIM>
struct BlendingWeights;


template<int DIM>
class LevelAdaptiveGrid {
public:
      LevelAdaptiveGrid() ;
      LevelAdaptiveGrid(int *L0, int *L, Wavelets *W) ;
     ~LevelAdaptiveGrid() ;

  void Init (int *L0, int *L, Wavelets *W)  ;
  void Clear() ;
  void Set  (Function *F , double c) ;
  void Print() ;
  size_t Size() ;
  
  void Attach(LevelAdaptiveData<DIM> *D) ;  
  void Attach(BlendingWeights<DIM>   *B) ;

  Matrix<bool,DIM> *a    ;
  int Level0[DIM],L[DIM] ;
  Wavelets         *W    ;
  // list of LevelAdaptiveData attached to *this
  int nLAD,nBW ;
  LevelAdaptiveData<DIM> *LAD[2000] ;
  BlendingWeights<DIM>   *BW [2000] ;
} ;



template<int DIM>
class LevelAdaptiveData {
public:
        LevelAdaptiveData() ;
        LevelAdaptiveData(LevelAdaptiveGrid<DIM> *AG) ;
       ~LevelAdaptiveData() ;

  void   Init    (LevelAdaptiveGrid<DIM> *AG) ;
  void   Allocate() ;
  void   Clear() ;
  size_t Size() ;

  void   SetBoundaryConditions(int BCs[DIM][2])     ;
  void   SetBoundaryConditions(LevelAdaptiveData<DIM> *S)   ; // Copy from (*S)
  void   CheckBC(int BC0[DIM][2] , int BC1[DIM][2]) ;
  
  void   Compatible(LevelAdaptiveData<DIM> *S) ;
  void   Compatible(LevelAdaptiveData<DIM> *S0,LevelAdaptiveData<DIM> *S1) ;
  void   Compatible(LevelAdaptiveData<DIM> *S0,LevelAdaptiveData<DIM> *S1,LevelAdaptiveData<DIM> *S2) ;
  void   Compatible(LevelAdaptiveData<DIM> *S0,LevelAdaptiveData<DIM> *S1,LevelAdaptiveData<DIM> *S2,LevelAdaptiveData<DIM> *S3) ;
  void   Compatible(LevelAdaptiveData<DIM> *S0,LevelAdaptiveData<DIM> *S1,LevelAdaptiveData<DIM> *S2,LevelAdaptiveData<DIM> *S3,LevelAdaptiveData<DIM> *S4) ;
  
  void   CopyExtensions(LevelAdaptiveData<DIM> *S) ;

  void   Blend (BlendingData<DIM> *B) ;
  
  void   PPlus(const double c0,LevelAdaptiveData<DIM> *S0) ;
  void   Add(const double c0,LevelAdaptiveData<DIM> *S0 , const double c1,LevelAdaptiveData<DIM> *S1) ;
  void   Add(const double c0,LevelAdaptiveData<DIM> *S0 , const double c1,LevelAdaptiveData<DIM> *S1 , const double c2,LevelAdaptiveData<DIM> *S2) ;
  void   Add(const double c0,LevelAdaptiveData<DIM> *S0 , const double c1,LevelAdaptiveData<DIM> *S1 , const double c2,LevelAdaptiveData<DIM> *S2 , 
             const double c3,LevelAdaptiveData<DIM> *S3) ;
  void   Sub(LevelAdaptiveData<DIM> *S0 , LevelAdaptiveData<DIM> *S1) ;
  void   Mul(double f) ;
  double InnerProd(LevelAdaptiveData<DIM> *S) ;
  double Max() ;
  double Min() ;
  double MaxAbs() ;
  
  void   Copy(LevelAdaptiveData<DIM> *From) ;
  void   Set(const double c) ;

  void   ToUniform  (UniformData<DIM>   *M) ;
  void   ToUniform  (Matrix<double,DIM> *M) ;
  void   FromUniform(UniformData<DIM>   *M) ;
  void   FromUniform(Matrix<double,DIM> *M) ;

  void   Multiply(LevelAdaptiveData<2*DIM> *K , LevelAdaptiveData<DIM> *S) ; 
  void   ApplyOp (int *BCT , LevelAdaptiveData<DIM> *X , int dir , unsigned int op) ;
  void   ApplyOp (            LevelAdaptiveData<DIM> *X , int dir , unsigned int op) ;
  
  void   WriteUDF(const char *name, UniformData<DIM> *M, LevelAdaptiveData<DIM> *Tmp=NULL, bool CastToFloat=false) ;
  void   WriteUDF(const char *name,             int  *L, LevelAdaptiveData<DIM> *Tmp=NULL, bool CastToFloat=false) ;
  void   WriteUDF(const char *name,             int   L, LevelAdaptiveData<DIM> *Tmp=NULL, bool CastToFloat=false) ;
  
  Matrix< Matrix<double,DIM>* , DIM> *a ;
  LevelAdaptiveGrid<DIM> *G             ;
  Extensions<DIM>       Ext             ;
} ;


template<int DIM>
class BlendingWeights{
public:
      BlendingWeights() ;
      BlendingWeights(LevelAdaptiveGrid<DIM> *G) ;
     ~BlendingWeights() ;

  void Init (LevelAdaptiveGrid<DIM> *G) ;
  void Clear() ;
  void Print(int all = 3) ;
  size_t Size() ;
  //  compute blending wheights
  void Set  (LevelAdaptiveGrid<DIM> *A) ;
  //  attach an instance of BlendingData to *this
  void Attach(BlendingData<DIM> *B) ;

  Matrix<int , DIM> *a ;
  LevelAdaptiveGrid<DIM> *G ;
  
  int nBD ;
  BlendingData<DIM> *BD[2000] ;
} ;


template<int DIM>
class BlendingData {
public: 
        BlendingData() ;
        BlendingData(BlendingWeights<DIM> *A) ;
       ~BlendingData() ;

  void   Init (BlendingWeights<DIM> *A) ;
  void   Clear() ;
  void   Print() ;
  void   Allocate() ;

  void   SetBoundaryConditions(int BCs[DIM][2]) ;
  void   CheckBC(int BC0[DIM][2] , int BC1[DIM][2]) ;

  void   Compatible(BlendingData<DIM> *S) ;
  void   Compatible(BlendingData<DIM> *S0, BlendingData<DIM> *S1) ;
  void   Compatible(BlendingData<DIM> *S0, BlendingData<DIM> *S1, BlendingData<DIM> *S2) ;
  void   Compatible(BlendingData<DIM> *S0, BlendingData<DIM> *S1, BlendingData<DIM> *S2, BlendingData<DIM> *S3) ;
  
  void   CopyExtensions(BlendingData<DIM> *S) ;
  
  void   Get(LevelAdaptiveData<DIM> *S) ;
  void   SetFunction(Function *F) ;
  void   Set(const double c) ;
  
  void   Add(const double c0,BlendingData<DIM> *S0 , const double c1,BlendingData<DIM> *S1) ;
  void   Add(const double c0,BlendingData<DIM> *S0 , const double c1,BlendingData<DIM> *S1 , 
             const double c2,BlendingData<DIM> *S2) ;
  void   Add(const double c0,BlendingData<DIM> *S0 , const double c1,BlendingData<DIM> *S1 , 
             const double c2,BlendingData<DIM> *S2 , const double c3,BlendingData<DIM> *S3) ;

  void   Sub(BlendingData<DIM> *S0 , BlendingData<DIM> *S1) ;
  double Max() ;
  double Min() ;
  double MaxAbs() ;
  
  void   ApplyOp(int *BCT , BlendingData<DIM> *S , int dir , unsigned int op) ;
  
  Matrix<Matrix<double,DIM> * , DIM> *a ;
  BlendingWeights<DIM> *A ;
  Extensions<DIM>       Ext ;
} ;


# include "LevelAdaptiveGrid.hpp"
# include "BlendingWeights.hpp"
# include "LevelAdaptiveData.hpp"
# include "BlendingData.hpp"
#endif

---