SelfSimilarInfallDG.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 "AbstractSelfSimilarInfallDG.h"
 
#include "tarch/logging/Log.h"
#include "tarch/accelerator/accelerator.h"
 
#include "spherical-accretion/MassAccumulator.h"
#include "spherical-accretion/GravityModel.h"
 
#include "EulerKernels.h"
 
 
namespace benchmarks {namespace exahype2 {namespace euler {namespace sphericalaccretionupscaling {
  class SelfSimilarInfallDG;
 
}}}}
 
 
class benchmarks::exahype2::euler::sphericalaccretionupscaling::SelfSimilarInfallDG: public AbstractSelfSimilarInfallDG {
  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;
 
    SelfSimilarInfallDG();
 
    void addDensity(
      const tarch::la::Vector<Dimensions,double>&  x,
      const tarch::la::Vector<3,double>&           cellSize,
      const tarch::la::Vector<3,int>&              index,
      double                                       density
    );
 
    
    void initialCondition(
      double * __restrict__ Q,
      const tarch::la::Vector<Dimensions,double>&  x,
      const tarch::la::Vector<Dimensions,double>&  h,
      const tarch::la::Vector<Dimensions,int>&     index,
      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>&  x,
      double                                       t,
      int                                          normal
    )  override;
    
 
    void sourceTerm(
      const double * __restrict__ Q,
      const tarch::la::Vector<Dimensions,double>&  x,
      double                                       t,
      double                                       dt,
      double * __restrict__ S
    ) override;
    
 
    virtual double maxEigenvalue(
      const double * __restrict__                  Q, // Q[5+0],
      const tarch::la::Vector<Dimensions,double>&  x,
      double                                       t,
      double                                       dt,
      int                                          normal
    ) override;
    
    
    virtual void flux(
      const double * __restrict__                  Q, // Q[5+0],
      const tarch::la::Vector<Dimensions,double>&  x,
      double                                       t,
      double                                       dt,
      int                                          normal,
      double * __restrict__ F // F[5]
    ) override;
    
     
    ::exahype2::RefinementCommand refinementCriterion(
      const double * __restrict__                  Q,    // Q[5+0]
      const tarch::la::Vector<Dimensions,double>&  x,
      const tarch::la::Vector<Dimensions,double>&  h,
      double                                       t
    ) override;
 
 
    void finishTimeStep() override;
 
    bool cellCanUseStatelessPDETerms(
      const tarch::la::Vector<Dimensions,double>&  patchCentre,
      const tarch::la::Vector<Dimensions,double>&  patchH,
      double                                       t,
      double                                       dt
    ) const override;
    
    
    #if defined(GPUOffloadingOMP)
    #pragma omp declare target
    #endif
    static GPUCallableMethod void sourceTerm(
       const double * __restrict__                  Q,
       const tarch::la::Vector<Dimensions,double>&  x,
       double                                       t,
       double                                       dt,
       double * __restrict__                        S,
       Offloadable
    ) {
      for (int i=0; i<NumberOfUnknowns; i++) {
        S[i] = 0.0;
      }
 
      applications::exahype2::euler::sphericalaccretion::addGravitationalSource_AlphaCDM(
        S,
        x,
        Q,
        TotalMassInPreviousTimeStep,
        aInitial,
        t
      );
    }
    #if defined(GPUOffloadingOMP)
    #pragma omp end declare target
    #endif
    
    #if defined(GPUOffloadingOMP)
    #pragma omp declare target
    #endif
    static GPUCallableMethod void flux(
      const double * __restrict__ Q,               // Q[5+0],
      const tarch::la::Vector<Dimensions,double>&  x,
      double                                       t,
      double                                       dt,
      int                                          normal,
      double * __restrict__                        F, // F[5]
      Offloadable
    ) {
      ::applications::exahype2::euler::flux(
        Q,
        x,
        0.0, // no volume, as point-wise
        t,
        dt,
        normal,
        NumberOfUnknowns,
        NumberOfAuxiliaryVariables,
        F,
        Gamma
      );
    }
    #if defined(GPUOffloadingOMP)
    #pragma omp end declare target
    #endif
};