Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <cstdint>
       7             : #include <vector>
       8             : 
       9             : #include "Evolution/Actions/RunEventsAndTriggers.hpp"
      10             : #include "Evolution/DiscontinuousGalerkin/EqualRateLts/ChangeFixedLtsRatio.hpp"
      11             : #include "Evolution/DiscontinuousGalerkin/Initialization/ProjectSpectralFilters.hpp"
      12             : #include "Evolution/Executables/GeneralizedHarmonic/GeneralizedHarmonicBase.hpp"
      13             : #include "Options/Protocols/FactoryCreation.hpp"
      14             : #include "Options/String.hpp"
      15             : #include "Parallel/ArrayCollection/DgElementCollection.hpp"
      16             : #include "Parallel/MemoryMonitor/MemoryMonitor.hpp"
      17             : #include "Parallel/PhaseControl/PhaseControlTags.hpp"
      18             : #include "Parallel/Protocols/RegistrationMetavariables.hpp"
      19             : #include "ParallelAlgorithms/Actions/MutateApply.hpp"
      20             : #include "ParallelAlgorithms/Amr/Projectors/CopyFromCreatorOrLeaveAsIs.hpp"
      21             : #include "PointwiseFunctions/AnalyticSolutions/GeneralRelativity/GaugeWave.hpp"
      22             : #include "Time/Actions/SelfStartActions.hpp"
      23             : #include "Time/AdvanceTime.hpp"
      24             : #include "Time/ChangeSlabSize/Action.hpp"
      25             : #include "Time/ChangeSlabSize/Tags.hpp"
      26             : #include "Time/Tags/StepperErrors.hpp"
      27             : #include "Utilities/ErrorHandling/Error.hpp"
      28             : #include "Utilities/ProtocolHelpers.hpp"
      29             : 
      30             : template <size_t VolumeDim>
      31             : struct EvolutionMetavars : public GeneralizedHarmonicTemplateBase<VolumeDim> {
      32             :   static constexpr size_t volume_dim = VolumeDim;
      33           0 :   using gh_base = GeneralizedHarmonicTemplateBase<volume_dim>;
      34             :   using typename gh_base::const_global_cache_tags;
      35             :   using typename gh_base::dg_registration_list;
      36           0 :   using initialization_actions =
      37             :       typename gh_base::template initialization_actions<EvolutionMetavars,
      38             :                                                         false>;
      39             :   using typename gh_base::initialize_initial_data_dependent_quantities_actions;
      40             :   using typename gh_base::observed_reduction_data_tags;
      41             :   using typename gh_base::system;
      42             :   static constexpr bool use_dg_element_collection =
      43             :       gh_base::use_dg_element_collection;
      44             : 
      45           0 :   using step_actions =
      46             :       typename gh_base::template step_actions<EvolutionMetavars, tmpl::list<>>;
      47             : 
      48           0 :   using gh_dg_element_array = DgElementArray<
      49             :       EvolutionMetavars,
      50             :       tmpl::flatten<tmpl::list<
      51             :           Parallel::PhaseActions<Parallel::Phase::Initialization,
      52             :                                  initialization_actions>,
      53             :           Parallel::PhaseActions<
      54             :               Parallel::Phase::RegisterWithElementDataReader,
      55             :               tmpl::list<importers::Actions::RegisterWithElementDataReader,
      56             :                          Parallel::Actions::TerminatePhase>>,
      57             :           Parallel::PhaseActions<
      58             :               Parallel::Phase::ImportInitialData,
      59             :               tmpl::list<
      60             :                   gh::Actions::SetInitialData,
      61             :                   tmpl::conditional_t<VolumeDim == 3,
      62             :                                       gh::Actions::ReceiveNumericInitialData,
      63             :                                       tmpl::list<>>,
      64             :                   Parallel::Actions::TerminatePhase>>,
      65             :           Parallel::PhaseActions<
      66             :               Parallel::Phase::InitializeInitialDataDependentQuantities,
      67             :               initialize_initial_data_dependent_quantities_actions>,
      68             :           Parallel::PhaseActions<
      69             :               Parallel::Phase::InitializeTimeStepperHistory,
      70             :               SelfStart::self_start_procedure<step_actions, system>>,
      71             :           Parallel::PhaseActions<Parallel::Phase::Register,
      72             :                                  tmpl::list<dg_registration_list,
      73             :                                             Parallel::Actions::TerminatePhase>>,
      74             :           Parallel::PhaseActions<Parallel::Phase::Restart,
      75             :                                  tmpl::list<dg_registration_list,
      76             :                                             Parallel::Actions::TerminatePhase>>,
      77             :           Parallel::PhaseActions<
      78             :               Parallel::Phase::WriteCheckpoint,
      79             :               tmpl::list<evolution::Actions::RunEventsAndTriggers<
      80             :                              Triggers::WhenToCheck::AtCheckpoints>,
      81             :                          Parallel::Actions::TerminatePhase>>,
      82             :           Parallel::PhaseActions<Parallel::Phase::CheckDomain,
      83             :                                  tmpl::list<::amr::Actions::SendAmrDiagnostics,
      84             :                                             Parallel::Actions::TerminatePhase>>,
      85             :           Parallel::PhaseActions<
      86             :               Parallel::Phase::Evolve,
      87             :               tmpl::flatten<
      88             :                   tmpl::list<evolution::Actions::RunEventsAndTriggers<
      89             :                                  Triggers::WhenToCheck::AtSteps>,
      90             :                              ::evolution::Actions::RunEventsAndTriggers<
      91             :                                  Triggers::WhenToCheck::AtSlabs>,
      92             :                              Actions::ChangeSlabSize,
      93             :                              evolution::dg::Actions::ChangeFixedLtsRatio,
      94             :                              step_actions, Actions::MutateApply<AdvanceTime<>>,
      95             :                              PhaseControl::Actions::ExecutePhaseChange>>>>>>;
      96             : 
      97             :   struct amr : tt::ConformsTo<::amr::protocols::AmrMetavariables> {
      98           0 :     using element_array = gh_dg_element_array;
      99           0 :     using projectors = tmpl::list<
     100             :         Initialization::ProjectTimeStepping<volume_dim>,
     101             :         evolution::dg::Initialization::ProjectDomain<volume_dim>,
     102             :         ::amr::projectors::ProjectVariables<volume_dim,
     103             :                                             typename system::variables_tag>,
     104             :         evolution::dg::Initialization::ProjectMortars<volume_dim>,
     105             :         Initialization::ProjectTimeStepperHistory<EvolutionMetavars>,
     106             :         evolution::Actions::ProjectRunEventsAndDenseTriggers,
     107             :         evolution::dg::Initialization::ProjectSpectralFilters<
     108             :             volume_dim, typename system::variables_tag::tags_list>,
     109             :         ::amr::projectors::DefaultInitialize<
     110             :             Initialization::Tags::InitialTimeDelta,
     111             :             Initialization::Tags::InitialSlabSize,
     112             :             ::domain::Tags::InitialExtents<volume_dim>,
     113             :             ::domain::Tags::InitialRefinementLevels<volume_dim>,
     114             :             evolution::dg::Tags::Quadrature,
     115             :             Tags::StepperErrors<typename system::variables_tag>,
     116             :             SelfStart::Tags::InitialValue<typename system::variables_tag>,
     117             :             SelfStart::Tags::InitialValue<Tags::TimeStep>>,
     118             :         ::amr::projectors::CopyFromCreatorOrLeaveAsIs<
     119             :             evolution::dg::Tags::ChangeFixedLtsRatio::NumberOfExpectedMessages,
     120             :             evolution::dg::Tags::ChangeFixedLtsRatio::NewStepSize,
     121             :             Tags::FixedLtsRatio,
     122             :             Parallel::Tags::Section<gh_dg_element_array,
     123             :                                     evolution::dg::Tags::EqualRateRegionId>,
     124             :             Tags::ChangeSlabSize::NumberOfExpectedMessages,
     125             :             Tags::ChangeSlabSize::NewSlabSize>>;
     126             :     static constexpr bool keep_coarse_grids = false;
     127             :     static constexpr bool p_refine_only_in_event = true;
     128             :   };
     129             : 
     130             :   struct registration
     131             :       : tt::ConformsTo<Parallel::protocols::RegistrationMetavariables> {
     132           0 :     using element_registrars =
     133             :         tmpl::map<tmpl::pair<gh_dg_element_array, dg_registration_list>>;
     134             :   };
     135             : 
     136           0 :   using component_list =
     137             :       tmpl::flatten<tmpl::list<::amr::Component<EvolutionMetavars>,
     138             :                                observers::Observer<EvolutionMetavars>,
     139             :                                observers::ObserverWriter<EvolutionMetavars>,
     140             :                                mem_monitor::MemoryMonitor<EvolutionMetavars>,
     141             :                                importers::ElementDataReader<EvolutionMetavars>,
     142             :                                gh_dg_element_array>>;
     143             : 
     144             :   static constexpr Options::String help{
     145             :       "Evolve the Einstein field equations using the Generalized Harmonic "
     146             :       "formulation\n"};
     147             : };