Riemann.h

Go to the documentation of this file.

// This file is part of the ExaHyPE2 project. For conditions of distribution and
// use, please see the copyright notice at www.peano-framework.org
#pragma once
 
#include <vector>
 
#include "DGUtils.h"
#include "Functors.h"
 
#include "tarch/la/Vector.h"
#include "tarch/multicore/multicore.h"
 
#include "peano4/grid/GridControlEvent.h"
#include "peano4/utils/Loop.h"
#include "peano4/utils/Globals.h"
 
namespace exahype2 {
  namespace dg {
    void interpolateRiemannSolution(
      const peano4::datamanagement::FaceMarker&  marker,
      int                                        order,
      int                                        unknowns,
      const double* __restrict__                 InterpolationMatrix1d,
      const double* __restrict__                 coarseGridFaceQ,
      double* __restrict__                       fineGridFaceQ
    );
 
    void restrictAndAccumulateProjectedFacePolynomial(
      const peano4::datamanagement::FaceMarker&  marker,
      int                                        order,
      int                                        numberOfProjectedQuantities,
      int                                        unknowns,
      int                                        auxiliaryVariables,
      const double* __restrict__                 RestrictionMatrix1d,
      const double* __restrict__                 fineGridFaceQ,
      double* __restrict__                       coarseGridFaceQ
    );
 
    void projectVolumetricDataOntoFaces(
      const double* __restrict__          cellQ,
      int                                 order,
      int                                 unknowns,
      int                                 auxiliaryVariables,
      const double* const __restrict__    BasisFunctionValuesLeft,
      double* const __restrict__          faceQLeft,
      double* const __restrict__          faceQRight,
      double* const __restrict__          faceQBottom,
      double* const __restrict__          faceQUp
    );
 
    void clearSolutionProjection(
      int                           order,
      int                           unknowns,
      int                           auxiliaryVariables,
      int                           numberOfProjectedQuantities,
      double* __restrict__          faceQ
    );
 
    void clearRiemannResult(
      int                           order,
      int                           unknowns,
      double* __restrict__          faceQ
    );
 
    void projectVolumetricDataOntoFaces(
      const double* __restrict__          cellQ,
      int                                 order,
      int                                 unknowns,
      int                                 auxiliaryVariables,
      const double* const __restrict__    BasisFunctionValuesLeft,
      double* const __restrict__          faceQLeft,
      double* const __restrict__          faceQRight,
      double* const __restrict__          faceQBottom,
      double* const __restrict__          faceQUp,
      double* const __restrict__          faceQFront,
      double* const __restrict__          faceQBack
    );
 
    void projectVolumetricDataAndGradientOntoFaces(
      const double* __restrict__          cellQ,
      int                                 order,
      int                                 unknowns,
      int                                 auxiliaryVariables,
      const double* const __restrict__    BasisFunctionValuesLeft,
      double* const __restrict__          faceQLeft,
      double* const __restrict__          faceQRight,
      double* const __restrict__          faceQBottom,
      double* const __restrict__          faceQUp
    );
 
    void projectVolumetricDataAndGradientOntoFaces(
      const double* __restrict__          cellQ,
      int                                 order,
      int                                 unknowns,
      int                                 auxiliaryVariables,
      const double* const __restrict__    BasisFunctionValuesLeft,
      double* const __restrict__          faceQLeft,
      double* const __restrict__          faceQRight,
      double* const __restrict__          faceQBottom,
      double* const __restrict__          faceQUp,
      double* const __restrict__          faceQFront,
      double* const __restrict__          faceQBack
    );
 
    void integrateOverRiemannSolutionsAndAddToVolume_GaussLegendre(
      const double* const __restrict__    faceQLeft,
      const double* const __restrict__    faceQRight,
      const double* const __restrict__    faceQBottom,
      const double* const __restrict__    faceQUp,
      int                                 order,
      int                                 unknowns,
      const int                           auxiliaryVariables,
      const tarch::la::Vector<2,double>&  cellSize,
      const double* const __restrict__    BasisFunctionValuesLeft,
      const double* __restrict__          MassMatrixDiagonal1d,
      double* __restrict__                cellQ
    );
 
    void integrateOverRiemannSolutionsAndAddToVolume_GaussLegendre(
      const double* const __restrict__    faceQLeft,
      const double* const __restrict__    faceQRight,
      const double* const __restrict__    faceQBottom,
      const double* const __restrict__    faceQUp,
      const double* const __restrict__    faceQFront,
      const double* const __restrict__    faceQBack,
      int                                 order,
      int                                 unknowns,
      const int                           auxiliaryVariables,
      const tarch::la::Vector<3,double>&  cellSize,
      const double* const __restrict__    BasisFunctionValuesLeft,
      const double* __restrict__          MassMatrixDiagonal1d,
      double* __restrict__                cellQ
    );
  }
}