hdsim.hpp

Go to the documentation of this file.

#ifndef HDSIM_1D_HPP
#define HDSIM_1D_HPP 1

#include "../common/equation_of_state.hpp"
#include "spatial_distribution1d.hpp"
#include "../common/riemann_solver.hpp"
#include "vertex_motion.hpp"
#include "spatial_reconstruction1d.hpp"
#include "boundary_conditions_1d.hpp"
#include "source_term_1d.hpp"
#include "physical_geometry_1d.hpp"
#include "cell_updater_1d.hpp"
#include "simulation_state_1d.hpp"
#include "time_step_function_1d.hpp"
#include "../two_dimensional/extensive.hpp"
#include "extensive_updater_1d.hpp"
#include "flux_calculator_1d.hpp"

class hdsim1D
{
private:

  const PhysicalGeometry1D& pg_;

  SimulationState1D ss_;

  EquationOfState const& eos_;

  vector<Extensive> extensives_;

  VertexMotion const& vm_;

  SourceTerm1D const& force_;

  const TimeStepFunction1D& tsf_;

  const FluxCalculator1D& fc_;

  const ExtensiveUpdater1D& eu_;
  
  const CellUpdater1D& cu_;

  double time_;

  int cycle_;

  vector<vector<double> > tracers_intensive_;
  
  vector<vector<double> > tracers_extensive_;

public:

  hdsim1D
  (const PhysicalGeometry1D& pg,
   const SimulationState1D& ss,
   const EquationOfState& eos,
   const VertexMotion& vm,
   const SourceTerm1D& force,
   const TimeStepFunction1D& tsf,
   const FluxCalculator1D& fc,
   const ExtensiveUpdater1D& eu,
   const CellUpdater1D& cu);

  void TimeAdvance(void);

  void TimeAdvance2(void);

  const SimulationState1D& getState(void) const;

  // Diagnostics

  double GetTime(void) const;

  int GetCycle(void) const;

  //  vector<Conserved> const& getFluxes(void) const;
};

#endif // HDSIM_1D_HPP