PatchUtils.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 "tarch/la/Vector.h"
#include "tarch/multicore/multicore.h"
 
#include "peano4/grid/GridControlEvent.h"
#include "peano4/utils/Globals.h"
 
#include "exahype2/fv/PatchUtils.h"
 
 
namespace exahype2 {
  namespace fd {
    static tarch::la::Vector<2,double>  getGridCellSize(
      const tarch::la::Vector<2,double>&  h,
      int                                 numberOfGridCellsPerPatchPerAxis
    )
    //InlineMethod
    {
      return ::exahype2::fv::getVolumeSize(h,numberOfGridCellsPerPatchPerAxis);
    }
 
 
    static tarch::la::Vector<3,double>  getGridCellSize(
      const tarch::la::Vector<3,double>&  h,
      int                                 numberOfGridCellsPerPatchPerAxis
    )
    //InlineMethod
    {
      return ::exahype2::fv::getVolumeSize(h,numberOfGridCellsPerPatchPerAxis);
    }
 
 
    static tarch::la::Vector<2,double>  getGridFaceSize(
      const tarch::la::Vector<2,double>&  h,
      int                                 numberOfGridCellsPerPatchPerAxis
    )
    //InlineMethod
    {
      return ::exahype2::fv::getFaceSize(h,numberOfGridCellsPerPatchPerAxis);
    }
 
 
    static tarch::la::Vector<3,double>  getGridFaceSize(
      const tarch::la::Vector<3,double>&  h,
      int                                 numberOfGridCellsPerPatchPerAxis
    )
    //InlineMethod
    {
      return ::exahype2::fv::getFaceSize(h,numberOfGridCellsPerPatchPerAxis);
    }
 
 
    static tarch::la::Vector<2,double>  getGridCellCentre(
      const tarch::la::Vector<2,double>&  x,
      const tarch::la::Vector<2,double>&  h,
      int                                 numberOfGridCellsPerPatchPerAxis,
      const tarch::la::Vector<2,int>&     index
    )
    //InlineMethod
    {
      return ::exahype2::fv::getVolumeCentre(x,h,numberOfGridCellsPerPatchPerAxis,index);
    }
 
 
    static tarch::la::Vector<3,double>  getGridCellCentre(
      const tarch::la::Vector<3,double>&  x,
      const tarch::la::Vector<3,double>&  h,
      int                                 numberOfGridCellsPerPatchPerAxis,
      const tarch::la::Vector<3,int>&     index
    )
    //InlineMethod
    {
      return ::exahype2::fv::getVolumeCentre(x,h,numberOfGridCellsPerPatchPerAxis,index);
    }
 
 
    static tarch::la::Vector<2,double>  getGridFaceCentre(
      const tarch::la::Vector<2,double>&  x,
      const tarch::la::Vector<2,double>&  h,
      int                                 numberOfGridCellsPerPatchPerAxis,
      int                                 overlap,
      int                                 normal,
      const tarch::la::Vector<2,int>&     index
    )
    //InlineMethod
    {
      return ::exahype2::fv::getFaceCentre(x,h,numberOfGridCellsPerPatchPerAxis,overlap,normal,index);
    }
 
 
    static tarch::la::Vector<3,double>  getGridFaceCentre(
      const tarch::la::Vector<3,double>&  x,
      const tarch::la::Vector<3,double>&  h,
      int                                 numberOfGridCellsPerPatchPerAxis,
      int                                 overlap,
      int                                 normal,
      const tarch::la::Vector<3,int>&     index
    )
    //InlineMethod
    {
      return ::exahype2::fv::getFaceCentre(x,h,numberOfGridCellsPerPatchPerAxis,overlap,normal,index);
    }
 
 
    std::string plotGridCell(
      const double* __restrict__ Q,
      int    unknowns
    );
 
    void validatePatch(
      const double* __restrict__ Q,
      int    unknowns,
      int    auxiliaryVariables,
      int    numberOfGridCellsPerPatchPerAxis,
      int    haloSize,
      const std::string& location = "",
      bool   triggerNonCriticalAssertion = true,
      double* minValues = nullptr,
      double* maxValues = nullptr
    );
 
    std::string plotPatch(
      const double* __restrict__ Q,
      int    unknowns,
      int    auxiliaryVariables,
      int    numberOfGridCellsPerPatchPerAxis,
      int    haloSize,
      bool   prettyPrint = false
    );
 
 
    std::string plotPatchOverlap(
      const double* __restrict__ Q,
      int    unknowns,
      int    auxiliaryVariables,
      int    numberOfGridCellsPerPatchPerAxis,
      int    haloSize,
      int    normal,
      bool   prettyPrint = false
    );
  }
}