SelfSimilarInfallFV.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 "../../../../applications/exahype2/euler/EulerKernels.h"
#include "AbstractSelfSimilarInfallFV.h"
#include "tarch/logging/Log.h"
 
#include "spherical-accretion/MassAccumulator.h"
#include "spherical-accretion/GravityModel.h"
 
 
 
namespace benchmarks {
  namespace exahype2 {
    namespace euler {
      namespace sphericalaccretionupscaling {
        class SelfSimilarInfallFV;
      }
    }
  }
}
 
 
class benchmarks::exahype2::euler::sphericalaccretionupscaling::SelfSimilarInfallFV: public benchmarks::exahype2::euler::sphericalaccretionupscaling::AbstractSelfSimilarInfallFV {
  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;
 
    SelfSimilarInfallFV();
 
    void addDensity(
      const tarch::la::Vector<Dimensions,double>&  volumeCentre,
      const tarch::la::Vector<Dimensions,double>&  volumeH,
      double                                       density
    );
 
 
    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,
      double * __restrict__                        Qoutside,
      const tarch::la::Vector<Dimensions,double>&  faceCentre,
      const tarch::la::Vector<Dimensions,double>&  volumeH,
      double                                       t,
      int                                          normal
    )  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;
    
 
    virtual double maxEigenvalue(
      const double * __restrict__                  Q,
      const tarch::la::Vector<Dimensions,double>&  faceCentre,
      const tarch::la::Vector<Dimensions,double>&  volumeH,
      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>&  volumeH,
      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>&  volumeCentre,
      const tarch::la::Vector<Dimensions,double>&  volumeH,
      double                                       t
    ) override;
 
 
    void finishTimeStep() override;
 
 
    #if defined(GPUOffloadingOMP)
    #pragma omp declare target
    #endif
    static GPUCallableMethod 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
    ) {
      for (int i=0; i<NumberOfUnknowns; i++) {
        S[i] = 0.0;
      }
 
      applications::exahype2::euler::sphericalaccretion::addGravitationalSource_AlphaCDM(
        S,
        volumeCentre,
        Q,
        TotalMassInPreviousTimeStep,
        aInitial,
        t
      );
    }
    #if defined(GPUOffloadingOMP)
    #pragma omp end declare target
    #endif
 
    
    #if defined(GPUOffloadingOMP)
    #pragma omp declare target
    #endif
    static GPUCallableMethod 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
    ) {
      return ::applications::exahype2::euler::maxEigenvalue(
          Q,
          faceCentre,
          volumeH,
          t,
          dt,
          normal,
          NumberOfUnknowns,
          NumberOfAuxiliaryVariables,
          Gamma
        );
    }
    #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>&  faceCentre,
      const tarch::la::Vector<Dimensions,double>&  volumeH,
      double                                       t,
      double                                       dt,
      int                                          normal,
      double * __restrict__                        F, // F[5]
      Offloadable
    ) {
      ::applications::exahype2::euler::flux(
        Q,
        faceCentre,
        volumeH,
        t,
        dt,
        normal,
        NumberOfUnknowns,
        NumberOfAuxiliaryVariables,
        F,
        Gamma
      );
    }
    #if defined(GPUOffloadingOMP)
    #pragma omp end declare target
    #endif
 
 
    bool patchCanUseStatelessPDETerms(
      const tarch::la::Vector<Dimensions,double>&  patchCentre,
      const tarch::la::Vector<Dimensions,double>&  patchH,
      double                                       t,
      double                                       dt
    ) const override;
 
//    Die Information was ausserhalb liegt, wird noch net im guard genutzt
};