Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include "DataStructures/Variables.hpp"
       7             : #include "Evolution/Systems/GrMhd/ValenciaDivClean/ConservativeFromPrimitive.hpp"
       8             : #include "Evolution/Systems/GrMhd/ValenciaDivClean/Fluxes.hpp"
       9             : #include "Utilities/Gsl.hpp"
      10             : 
      11             : namespace grmhd::ValenciaDivClean {
      12             : namespace detail {
      13             : template <typename FluxTags, typename BoundaryVarsType, typename... ReturnTags,
      14             :           typename... ConservativeTags, typename... ArgumentTags,
      15             :           typename... FluxArgumentTags>
      16             : void compute_fluxes_from_primitives_impl(
      17             :     const gsl::not_null<Variables<FluxTags>*> flux_vars,
      18             :     const BoundaryVarsType& boundary_vars, tmpl::list<ReturnTags...> /*meta*/,
      19             :     tmpl::list<ConservativeTags...> /*meta*/,
      20             :     tmpl::list<ArgumentTags...> /*meta*/,
      21             :     tmpl::list<FluxArgumentTags...> /*meta*/) {
      22             :   Variables<typename tmpl::list<ConservativeTags...>> conserved_vars{
      23             :       flux_vars->number_of_grid_points()};
      24             :   grmhd::ValenciaDivClean::ConservativeFromPrimitive::apply(
      25             :       make_not_null(&get<ConservativeTags>(conserved_vars))...,
      26             :       get<ArgumentTags>(boundary_vars)...);
      27             :   grmhd::ValenciaDivClean::ComputeFluxes::apply(
      28             :       make_not_null(&get<ReturnTags>(*flux_vars))...,
      29             :       get<ConservativeTags>(conserved_vars)...,
      30             :       get<FluxArgumentTags>(boundary_vars)...);
      31             : }
      32             : }  // namespace detail
      33             : 
      34             : /*!
      35             :  * \brief Helper function that computes fluxes from a given set of primitive
      36             :  * variables. Primarily used for filling neighbor data when flux information is
      37             :  * needed for boundaries.
      38             :  *
      39             :  * \warning This computes the conservative variables internally. If you also
      40             :  * need those, you shouldn't call this function.
      41             :  */
      42             : template <typename FluxTags, typename BoundaryVarsType>
      43           1 : void compute_fluxes_from_primitives(
      44             :     const gsl::not_null<Variables<FluxTags>*> flux_vars,
      45             :     const BoundaryVarsType& boundary_vars) {
      46             :   using ConservativeTags =
      47             :       typename grmhd::ValenciaDivClean::ConservativeFromPrimitive::return_tags;
      48             :   using FluxArgumentTags =
      49             :       typename grmhd::ValenciaDivClean::ComputeFluxes::argument_tags;
      50             :   using flux_arg_tags = tmpl::back<tmpl::split_at<
      51             :       FluxArgumentTags, tmpl::next<tmpl::index_of<
      52             :                             FluxArgumentTags, tmpl::back<ConservativeTags>>>>>;
      53             :   detail::compute_fluxes_from_primitives_impl(
      54             :       flux_vars, boundary_vars,
      55             :       typename grmhd::ValenciaDivClean::ComputeFluxes::return_tags{},
      56             :       ConservativeTags{},
      57             :       typename grmhd::ValenciaDivClean::ConservativeFromPrimitive::
      58             :           argument_tags{},
      59             :       flux_arg_tags{});
      60             : }
      61             : }  // namespace grmhd::ValenciaDivClean