Blockstructured Adaptive Mesh Refinement in object-oriented C++

---

       

---

Flagging.h

Go to the documentation of this file.

#ifndef AMROC_FLAGGING_H
#define AMROC_FLAGGING_H

#include "FlaggingInterface.h"
#include "Integrator.h"

#define GoodPoint (0)
#define BadPoint  (1)
#define VizServer 0                                               

#ifndef FlaggingName
#define Flagging(dim)      name2(Flagging,dim)
#define FlaggingName
#endif

#ifndef SolverName
#define Solver(dim)      name2(Solver,dim)
#define SolverName
#endif

template <class VectorType, class FixupType, class FlagType>
class Solver(DIM);

template <class VectorType, class FixupType, class FlagType>
class Flagging(DIM) : public FlaggingInterface(DIM)<VectorType,FlagType> {
  typedef Solver(DIM)<VectorType,FixupType,FlagType> solver_type;

  typedef GridFunction(DIM)<VectorType> vec_grid_fct_type;  
  typedef GridFunction(DIM)<typename VectorType::InternalDataType>   grid_fct_type;
  typedef GridFunction(DIM)<FlagType> flag_fct_type; 

public:
  typedef Integrator(DIM)<VectorType> integrator_type;

  Flagging(DIM)(solver_type& solver) : _Solver(solver), _EstimateError(false) {    
    _Hierarchy = (GridHierarchy*) 0;
    _u = (vec_grid_fct_type *) 0;
    _u_sh = (vec_grid_fct_type *) 0;
    _work = (grid_fct_type *) 0;
    _work_sh = (grid_fct_type *) 0;
    _flag = (flag_fct_type *) 0; 
  }

  ~Flagging(DIM)() {
    if (_flag)
      delete _flag;
  }

  virtual void register_at(ControlDevice& Ctrl) {}
  virtual void register_at(ControlDevice& Ctrl,const string& prefix) {
    LocCtrl = Ctrl.getSubDevice(prefix+"Flagging");
  }
  virtual void init() {}
  virtual void update() {}
  virtual void finish() {
    if (_flag) {
      delete _flag;
      _flag = (flag_fct_type *) 0; 
    }
  } 

  void SetBufferwidth(const int buffw) { _Bufferwidth = buffw; }
  inline const int& Bufferwidth() const { return _Bufferwidth; }
  inline int Bufferwidth() { return _Bufferwidth; }

  void SetupData() {
    int t_sten = 0;            
    int s_sten = Bufferwidth(); 
    sprintf(flagName,"FlagPlaceholder");
    _flag = new flag_fct_type(flagName, t_sten, s_sten, GH(), DAGHCommSimple);
    SetCheckpointFlag(*_flag, DAGHFalse);   
    SetProlongFlag(*_flag, DAGHFalse);   
    _flag->GF_SetMaxRecomposeLevel(DAGHNull);
  }    

  void FlagLevel(const int Time, const int Level, double t, double dt) {
    forall (Flags(),Time,Level,c)  
      Flags()(Time,Level,c) = GoodPoint;
    end_forall    
    SetFlags(Time, Level, t, dt);
  }   

  inline void SetGridHierarchy(GridHierarchy* gh) { _Hierarchy = gh; }
  inline GridHierarchy& GH() { return *_Hierarchy; }
  inline const GridHierarchy& GH() const { return *_Hierarchy; }

  inline void SetGridFunctions(vec_grid_fct_type* u, vec_grid_fct_type* ush, 
                               grid_fct_type* work, grid_fct_type* worksh) {
    _u = u; _u_sh = ush;
    _work = work; _work_sh = worksh; 
  }
  inline vec_grid_fct_type& U() { return *_u; }
  inline const vec_grid_fct_type& U() const { return *_u; }
  inline vec_grid_fct_type& Ush() { return *_u_sh; }
  inline const vec_grid_fct_type& Ush() const { return *_u_sh; }
  inline grid_fct_type& Work() { return *_work; }
  inline const grid_fct_type& Work() const { return *_work; }
  inline flag_fct_type& Flags() { return *_flag; }
  inline grid_fct_type& Worksh() { return *_work_sh; }
  inline const grid_fct_type& Worksh() const { return *_work_sh; }
  inline const flag_fct_type& Flags() const { return *_flag; }
  inline const bool& ErrorEstimation() const { return _EstimateError;}
  void SetErrorEstimation(const bool est) { _EstimateError = est; }

  inline solver_type& Solver_() { return _Solver; }
  inline const solver_type& Solver_() const { return _Solver; }
  inline integrator_type& Integrator_() { return Solver_().Integrator_(); }
  inline const integrator_type& Integrator_() const { return Solver_().Integrator_(); }

  inline const int& NEquations() const { return Integrator_().NEquations(); }
  inline const int& NGhosts() const { return Integrator_().NGhosts(); }

protected:
  solver_type& _Solver;
  GridHierarchy* _Hierarchy;
  vec_grid_fct_type *_u, *_u_sh;
  grid_fct_type *_work, *_work_sh;
  flag_fct_type* _flag;
  char flagName[DAGHBktGFNameWidth];
  int _Bufferwidth;
  bool _EstimateError;
  ControlDevice LocCtrl;
};


#endif

---

Quickstart     Users Guide     Programmers Reference     Installation      Examples     Download

---

AMROC Main      Home      Contact

last update: 06/01/04