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/System.hpp"
      14              : #include "Evolution/Systems/GrMhd/GhValenciaDivClean/Tags.hpp"
      15              : #include "Evolution/Systems/GrMhd/ValenciaDivClean/ConservativeFromPrimitive.hpp"
      16              : #include "Evolution/Systems/GrMhd/ValenciaDivClean/NewmanHamlin.hpp"
      17              : #include "Evolution/Systems/GrMhd/ValenciaDivClean/PrimitiveFromConservative.hpp"
      18              : #include "Evolution/Systems/GrMhd/ValenciaDivClean/System.hpp"
      19              : #include "Evolution/Systems/GrMhd/ValenciaDivClean/Tags.hpp"
      20              : #include "Evolution/Systems/GrMhd/ValenciaDivClean/TimeDerivativeTerms.hpp"
      21              : #include "Evolution/Systems/RadiationTransport/M1Grey/System.hpp"
      22              : #include "PointwiseFunctions/GeneralRelativity/Tags.hpp"
      23              : #include "PointwiseFunctions/Hydro/Tags.hpp"
      24              : #include "Utilities/TMPL.hpp"
      25              : 
      26              : /// Namespace associated with utilities for the combined Generalized Harmonic,
      27              : /// Valencia formulation of ideal GRMHD with divergence cleaning, and M1Grey
      28              : /// neutrino systems.
      29              : namespace grmhd::GhValenciaDivClean {
      30              : /// \cond
      31              : struct TimeDerivativeTerms;
      32              : /// \endcond
      33              : 
      34              : template <typename... NeutrinoSpecies>
      35            0 : struct System<
      36              :     RadiationTransport::M1Grey::System<tmpl::list<NeutrinoSpecies...>>> {
      37              :   // Will need M1 BCs on top of existing
      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 =
      44              :       RadiationTransport::M1Grey::System<tmpl::list<NeutrinoSpecies...>>;
      45              : 
      46              :   static_assert(std::is_same_v<Tags::spacetime_reconstruction_tags,
      47              :                                typename gh_system::variables_tag::tags_list>);
      48              : 
      49            0 :   using variables_tag = ::Tags::Variables<tmpl::append<
      50              :       typename gh_system::variables_tag::tags_list,
      51              :       typename grmhd_system::variables_tag::tags_list,
      52              :       typename neutrino_transport_system::variables_tag::tags_list>>;
      53            0 :   using non_conservative_variables =
      54              :       typename gh_system::variables_tag::tags_list;
      55            0 :   using flux_variables =
      56              :       tmpl::append<typename gh_system::flux_variables,
      57              :                    typename grmhd_system::flux_variables,
      58              :                    typename neutrino_transport_system::flux_variables>;
      59            0 :   using gradient_variables =
      60              :       tmpl::append<typename gh_system::gradient_variables,
      61              :                    typename grmhd_system::gradient_variables,
      62              :                    typename neutrino_transport_system::gradient_variables>;
      63            0 :   using gradients_tags = gradient_variables;
      64            0 :   static constexpr bool is_in_flux_conservative_form = false;
      65              : 
      66              :   // might need to filter duplicate tags
      67            0 :   using primitive_variables_tag =
      68              :       tmpl::append<typename grmhd_system::primitive_variables_tag,
      69              :                    typename neutrino_transport_system::primitive_variables_tag>;
      70              : 
      71              :   // contains required metric and metric derivatives for Gh system as well
      72            0 :   using spacetime_variables_tags =
      73              :       typename neutrino_transport_system::spacetime_variables_tag;
      74              : 
      75            0 :   using compute_volume_time_derivative_terms = TimeDerivativeTerms;
      76              : 
      77            0 :   using conservative_from_primitive =
      78              :       typename grmhd_system::conservative_from_primitive;
      79              :   template <typename OrderedListOfPrimitiveRecoverySchemes>
      80            0 :   using primitive_from_conservative =
      81              :       typename grmhd_system::template primitive_from_conservative<
      82              :           OrderedListOfPrimitiveRecoverySchemes>;
      83              : 
      84            0 :   using compute_largest_characteristic_speed =
      85              :       Tags::ComputeLargestCharacteristicSpeed<>;
      86              : 
      87            0 :   using inverse_spatial_metric_tag =
      88              :       typename gh_system::inverse_spatial_metric_tag;
      89              : };
      90              : }  // namespace grmhd::GhValenciaDivClean