Blockstructured Adaptive Mesh Refinement in object-oriented C++

---

       

---

Integrator.h

Go to the documentation of this file.

#ifndef AMROC_INTEGRATOR_H
#define AMROC_INTEGRATOR_H

#include "IntegratorInterface.h"
#include "TSSrcIntegrator.h"

#ifndef temporary_fluxes
#define temporary_fluxes(f, u) { \
  vec_grid_data_type* (f)[2*DIM]; \
  register int name2(f,_CU) = 0; \
  for (; name2(f,_CU)<2*DIM; name2(f,_CU)++) { \
    (f)[name2(f,_CU)] = new vec_grid_data_type((u).bbox()); }
#endif

#ifndef end_temporary_fluxes
#define end_temporary_fluxes(f) for (name2(f,_CU)=0; \
  name2(f,_CU)<2*DIM; name2(f,_CU)++) { \
    delete (f)[name2(f,_CU)]; } \
} 
#endif


#ifndef IntegratorName
#define Integrator(dim)      name2(Integrator,dim)
#define IntegratorName
#endif

template <class VectorType>
class Integrator(DIM) : public IntegratorInterface(DIM)<VectorType> {
  typedef GridData(DIM)<VectorType> vec_grid_data_type;
  typedef GridFunction(DIM)<VectorType> vec_grid_fct_type;  
  typedef TSSrcIntegrator(DIM)<VectorType> src_integrator_type;

public:
  Integrator(DIM)(const int gh) : 
    _Equations(VectorType::Length()), _Ghosts(gh), _SrcTreatment(0) {    
    _Hierarchy = (GridHierarchy*) 0;
    _SrcIntegrator = (src_integrator_type *) 0;
    _abort = DAGHFalse;
  }

  virtual void register_at(ControlDevice& Ctrl,const string& prefix) {
    ControlDevice LocCtrl = Ctrl.getSubDevice(prefix+"Integrator");
    if (_SrcIntegrator)
      SrcIntegrator().register_at(LocCtrl,prefix);
  }
  virtual void register_at(ControlDevice& Ctrl) {
    register_at(Ctrl, "");
  }

  virtual void init() {
    if (_SrcIntegrator)
      SrcIntegrator().init();
  }
  virtual void update() {
    if (_SrcIntegrator)
      SrcIntegrator().update();
  }
  virtual void finish() {
    if (_SrcIntegrator)
      SrcIntegrator().finish();
  } 

  virtual void SetupData() {
    if (_SrcIntegrator)
      SrcIntegrator().SetupData();
  }
   
  inline void SetGridHierarchy(GridHierarchy* gh) { 
    _Hierarchy = gh; 
    if (_SrcIntegrator)
      SrcIntegrator().SetGridHierarchy(gh);
  }
  inline GridHierarchy& GH() { return *_Hierarchy; }
  inline const GridHierarchy& GH() const { return *_Hierarchy; }

  inline void SetSrcIntegrator(src_integrator_type* src_integ) 
    { _SrcIntegrator = src_integ; }
  inline src_integrator_type& SrcIntegrator() { return *_SrcIntegrator; }
  inline const src_integrator_type& SrcIntegrator() const { return *_SrcIntegrator; }
  inline void SetSrcTreatment(int st) { _SrcTreatment = st; }
  inline const int& SrcTreatment() { return _SrcTreatment; }
  inline void SetMaxIntegratorPasses(int mp) { _MaxPass = mp; }
  inline const int& MaxIntegratorPasses() { return _MaxPass; }

  inline const int& Abort() const { return _abort; }
  void CheckGrid(vec_grid_data_type& StateVec, const BBox& where, const char *text) {
    if (!ControlGrid(StateVec, where) && _abort) {
      char Proc[32];
#ifndef DAGH_NO_MPI
      sprintf(Proc,"On Processor %d: ",comm_service::proc_me());
#else
      Proc[0] = '\0';
#endif 
      string ErrorText = "\n";
      ErrorText += Proc; ErrorText += " CheckGrid() at ";
      ErrorText += text; ErrorText += " not successful for ";
      cerr << ErrorText << StateVec.bbox() << "\nAborting programm...\n" << flush;
      cout << ErrorText << StateVec.bbox() << "\nAborting programm...\n" << flush;
#ifdef DEBUG_PRINT
      ( comm_service::log() << ErrorText << StateVec.bbox() 
        << "\nAborting programm...\n\n" ).flush();
      ( comm_service::log() << StateVec << "\n" ).flush();     
#endif
      assert(0);
    } 
  }

  inline const int& NEquations() const { return _Equations; }
  inline const int& NGhosts() const { return _Ghosts; }

protected:
  int _Equations;
  int _Ghosts;
  int _SrcTreatment;
  int _MaxPass;
  int _abort;
  GridHierarchy* _Hierarchy;
  vec_grid_fct_type* _u;
  src_integrator_type* _SrcIntegrator;
};


#endif

---

Quickstart     Users Guide     Programmers Reference     Installation      Examples     Download

---

AMROC Main      Home      Contact

last update: 06/01/04