SelfSimilarInfallFD4.h

Go to the documentation of this file.

//
// ExaHyPE2 solver file
// Generated by Peano's Python API
// www.peano-framework.org
//
// This is generated. If you change fundamental properties, you will have to 
// generate this file. Backup your manual changes before you do so.
//
#pragma once
 
#include "AbstractSelfSimilarInfallFD4.h"
#include "tarch/logging/Log.h"
 
#include "spherical-accretion/MassAccumulator.h"
#include "spherical-accretion/GravityModel.h"
 
#include "EulerKernels.h"
 
namespace benchmarks {
  namespace exahype2 {
    namespace euler {
      namespace sphericalaccretionupscaling {
        class SelfSimilarInfallFD4;
      }
    }
  }
}
 
 
class benchmarks::exahype2::euler::sphericalaccretionupscaling::SelfSimilarInfallFD4: public benchmarks::exahype2::euler::sphericalaccretionupscaling::AbstractSelfSimilarInfallFD4 {
  private:
    static tarch::logging::Log   _log;
 
    ::applications::exahype2::euler::sphericalaccretion::MassAccumulator  _accumulator;
 
    static double TotalMassInPreviousTimeStep;
 
  public:
    static constexpr double BaseDensity         = 0.1;
    static constexpr double Gamma               = 5.0/3.0;
    static constexpr double pInitial            = 1e-6;
    static constexpr double aInitial            = 0.001;
    static constexpr double AdditionalMass      = 0.15;
    static constexpr double InitialTopHatRadius = 0.05;
 
    SelfSimilarInfallFD4();
 
 
    void addDensity(
      const tarch::la::Vector<Dimensions,double>&  volumeCentre,
      const tarch::la::Vector<Dimensions,double>&  volumeH,
      double                                       density
    );
 
    
    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;
    
 
    
    
    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;
    
 
    
  public:
    
    void sourceTerm(
      const double * __restrict__ Q,
      const tarch::la::Vector<Dimensions,double>&  gridCellCentre,
      const tarch::la::Vector<Dimensions,double>&  gridCellH,
      double                                       t,
      double                                       dt,
      double * __restrict__ S
    ) override;
    
 
    
    virtual double maxEigenvalue(
      const double * __restrict__ Q, // Q[5+0],
      const tarch::la::Vector<Dimensions,double>&  faceCentre,
      const tarch::la::Vector<Dimensions,double>&  gridCellH,
      double                                       t,
      double                                       dt,
      int                                          normal
    ) override;
    
 
 
    
    virtual void flux(
      const double * __restrict__ Q, // Q[5+0],
      const tarch::la::Vector<Dimensions,double>&  faceCentre,
      const tarch::la::Vector<Dimensions,double>&  gridCellH,
      double                                       t,
      double                                       dt,
      int                                          normal,
      double * __restrict__ F // F[5]
    ) override;
    
     
     
    
    
   
    
 
 
    
 
 
    
     
   
        
};