Department of Scientific Computing    Institute for Numerical Simulation    University of Bonn

Documentation

Download

Programming References

Bug Reports / Suggestions

FAQ

Authors

           

---

AdaptiveLevels.hpp

//
//  Univariate classes 
//

#ifndef _UNI_VARIAT_
# define _UNI_VARIAT_

# include<map>
# include<pthread.h>

#include "Adaptive1D.hpp"

struct AdaptiveLevels {

        AdaptiveLevels() {;}
       ~AdaptiveLevels() {;}
 
         void Init() ;
         void Clear() ;
         void Delete() ;
         size_t     size() ;

  typedef map <int,long,intlt> AdaptiveL ;
  inline AdaptiveL *operator[](const unsigned int t)  { return a[t] ;}  
         AdaptiveL *a[LMAX+1] ;
} ;

struct AdaptiveLevelsIndex {
  int           dim  ;
  unsigned char l[DMAX-1] ; // DIM-1
  unsigned int  t[DMAX-1] ; // DIM-1
  
  AdaptiveLevelsIndex() {dim =-1;}
  void Init () {dim=-1;}
  void Clear() {dim=-1;}
  void Delete() {dim=-1;}

  void Print() ;
} ;

struct AdaptiveLevelsIndexCmp {
  bool operator()(const AdaptiveLevelsIndex *a,const AdaptiveLevelsIndex *b) const ;
} ;

struct AdaptiveLevelsPointer {
  size_t               s  ; // size of AdaptiveLevels *a
  AdaptiveLevels      *a  ;
  AdaptiveLevelsIndex *lt ; // dimension reduced level/spatial index 
                  AdaptiveLevelsPointer() ;
  void            Copy(AdaptiveLevelsPointer *alp) ;
} ;

int   AdaptiveLevelsPointercmp(const void *pa,const void *pb) ;

struct SMPjob {
          SMPjob() ;
  size_t  leveloffset(int *l) ;
  void    set(int operation, int direction, int dim,
              AdaptiveLevelsPointer *alpstart, size_t num_alp,
              vector<double> **input,
              vector<double> **output,
              vector<double> **roespeed=NULL) ;

  bool    isfinished() ; // called by main program to check whether the job has been processed
  void    setpoll() ;
  bool    getpoll() ;
  void    clrpoll() ;

  struct  AdaptiveLevelsPointer *alp ; // first AdaptiveLevel
  size_t  nalp                       ; // number of AdaptiveLevel

  AdaptiveB  *AB1,*AB2 ;
  AdaptiveGS *AGS,*AG2 ;

  int       op, dir, DIM ;
  Wavelets *WC ;

  vector<double> **in, **out, **wenoroe ;

  bool finished, terminate, poll ;

  size_t           jobcount,requests ;

} ;

size_t IndexSetWriteToAdaptiveL(AdaptiveLevels::AdaptiveL *al,size_t *counter,struct IndexSet *I) ;
void   ApplyUnivariateOperators(AdaptiveB *AB1, AdaptiveB *AB2, AdaptiveB **ABX, AdaptiveGS *AGS, AdaptiveGS *AG2, 
                                int op, Wavelets *WC) ;

void *ProcessSMPjob    (void *jb) ;
void *ProcessSMPjobWENO(void *jb) ;

void *PollProcessSMPjob(void *jb) ;

#endif

---