CCZ4.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 "AbstractCCZ4.h"
#include "tarch/logging/Log.h"
 
 
#include <cstring>
#include "InitialValues.h"
 
#include "CCZ4Kernels.h"
 
 
namespace benchmarks {
  namespace exahype2 {
    namespace ccz4 {
      class CCZ4;
    }
  }
}
 
 
class benchmarks::exahype2::ccz4::CCZ4: public AbstractCCZ4 {
  private:
    static tarch::logging::Log   _log;
 
  public:
    CCZ4();
    
    ::exahype2::RefinementCommand refinementCriterion(
      const double * NOALIAS Q, // Q[59+0],
      const tarch::la::Vector<DIMENSIONS,double>&  volumeCentre,
      const tarch::la::Vector<DIMENSIONS,double>&  volumeH,
      double                                       t
    ) override;
    
#ifdef IncludeTwoPunctures
    void prepare(TP::TwoPunctures* tp); 
#endif
 
    void initialCondition(
      double * NOALIAS Q,
      const tarch::la::Vector<DIMENSIONS,double>&  volumeCentre,
      const tarch::la::Vector<DIMENSIONS,double>&  volumeH
    ) override;
    
 
    
    
    virtual void boundaryConditions(
      const double * NOALIAS Qinside, // Qinside[59+0]
      double * NOALIAS Qoutside, // Qoutside[59+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 * NOALIAS Q,
      const tarch::la::Vector<DIMENSIONS,double>&  volumeCentre,
      const tarch::la::Vector<DIMENSIONS,double>&  volumeH,
      double                                       t,
      double                                       dt,
      double * NOALIAS S
    ) override;
    
 
    
    virtual double maxEigenvalue(
      const double * NOALIAS Q, // Q[59+0],
      const tarch::la::Vector<DIMENSIONS,double>&  faceCentre,
      const tarch::la::Vector<DIMENSIONS,double>&  volumeH,
      double                                       t,
      double                                       dt,
      int                                          normal
    ) override;
    
 
 
    
     
     
    
    virtual void nonconservativeProduct(
      const double * NOALIAS Q, // Q[59+0],
      const double * NOALIAS             deltaQ, // [59+0]
      const tarch::la::Vector<DIMENSIONS,double>&  faceCentre,
      const tarch::la::Vector<DIMENSIONS,double>&  volumeH,
      double                                       t,
      double                                       dt,
      int                                          normal,
      double * NOALIAS BgradQ // BgradQ[59]
    ) override;
    
    
   
    
    static GPU_CALLABLE_METHOD void sourceTerm(
       const double * NOALIAS Q,
       const tarch::la::Vector<DIMENSIONS,double>&  volumeCentre,
       const tarch::la::Vector<DIMENSIONS,double>&  volumeH,
       double                                       t,
       double                                       dt,
       double * NOALIAS S,
       ::exahype2::Solver::Offloadable
     ) {
      memset(S, 0, NumberOfUnknowns*sizeof(double));
      applications::exahype2::ccz4::source(S,Q, CCZ4LapseType, CCZ4ds, CCZ4c, CCZ4e, CCZ4f, CCZ4bs, CCZ4sk, CCZ4xi, CCZ4itau, CCZ4eta, CCZ4k1, CCZ4k2, CCZ4k3, CCZ4SO);
    }
 
    static GPU_CALLABLE_METHOD double maxEigenvalue(
      const double * NOALIAS Q, // Q[59+0],
      const tarch::la::Vector<DIMENSIONS,double>&  faceCentre,
      const tarch::la::Vector<DIMENSIONS,double>&  volumeH,
      double                                       t,
      double                                       dt,
      int                                          normal,
      ::exahype2::Solver::Offloadable
    ) {
      const double np = 1.0;
      #if defined(CCZ4EINSTEIN)
      const double qmin = std::min({Q[0],Q[3],Q[5]});
      const double alpha = std::max({1.0, Q[16]}) * std::max({1.0, pow(Q[54],1.0/np)}) / std::sqrt(qmin);
      #else
      const double alpha = 1.0;
      #endif
 
      constexpr double sqrtwo = 1.4142135623730951;
      // NOTE parameters are stored in abstract superclass
      const double tempA = alpha * std::max({sqrtwo, CCZ4e, CCZ4ds, CCZ4GLMc/alpha, CCZ4GLMd/alpha});
      const double tempB = Q[17+normal];//DOT_PRODUCT(Q(18:20),nv(:))
      double tem = std::max({1.0, std::abs(-tempA-tempB), std::abs(tempA-tempB)});
      bool flag = false;
      return tem;
    }
 
    static GPU_CALLABLE_METHOD void nonconservativeProduct(
      const double * NOALIAS                  Q,         // Q[5+0],
      const double * NOALIAS                  deltaQ,    // [5+0]
      const tarch::la::Vector<DIMENSIONS,double>&  faceCentre,
      const tarch::la::Vector<DIMENSIONS,double>&  volumeH,
      double                                       t,
      double                                       dt,
      int                                          normal,
      double * NOALIAS                        BgradQ,     // BgradQ[5]
      ::exahype2::Solver::Offloadable
    );
};