SecondOrderAuxiliaryVariablesReconstruction.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 "exahype2/CellData.h"
#include "exahype2/enumerator/AoSLexicographicEnumerator.h"
 
 
namespace applications {
  namespace exahype2 {
    namespace ccz4 {
      void recomputeAuxiliaryVariablesFD4_4thOrder(
        ::exahype2::CellData<double, double>&   patchData,
        int                                     numberOfGridCellsPerPatchPerAxis,
        int                                     haloSize,
        int                                     unknowns,
        int                                     auxiliaryVariables
      );
 
      void recomputeAuxiliaryVariablesFD4_centralDifferences(
        ::exahype2::CellData<double, double>&   patchData,
        int                                     numberOfGridCellsPerPatchPerAxis,
        int                                     haloSize,
        int                                     unknowns,
        int                                     auxiliaryVariables
      );
 
      void recomputeAuxiliaryVariablesFD4_leftDifferences(
        ::exahype2::CellData<double, double>&   patchData,
        int                                     numberOfGridCellsPerPatchPerAxis,
        int                                     haloSize,
        int                                     unknowns,
        int                                     auxiliaryVariables
      );
 
      void recomputeAuxiliaryVariablesFD4_rightDifferences(
        ::exahype2::CellData<double, double>&   patchData,
        int                                     numberOfGridCellsPerPatchPerAxis,
        int                                     haloSize,
        int                                     unknowns,
        int                                     auxiliaryVariables
      );
 
      namespace internal {
        void recomputeAuxiliaryVariablesFD4_4thOrder_LoopBody(
          double* NOALIAS                                       QIn,
          const ::exahype2::enumerator::AoSLexicographicEnumerator& QInEnumerator,
          const tarch::la::Vector<DIMENSIONS,double>&  patchCentre,
          const tarch::la::Vector<DIMENSIONS,double>&  patchSize,
          int                                          patchIndex,
          const tarch::la::Vector<DIMENSIONS,int>&     volumeIndex,
          int                                          normal
        );
 
 
        void recomputeAuxiliaryVariablesFD4_centralDifferences_LoopBody(
          double* NOALIAS                                       QIn,
          const ::exahype2::enumerator::AoSLexicographicEnumerator& QInEnumerator,
          const tarch::la::Vector<DIMENSIONS,double>&  patchCentre,
          const tarch::la::Vector<DIMENSIONS,double>&  patchSize,
          int                                          patchIndex,
          const tarch::la::Vector<DIMENSIONS,int>&     volumeIndex,
          int                                          normal
        );
 
 
        void recomputeAuxiliaryVariablesFD4_leftDifferences_LoopBody(
          double* NOALIAS                                       QIn,
          const ::exahype2::enumerator::AoSLexicographicEnumerator& QInEnumerator,
          const tarch::la::Vector<DIMENSIONS,double>&  patchCentre,
          const tarch::la::Vector<DIMENSIONS,double>&  patchSize,
          int                                          patchIndex,
          const tarch::la::Vector<DIMENSIONS,int>&     volumeIndex,
          int                                          normal
        );
 
 
        void recomputeAuxiliaryVariablesFD4_rightDifferences_LoopBody(
          double* NOALIAS                                       QIn,
          const ::exahype2::enumerator::AoSLexicographicEnumerator& QInEnumerator,
          const tarch::la::Vector<DIMENSIONS,double>&  patchCentre,
          const tarch::la::Vector<DIMENSIONS,double>&  patchSize,
          int                                          patchIndex,
          const tarch::la::Vector<DIMENSIONS,int>&     volumeIndex,
          int                                          normal
        );
      }
    }
  }
}