CD_TimeStepper.H

Go to the documentation of this file.

/* chombo-discharge
 * Copyright © 2021 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
#ifndef CD_TimeStepper_H
#define CD_TimeStepper_H
 
// Chombo includes
#include <CH_HDF5.H>
 
// Our includes
#include <CD_ComputationalGeometry.H>
#include <CD_MultiFluidIndexSpace.H>
#include <CD_AmrMesh.H>
#include <CD_NamespaceHeader.H>
 
class TimeStepper
{
public:
  TimeStepper();
 
  virtual ~TimeStepper();
 
  void
  setAmr(const RefCountedPtr<AmrMesh>& a_amr);
 
  void
  setComputationalGeometry(const RefCountedPtr<ComputationalGeometry>& a_computationalGeometry);
 
  virtual void
  setupSolvers() = 0;
 
  virtual void
  allocate() = 0;
 
  virtual void
  initialData() = 0;
 
  virtual void
  postInitialize() = 0;
 
  virtual void
  postCheckpointSetup() = 0;
 
  virtual void
  registerRealms() = 0;
 
  virtual void
  registerOperators() = 0;
 
  virtual void
  parseRuntimeOptions();
 
#ifdef CH_USE_HDF5
  virtual void
  readCheckpointHeader(HDF5HeaderData& a_header);
#endif
 
#ifdef CH_USE_HDF5
  virtual void
  writeCheckpointHeader(HDF5HeaderData& a_header) const;
#endif
 
#ifdef CH_USE_HDF5
  virtual void
  writeCheckpointData(HDF5Handle& a_handle, const int a_lvl) const = 0;
#endif
 
#ifdef CH_USE_HDF5
  virtual void
  readCheckpointData(HDF5Handle& a_handle, const int a_lvl) = 0;
#endif
 
  virtual int
  getNumberOfPlotVariables() const = 0;
 
  virtual Vector<std::string>
  getPlotVariableNames() const = 0;
 
  virtual void
  writePlotData(LevelData<EBCellFAB>& a_output,
                int&                  a_icomp,
                const std::string     a_outputRealm,
                const int             a_level) const = 0;
 
  virtual void
  prePlot();
 
  virtual void
  postPlot();
 
  virtual Vector<long int>
  getCheckpointLoads(const std::string a_realm, const int a_level) const;
 
  virtual Real
  computeDt() = 0;
 
  virtual Real
  advance(const Real a_dt) = 0;
 
  virtual void
  synchronizeSolverTimes(const int a_step, const Real a_time, const Real a_dt) = 0;
 
  virtual void
  printStepReport() = 0;
 
  virtual void
  preRegrid(const int a_lmin, const int a_oldFinestLevel) = 0;
 
  virtual void
  regrid(const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel) = 0;
 
  virtual void
  postRegrid() = 0;
 
  virtual bool
  needToRegrid();
 
  virtual bool
  loadBalanceThisRealm(const std::string a_realm) const;
 
  virtual void
  loadBalanceBoxes(Vector<Vector<int>>&             a_procs,
                   Vector<Vector<Box>>&             a_boxes,
                   const std::string                a_realm,
                   const Vector<DisjointBoxLayout>& a_grids,
                   const int                        a_lmin,
                   const int                        a_finestLevel);
 
protected:
  int m_verbosity;
 
  int m_timeStep;
 
  Real m_time;
 
  Real m_dt;
 
  RefCountedPtr<AmrMesh> m_amr;
 
  RefCountedPtr<ComputationalGeometry> m_computationalGeometry;
};
 
#include <CD_NamespaceFooter.H>
 
#endif