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/Variables.hpp"
       9             : #include "Evolution/Systems/ScalarAdvection/Fluxes.hpp"
      10             : #include "Utilities/Gsl.hpp"
      11             : #include "Utilities/TMPL.hpp"
      12             : 
      13           1 : namespace ScalarAdvection::subcell {
      14             : namespace detail {
      15             : /*!
      16             :  * \brief Helper function that calls `ScalarAdvection::Fluxes::apply` by
      17             :  * retrieving the return and argument tags from the Variables object `vars`.
      18             :  */
      19             : template <size_t Dim, typename TagsList, typename... ReturnTags,
      20             :           typename... ArgumentTags>
      21             : void compute_fluxes_impl(const gsl::not_null<Variables<TagsList>*> vars,
      22             :                          tmpl::list<ReturnTags...> /*meta*/,
      23             :                          tmpl::list<ArgumentTags...> /*meta*/) {
      24             :   Fluxes<Dim>::apply(make_not_null(&get<ReturnTags>(*vars))...,
      25             :                      get<ArgumentTags>(*vars)...);
      26             : }
      27             : }  // namespace detail
      28             : 
      29             : /*!
      30             :  * \brief Compute fluxes for subcell variables.
      31             :  */
      32             : template <size_t Dim, typename TagsList>
      33           1 : void compute_fluxes(const gsl::not_null<Variables<TagsList>*> vars) {
      34             :   detail::compute_fluxes_impl<Dim>(vars, typename Fluxes<Dim>::return_tags{},
      35             :                                    typename Fluxes<Dim>::argument_tags{});
      36             : }
      37             : }  // namespace ScalarAdvection::subcell