SSInfall.h

Go to the documentation of this file.

#pragma once
 
#include <iomanip>
 
#include "AbstractSSInfall.h"
#include "tarch/logging/Log.h"
 
namespace applications::exahype2::euler::sphericalaccretion {
  class SSInfall;
} // namespace applications::exahype2::euler::sphericalaccretion
 
 
class applications::exahype2::euler::sphericalaccretion::SSInfall: public AbstractSSInfall {
private:
  static tarch::logging::Log _log;
 
  std::string plotGlobalMTot(bool plotCopy) const;
  std::string plotRs() const;
 
#if defined(GPUOffloadingOMP)
#pragma omp declare target
#endif
  static double _rMax;
 
  static double _rhoMax;
 
  static double _mTotMax;
 
  // @todo I've just removed this field. Couldn't find out where it is used.
  //    double t_record=0.0;
 
#if defined(GPUOffloadingOMP)
#pragma omp end declare target
#endif
 
  double getInitialMIn(double r_coord) const;
 
  static double getForceDensityNorm(double r_coord, double m_in, double t, double density, double a_input);
 
public:
  static const tarch::la::Vector<Dimensions, double> OverDensityCentre;
 
  void startTimeStep(
    double globalMinTimeStamp, double globalMaxTimeStamp, double globalMinTimeStepSize, double globalMaxTimeStepSize
  ) override;
 
  void finishTimeStep() override;
 
  ::exahype2::RefinementCommand refinementCriterion(
    const double* __restrict__ Q, // Q[5+0],
    const tarch::la::Vector<Dimensions, double>& volumeCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t
  ) override;
 
 
  void initialCondition(
    double* __restrict__ Q,
    const tarch::la::Vector<Dimensions, double>& volumeCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    bool                                         gridIsConstructed
  ) override;
 
 
  void sourceTerm(
    const double* __restrict__ Q,
    const tarch::la::Vector<Dimensions, double>& volumeCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    double                                       dt,
    double* __restrict__ S
  ) override;
 
 
#if defined(GPUOffloadingOMP)
#pragma omp declare target
#endif
  static void sourceTerm(
    const double* __restrict__ Q,
    const tarch::la::Vector<Dimensions, double>& volumeCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    double                                       dt,
    double* __restrict__ S,
    Offloadable
  );
#if defined(GPUOffloadingOMP)
#pragma omp end declare target
#endif
 
 
  virtual void boundaryConditions(
    const double* __restrict__ Qinside, // Qinside[5+0]
    double* __restrict__ Qoutside,      // Qoutside[5+0]
    const tarch::la::Vector<Dimensions, double>& faceCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    int                                          normal
  ) override;
 
 
  virtual double maxEigenvalue(
    const double* __restrict__ Q, // Q[5+0],
    const tarch::la::Vector<Dimensions, double>& faceCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    double                                       dt,
    int                                          normal
  ) override;
 
 
#if defined(GPUOffloadingOMP)
#pragma omp declare target
#endif
  static double maxEigenvalue(
    const double* __restrict__ Q, // Q[5+0],
    const tarch::la::Vector<Dimensions, double>& faceCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    double                                       dt,
    int                                          normal,
    Offloadable
  );
#if defined(GPUOffloadingOMP)
#pragma omp end declare target
#endif
 
  virtual void flux(
    const double* __restrict__ Q, // Q[5+0],
    const tarch::la::Vector<Dimensions, double>& faceCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    double                                       dt,
    int                                          normal,
    double* __restrict__ F // F[5]
  ) override;
 
 
#if defined(GPUOffloadingOMP)
#pragma omp declare target
#endif
  static void flux(
    const double* __restrict__ Q, // Q[5+0],
    const tarch::la::Vector<Dimensions, double>& faceCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    double                                       dt,
    int                                          normal,
    double* __restrict__ F, // F[5]
    Offloadable
  );
#if defined(GPUOffloadingOMP)
#pragma omp end declare target
#endif
 
 
  void add_mass(const double r_coor, const double rho, const double size);
 
  double mass_interpolate(const double r_coor, const int MassCal);
 
  virtual bool patchCanUseStatelessPDETerms(
    const tarch::la::Vector<Dimensions, double>& patchCentre,
    const tarch::la::Vector<Dimensions, double>& patchH,
    double                                       t,
    double                                       dt
  ) const;
 
 
  bool isOutsideOfLargestRadius(
    const tarch::la::Vector<Dimensions, double>& patchCentre, const tarch::la::Vector<Dimensions, double>& patchH
  ) const;
};