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 :   static std::string name() { return "GhValenciaDivClean"; }
      39              : 
      40            0 :   using boundary_conditions_base = BoundaryConditions::BoundaryCondition;
      41            0 :   static constexpr bool has_primitive_and_conservative_vars = true;
      42            0 :   static constexpr size_t volume_dim = 3;
      43            0 :   using grmhd_system = grmhd::ValenciaDivClean::System;
      44            0 :   using gh_system = gh::System<3_st>;
      45            0 :   using neutrino_transport_system = NeutrinoTransportSystem;
      46              : 
      47              :   static_assert(std::is_same_v<Tags::spacetime_reconstruction_tags,
      48              :                                typename gh_system::variables_tag::tags_list>);
      49              : 
      50            0 :   using variables_tag = ::Tags::Variables<
      51              :       tmpl::append<typename gh_system::variables_tag::tags_list,
      52              :                    typename grmhd_system::variables_tag::tags_list>>;
      53            0 :   using non_conservative_variables =
      54              :       typename gh_system::variables_tag::tags_list;
      55            0 :   using flux_variables = tmpl::append<typename gh_system::flux_variables,
      56              :                                       typename grmhd_system::flux_variables>;
      57            0 :   using gradient_variables =
      58              :       tmpl::append<typename gh_system::gradient_variables,
      59              :                    typename grmhd_system::gradient_variables>;
      60            0 :   using gradients_tags = gradient_variables;
      61            0 :   static constexpr bool is_in_flux_conservative_form = false;
      62              : 
      63            0 :   using primitive_variables_tag =
      64              :       typename grmhd_system::primitive_variables_tag;
      65            0 :   using spacetime_variables_tag = ::Tags::Variables<
      66              :       tmpl::list<::Tags::deriv<gr::Tags::Lapse<DataVector>, tmpl::size_t<3>,
      67              :                                Frame::Inertial>,
      68              :                  ::Tags::deriv<gr::Tags::Shift<DataVector, 3>, tmpl::size_t<3>,
      69              :                                Frame::Inertial>,
      70              :                  ::Tags::deriv<gr::Tags::SpatialMetric<DataVector, 3>,
      71              :                                tmpl::size_t<3>, Frame::Inertial>,
      72              :                  gr::Tags::ExtrinsicCurvature<DataVector, 3>>>;
      73              : 
      74            0 :   using compute_volume_time_derivative_terms = TimeDerivativeTerms;
      75              : 
      76            0 :   using conservative_from_primitive =
      77              :       typename grmhd_system::conservative_from_primitive;
      78              :   template <typename OrderedListOfPrimitiveRecoverySchemes>
      79            0 :   using primitive_from_conservative =
      80              :       typename grmhd_system::template primitive_from_conservative<
      81              :           OrderedListOfPrimitiveRecoverySchemes>;
      82              : 
      83            0 :   using compute_largest_characteristic_speed =
      84              :       Tags::ComputeLargestCharacteristicSpeed<>;
      85              : 
      86            0 :   using inverse_spatial_metric_tag =
      87              :       typename gh_system::inverse_spatial_metric_tag;
      88              : };
      89              : }  // namespace GhValenciaDivClean
      90              : }  // namespace grmhd