amroc/amr/Integrator.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2002 Ralf Deiterding
// Brandenburgische Universitaet Cottbus
//
// Copyright (C) 2003-2007 California Institute of Technology
// Ralf Deiterding, ralf@cacr.caltech.edu

#ifndef AMROC_INTEGRATOR_H
#define AMROC_INTEGRATOR_H

#include <iostream>
#include <string>
#include <cassert>
#include <cstdio>
#include "AMRBase.h"

template <class VectorType, int dim>
class Integrator : public AMRBase<VectorType,dim> {
  typedef AMRBase<VectorType,dim> base;
public:
  typedef typename base::vec_grid_fct_type vec_grid_fct_type;  
  typedef typename base::vec_grid_data_type vec_grid_data_type;

  Integrator() : base(), _MaxPass(0), _abort(DAGHFalse) {}

  virtual ~Integrator() {}

  //******************************************************************************
  // Abstract class interface
  //******************************************************************************
  virtual double CalculateGrid(vec_grid_data_type& NewStateVec, 
                               vec_grid_data_type& OldStateVec,
                               vec_grid_data_type* Flux[],
                               const int& level, 
                               const double& t, const double& dt, 
                               const int& mpass) = 0;
  virtual void AllocGridFluxes(const BBox &bb, vec_grid_data_type**& Flux)  = 0;
  virtual void DeAllocGridFluxes(vec_grid_data_type**& Flux) = 0;
  virtual void ResetGridFluxes(vec_grid_data_type**& Flux)  = 0;
  virtual int ControlGrid(vec_grid_data_type& StateVec, const int& level, 
                          const BBox& where, const double& time, const int verbose) = 0;  
  virtual int NMethodOrder() const = 0;
  //******************************************************************************

  virtual void register_at(ControlDevice& Ctrl, const std::string& prefix) {
    base::LocCtrl = Ctrl.getSubDevice(prefix+"Integrator");
  }
  virtual void register_at(ControlDevice& Ctrl) {
    register_at(Ctrl, "");
  }

  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 int& level, 
                 const BBox& where, const double& time, const char *text) {
    if (ControlGrid(StateVec, level, where, time, 1)==0 && _abort) {
      char Proc[32];
#ifndef DAGH_NO_MPI
      std::sprintf(Proc,"On Processor %d: ",comm_service::proc_me());
#else
      Proc[0] = '\0';
#endif 
      std::string ErrorText = "\n";
      ErrorText += Proc; ErrorText += " CheckGrid() at ";
      ErrorText += text; ErrorText += " not successful for ";
      std::cerr << ErrorText << StateVec.bbox()
                << "\nAborting program...\n" << std::flush;
      std::cout << ErrorText << StateVec.bbox()
                << "\nAborting program...\n" << std::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 void SetGridFunction(vec_grid_fct_type* u) { _u = u; }
  inline vec_grid_fct_type& U() { return *_u; }
  inline const vec_grid_fct_type& U() const { return *_u; }

protected:
  int _MaxPass;
  int _abort;
  vec_grid_fct_type *_u;
};


#endif

---