Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <array>
       7             : #include <cstddef>
       8             : #include <utility>
       9             : 
      10             : #include "Domain/Structure/DirectionalIdMap.hpp"
      11             : 
      12             : /// \cond
      13             : class DataVector;
      14             : template <typename TagsList>
      15             : class Variables;
      16             : namespace gsl {
      17             : template <typename>
      18             : class not_null;
      19             : }  // namespace gsl
      20             : template <size_t Dim>
      21             : class Direction;
      22             : template <size_t Dim>
      23             : class ElementId;
      24             : template <size_t Dim>
      25             : class Element;
      26             : template <size_t Dim>
      27             : class Mesh;
      28             : namespace EquationsOfState {
      29             : template <bool IsRelativistic, size_t ThermodynamicDim>
      30             : class EquationOfState;
      31             : }  // namespace EquationsOfState
      32             : namespace evolution::dg::subcell {
      33             : class GhostData;
      34             : }  // namespace evolution::dg::subcell
      35             : /// \endcond
      36             : 
      37             : namespace NewtonianEuler::fd {
      38             : /*!
      39             :  * \brief Reconstructs the mass density, velocity, and pressure, then computes
      40             :  * the specific internal energy and conserved variables. All results are written
      41             :  * into `vars_on_lower_face` and `vars_on_upper_face`.
      42             :  */
      43             : template <typename PrimsTags, typename TagsList, size_t Dim, typename F>
      44           1 : void reconstruct_prims_work(
      45             :     gsl::not_null<std::array<Variables<TagsList>, Dim>*> vars_on_lower_face,
      46             :     gsl::not_null<std::array<Variables<TagsList>, Dim>*> vars_on_upper_face,
      47             :     const F& reconstruct, const Variables<PrimsTags>& volume_prims,
      48             :     const EquationsOfState::EquationOfState<false, 2>& eos,
      49             :     const Element<Dim>& element,
      50             :     const DirectionalIdMap<Dim, evolution::dg::subcell::GhostData>& ghost_data,
      51             :     const Mesh<Dim>& subcell_mesh, size_t ghost_zone_size);
      52             : 
      53             : /*!
      54             :  * \brief Reconstructs the mass density, velocity, and pressure, then computes
      55             :  * the specific internal energy and conserved variables. All results are written
      56             :  * into `vars_on_face`.
      57             :  *
      58             :  * This is used on DG elements to reconstruct their subcell neighbors' solution
      59             :  * on the shared faces.
      60             :  */
      61             : template <typename TagsList, typename PrimsTags, size_t Dim, typename F0,
      62             :           typename F1>
      63           1 : void reconstruct_fd_neighbor_work(
      64             :     gsl::not_null<Variables<TagsList>*> vars_on_face,
      65             :     const F0& reconstruct_lower_neighbor, const F1& reconstruct_upper_neighbor,
      66             :     const Variables<PrimsTags>& subcell_volume_prims,
      67             :     const EquationsOfState::EquationOfState<false, 2>& eos,
      68             :     const Element<Dim>& element,
      69             :     const DirectionalIdMap<Dim, evolution::dg::subcell::GhostData>& ghost_data,
      70             :     const Mesh<Dim>& subcell_mesh,
      71             :     const Direction<Dim>& direction_to_reconstruct, size_t ghost_zone_size);
      72             : }  // namespace NewtonianEuler::fd