CCZ4Kernels.h

Go to the documentation of this file.

//
// ExaHyPE2 CCZ4 implementation
//
#pragma once
 
 
#include <cmath>
#include "tarch/tarch.h"
#include "tarch/CompilerSpecificSettings.h"
#include "tarch/multicore/multicore.h"
#include "tarch/accelerator/accelerator.h"
 
 
namespace applications {
  namespace exahype2 {
    namespace ccz4 {

      #if defined(WITH_OPENMP)
      #pragma omp declare simd
      #endif
      static inline void ncp(double* BgradQ, const double* const Q, const double* const gradQSerialised, const int normal,
        const int CCZ4LapseType,
        const double CCZ4ds,
        const double CCZ4c,
        const double CCZ4e,
        const double CCZ4f,
        const double CCZ4bs,
        const double CCZ4sk,
        const double CCZ4xi,
        const double CCZ4mu,
        const double CCZ4SO
      );

      #if defined(WITH_OPENMP)
      #pragma omp declare simd
      #endif
      static inline void source(double* S, const double* const Q,
        const int CCZ4LapseType,
        const double CCZ4ds,
        const double CCZ4c,
        const double CCZ4e,
        const double CCZ4f,
        const double CCZ4bs,
        const double CCZ4sk,
        const double CCZ4xi,
        const double CCZ4itau,
        const double CCZ4eta,
        const double CCZ4k1,
        const double CCZ4k2,
        const double CCZ4k3,
        const double CCZ4SO
      );
 
      #if defined(WITH_OPENMP)
      #pragma omp declare simd
      #endif
      static inline double maxEigenvalue(
        const double* const Q,
        int                 normal,
        const double        CCZ4e,
        const double        CCZ4ds,
        const double        CCZ4GLMc,
        const double        CCZ4GLMd
      );
 

      #if defined(WITH_OPENMP)
      #pragma omp declare simd
      #endif
      static inline double maxEigenvalue(
        const double* const Q,
        int                 normal,
        const double        CCZ4e,
        const double        CCZ4ds,
        const double        CCZ4GLMc,
        const double        CCZ4GLMd
      );

      void admconstraints(double* constraints, const double* const Q, const double* const gradQSerialised);

      void ThetaOutputNCP(double* NCPterm, const double* const Q, const double* const gradQSerialised, int normal);
      void TestingOutput(double* terms, const double* const Q, const double* const gradQSerialised);

      void Psi4Calc(double* Psi4, const double* const Q, const double* const gradQSerialised, double* coor);
 

      void enforceCCZ4constraints(double* NOALIAS  newQ);

      void enforceCCZ4constraints(
        const double* NOALIAS  oldQ,
        double* NOALIAS        dQdt,
        double                      timeStepSize
      );
    }
  }
}
 
#include "CCZ4Kernels.cpph"