Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <cstddef>
       7              : 
       8              : #include "DataStructures/Tensor/TypeAliases.hpp"
       9              : 
      10              : /// \cond
      11              : class DataVector;
      12              : template <size_t Dim>
      13              : class Mesh;
      14              : template <size_t Dim>
      15              : class Direction;
      16              : namespace gsl {
      17              : template <typename T>
      18              : class not_null;
      19              : }  // namespace gsl
      20              : /// \endcond
      21              : 
      22              : namespace evolution::dg::subcell::fd {
      23              : /// @{
      24              : /*!
      25              :  * \brief Computes the logical coordinates of ghost grid points for a given
      26              :  * direction and ghost zone size.
      27              :  *
      28              :  * Let `d` be the axis dimension of the `direction`. The returned coordinate has
      29              :  * extents that is same as the volume mesh extents but `[d]`-th value
      30              :  * replaced by the ghost zone size.
      31              :  *
      32              :  * For instance if the (volume) subcell mesh has extents \f$(6,6,6)\f$, ghost
      33              :  * zone size is 2, and the `direction` is along Xi axis, the resulting
      34              :  * coordinates computed by this function has extents \f$(2,6,6)\f$.
      35              :  *
      36              :  */
      37              : template <size_t Dim>
      38            1 : void ghost_zone_logical_coordinates(
      39              :     const gsl::not_null<tnsr::I<DataVector, Dim, Frame::ElementLogical>*>
      40              :         ghost_logical_coords,
      41              :     const Mesh<Dim>& subcell_mesh, const size_t ghost_zone_size,
      42              :     const Direction<Dim>& direction);
      43              : 
      44              : template <size_t Dim>
      45            1 : tnsr::I<DataVector, Dim, Frame::ElementLogical> ghost_zone_logical_coordinates(
      46              :     const Mesh<Dim>& subcell_mesh, const size_t ghost_zone_size,
      47              :     const Direction<Dim>& direction);
      48              : /// @}
      49              : }  // namespace evolution::dg::subcell::fd