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/EagerMath/Magnitude.hpp"
       9             : #include "DataStructures/VariablesTag.hpp"
      10             : #include "Evolution/Systems/GeneralizedHarmonic/System.hpp"
      11             : #include "Evolution/Systems/GrMhd/GhValenciaDivClean/BoundaryConditions/BoundaryCondition.hpp"
      12             : #include "Evolution/Systems/GrMhd/GhValenciaDivClean/Characteristics.hpp"
      13             : #include "Evolution/Systems/GrMhd/GhValenciaDivClean/Tags.hpp"
      14             : #include "Evolution/Systems/GrMhd/ValenciaDivClean/ConservativeFromPrimitive.hpp"
      15             : #include "Evolution/Systems/GrMhd/ValenciaDivClean/NewmanHamlin.hpp"
      16             : #include "Evolution/Systems/GrMhd/ValenciaDivClean/PrimitiveFromConservative.hpp"
      17             : #include "Evolution/Systems/GrMhd/ValenciaDivClean/System.hpp"
      18             : #include "Evolution/Systems/GrMhd/ValenciaDivClean/Tags.hpp"
      19             : #include "Evolution/Systems/GrMhd/ValenciaDivClean/TimeDerivativeTerms.hpp"
      20             : #include "PointwiseFunctions/GeneralRelativity/Tags.hpp"
      21             : #include "PointwiseFunctions/Hydro/Tags.hpp"
      22             : #include "Utilities/TMPL.hpp"
      23             : 
      24             : namespace grmhd {
      25             : 
      26             : /// Namespace associated with utilities for the combined Generalized Harmonic
      27             : /// and Valencia formulation of ideal GRMHD with divergence cleaning systems.
      28             : namespace GhValenciaDivClean {
      29             : /// \cond
      30             : struct TimeDerivativeTerms;
      31             : /// \endcond
      32             : 
      33             : template <typename NeutrinoTransportSystem>
      34             : struct System;
      35             : 
      36             : template <typename NeutrinoTransportSystem>
      37           0 : struct System {
      38           0 :   using boundary_conditions_base = BoundaryConditions::BoundaryCondition;
      39           0 :   static constexpr bool has_primitive_and_conservative_vars = true;
      40           0 :   static constexpr size_t volume_dim = 3;
      41           0 :   using grmhd_system = grmhd::ValenciaDivClean::System;
      42           0 :   using gh_system = gh::System<3_st>;
      43           0 :   using neutrino_transport_system = NeutrinoTransportSystem;
      44             : 
      45             :   static_assert(std::is_same_v<Tags::spacetime_reconstruction_tags,
      46             :                                typename gh_system::variables_tag::tags_list>);
      47             : 
      48           0 :   using variables_tag = ::Tags::Variables<
      49             :       tmpl::append<typename gh_system::variables_tag::tags_list,
      50             :                    typename grmhd_system::variables_tag::tags_list>>;
      51           0 :   using non_conservative_variables =
      52             :       typename gh_system::variables_tag::tags_list;
      53           0 :   using flux_variables = tmpl::append<typename gh_system::flux_variables,
      54             :                                       typename grmhd_system::flux_variables>;
      55           0 :   using gradient_variables =
      56             :       tmpl::append<typename gh_system::gradient_variables,
      57             :                    typename grmhd_system::gradient_variables>;
      58           0 :   using gradients_tags = gradient_variables;
      59           0 :   static constexpr bool is_in_flux_conservative_form = false;
      60             : 
      61           0 :   using primitive_variables_tag =
      62             :       typename grmhd_system::primitive_variables_tag;
      63           0 :   using spacetime_variables_tag = ::Tags::Variables<
      64             :       tmpl::list<::Tags::deriv<gr::Tags::Lapse<DataVector>, tmpl::size_t<3>,
      65             :                                Frame::Inertial>,
      66             :                  ::Tags::deriv<gr::Tags::Shift<DataVector, 3>, tmpl::size_t<3>,
      67             :                                Frame::Inertial>,
      68             :                  ::Tags::deriv<gr::Tags::SpatialMetric<DataVector, 3>,
      69             :                                tmpl::size_t<3>, Frame::Inertial>,
      70             :                  gr::Tags::ExtrinsicCurvature<DataVector, 3>>>;
      71             : 
      72           0 :   using compute_volume_time_derivative_terms = TimeDerivativeTerms;
      73             : 
      74           0 :   using conservative_from_primitive =
      75             :       typename grmhd_system::conservative_from_primitive;
      76             :   template <typename OrderedListOfPrimitiveRecoverySchemes>
      77           0 :   using primitive_from_conservative =
      78             :       typename grmhd_system::template primitive_from_conservative<
      79             :           OrderedListOfPrimitiveRecoverySchemes>;
      80             : 
      81           0 :   using compute_largest_characteristic_speed =
      82             :       Tags::ComputeLargestCharacteristicSpeed<>;
      83             : 
      84           0 :   using inverse_spatial_metric_tag =
      85             :       typename gh_system::inverse_spatial_metric_tag;
      86             : };
      87             : }  // namespace GhValenciaDivClean
      88             : }  // namespace grmhd