Blockstructured Adaptive Mesh Refinement in object-oriented C++

---

       

---

TSSrcIntegrator.h

Go to the documentation of this file.

#ifndef AMROC_TSSRCINTEGRATOR_H
#define AMROC_TSSRCINTEGRATOR_H

#include "ODEIntegrators/ODEIntegrator.h"
#include "TSSrcIntegratorInterface.h"

#ifndef TSSrcIntegratorName
#define TSSrcIntegrator(dim)      name2(TSSrcIntegrator,dim)
#define TSSrcIntegratorName
#endif

template <class VectorType>
class TSSrcIntegrator(DIM) : public TSSrcIntegratorInterface(DIM)<VectorType> {
  typedef ODEIntegrator(DIM)<VectorType> ode_integrator_type;
  typedef GridDataFunction(DIM)<VectorType> grid_data_fct_type;
  typedef GridData(DIM)<VectorType> vec_grid_data_type;

public:
  TSSrcIntegrator(DIM)(const int gh) : 
    _Equations(VectorType::Length()), _Ghosts(gh) {
    _Hierarchy = (GridHierarchy*) 0;
    _TimeIntegrator = (ode_integrator_type *)0;
    _SrcTerm = (grid_data_fct_type *)0;
  }

  virtual void register_at(ControlDevice& Ctrl,const string& prefix) {
    ControlDevice LocCtrl = Ctrl.getSubDevice(prefix+"TSSrcIntegrator");
   if (_TimeIntegrator)
      TimeIntegrator().register_at(LocCtrl,prefix);
   if (_TimeIntegrator)
      SrcTerm().register_at(LocCtrl,prefix);
  }
  virtual void register_at(ControlDevice& Ctrl) {
    register_at(Ctrl, "");
  }

  virtual void init() {
    if (_TimeIntegrator)
      TimeIntegrator().init();
    if (_SrcTerm)
      SrcTerm().init();
  }

  virtual void update() {
    if (_TimeIntegrator)
      TimeIntegrator().update();
    if (_SrcTerm)
      SrcTerm().update();
  }

  virtual void finish() {
    if (_TimeIntegrator)
      TimeIntegrator().finish();
    if (_SrcTerm)
      SrcTerm().finish();  
  } 

  virtual void SetupData() {
    if (_TimeIntegrator)
      TimeIntegrator().SetupData();
    if (_SrcTerm)
      SrcTerm().SetupData();     
  }

  inline void SetSrcTerm(grid_data_fct_type* src) { 
    _SrcTerm = src;
    if (_TimeIntegrator)
      TimeIntegrator().SetGridDataFct(src); 
  }
  inline grid_data_fct_type& SrcTerm() { return *_SrcTerm;  }
  inline const grid_data_fct_type& SrcTerm() const { return *_SrcTerm; }

  inline void SetTimeIntegrator(ode_integrator_type* integ) { _TimeIntegrator = integ; }
  inline ode_integrator_type& TimeIntegrator() { return *_TimeIntegrator; }
  inline const ode_integrator_type& TimeIntegrator() const { return *_TimeIntegrator; }

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

  /* Before the transport step, the source term has to be updated 
     even on the ghost cells! Use bnd==1 for this.
     See CLAWPACK-Manual Note #5 "Source Terms".
     Otherwise the boundary conditions are not consistent with the
     solution anymore. */

  virtual void CalculateSrcOnGrid(vec_grid_data_type& StateVec, double t, 
                                  double dt, DCoords& dx, const int& bnd) {
    if (_TimeIntegrator && _SrcTerm) {
      BBox OpBox = StateVec.bbox();
      if (bnd == 0) 
        for (int d=0; d<DIM; d++) 
          if (OpBox.extents(d) > 2*NGhosts()) {
            OpBox.growlower(d,-NGhosts()); 
            OpBox.growupper(d,-NGhosts()); 
          }
      TimeIntegrator().IntegrateGrid(StateVec, OpBox, t, dt, dx);
    }
  }

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

protected:
  int _Equations;
  int _Ghosts;
  GridHierarchy* _Hierarchy;
  ode_integrator_type* _TimeIntegrator;
  grid_data_fct_type* _SrcTerm;
};


#endif

---

Quickstart     Users Guide     Programmers Reference     Installation      Examples     Download

---

AMROC Main      Home      Contact

last update: 06/01/04