Kernels.h

Go to the documentation of this file.

// This file is part of the ExaGRyPE project. For conditions of distribution and
// use, please see the copyright notice at www.peano-framework.org
#pragma once
 
 
#include <cmath>
#include "tarch/tarch.h"
#include "tarch/compiler/CompilerSpecificSettings.h"
#include "tarch/multicore/multicore.h"
#include "tarch/accelerator/accelerator.h"
#include "tarch/la/Vector.h"
 
 
namespace exagrype {
  namespace ccz4 {
 
      #if defined(GPUOffloadingOMP)
      #pragma omp declare target
      #endif
      #if defined(SharedOMP)
      #pragma omp declare simd
      #endif
      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(SharedOMP)
      #pragma omp declare simd
      #endif
      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(SharedOMP)
      #pragma omp declare simd
      #endif
      void maxEigenvalue(
        const double* const Q,
        int                 normal,
        const double        CCZ4e,
        const double        CCZ4ds,
        const double        CCZ4GLMc,
        const double        CCZ4GLMd,
        double*             maxEigenvalue
      );
      #if defined(GPUOffloadingOMP)
      #pragma omp end declare target
      #endif
 
      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* __restrict__  newQ);
 
      void addKreissOligerDissipation(
        const double* __restrict__          QWithHalo,
        double* __restrict__                rightHandSide,
        const tarch::la::Vector<3, double>& patchSize,
        int                                 numberOfGridCellsPerPatchPerAxis,
        int                                 haloSize,
        int                                 unknowns,
        int                                 auxiliaryVariables,
        double                              sigma,
        const int*                          quantities,
        int                                 numberOfQuantities
      );
  }
}