SpECTRE Documentation Coverage Report
Current view: top level - Evolution/Executables/NewtonianEuler - EvolveNewtonianEuler.hpp Hit Total Coverage
Commit: ef8d89750d7b3d281ad4d2b7c8a0b5c5e0ddc774 Lines: 0 30 0.0 %
Date: 2026-07-14 17:10:18
Legend: Lines: hit not hit

          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             : #include <cstdint>
       8             : #include <vector>
       9             : 
      10             : #include "Domain/Creators/Factory1D.hpp"
      11             : #include "Domain/Creators/Factory2D.hpp"
      12             : #include "Domain/Creators/Factory3D.hpp"
      13             : #include "Domain/Tags.hpp"
      14             : #include "Evolution/Actions/RunEventsAndDenseTriggers.hpp"
      15             : #include "Evolution/Actions/RunEventsAndTriggers.hpp"
      16             : #include "Evolution/BoundaryCorrection.hpp"
      17             : #include "Evolution/ComputeTags.hpp"
      18             : #include "Evolution/Conservative/UpdateConservatives.hpp"
      19             : #include "Evolution/DgSubcell/Actions/Initialize.hpp"
      20             : #include "Evolution/DgSubcell/Actions/Labels.hpp"
      21             : #include "Evolution/DgSubcell/Actions/ReconstructionCommunication.hpp"
      22             : #include "Evolution/DgSubcell/Actions/SelectNumericalMethod.hpp"
      23             : #include "Evolution/DgSubcell/Actions/TakeTimeStep.hpp"
      24             : #include "Evolution/DgSubcell/Actions/TciAndRollback.hpp"
      25             : #include "Evolution/DgSubcell/Actions/TciAndSwitchToDg.hpp"
      26             : #include "Evolution/DgSubcell/GetTciDecision.hpp"
      27             : #include "Evolution/DgSubcell/NeighborReconstructedFaceSolution.hpp"
      28             : #include "Evolution/DgSubcell/NeighborTciDecision.hpp"
      29             : #include "Evolution/DgSubcell/PrepareNeighborData.hpp"
      30             : #include "Evolution/DgSubcell/SubcellEqualRateRegion.hpp"
      31             : #include "Evolution/DgSubcell/Tags/ObserverCoordinates.hpp"
      32             : #include "Evolution/DgSubcell/Tags/ObserverMesh.hpp"
      33             : #include "Evolution/DgSubcell/Tags/TciStatus.hpp"
      34             : #include "Evolution/DiscontinuousGalerkin/Actions/ApplyBoundaryCorrections.hpp"
      35             : #include "Evolution/DiscontinuousGalerkin/Actions/ComputeTimeDerivative.hpp"
      36             : #include "Evolution/DiscontinuousGalerkin/CleanMortarHistory.hpp"
      37             : #include "Evolution/DiscontinuousGalerkin/DgElementArray.hpp"
      38             : #include "Evolution/DiscontinuousGalerkin/EqualRateLts/ChangeFixedLtsRatio.hpp"
      39             : #include "Evolution/DiscontinuousGalerkin/EqualRateLts/FixedLtsRatio.hpp"
      40             : #include "Evolution/DiscontinuousGalerkin/EqualRateLts/NonconformingEqualRateRegions.hpp"
      41             : #include "Evolution/DiscontinuousGalerkin/Initialization/Mortars.hpp"
      42             : #include "Evolution/DiscontinuousGalerkin/Initialization/QuadratureTag.hpp"
      43             : #include "Evolution/DiscontinuousGalerkin/Initialization/SetupEqualRateRegions.hpp"
      44             : #include "Evolution/DiscontinuousGalerkin/Initialization/SpectralFilters.hpp"
      45             : #include "Evolution/Initialization/ConservativeSystem.hpp"
      46             : #include "Evolution/Initialization/DgDomain.hpp"
      47             : #include "Evolution/Initialization/Evolution.hpp"
      48             : #include "Evolution/Initialization/SetVariables.hpp"
      49             : #include "Evolution/Systems/NewtonianEuler/AllSolutions.hpp"
      50             : #include "Evolution/Systems/NewtonianEuler/BoundaryConditions/Factory.hpp"
      51             : #include "Evolution/Systems/NewtonianEuler/BoundaryCorrections/Factory.hpp"
      52             : #include "Evolution/Systems/NewtonianEuler/FiniteDifference/Factory.hpp"
      53             : #include "Evolution/Systems/NewtonianEuler/FiniteDifference/Tag.hpp"
      54             : #include "Evolution/Systems/NewtonianEuler/SoundSpeedSquared.hpp"
      55             : #include "Evolution/Systems/NewtonianEuler/Sources/Factory.hpp"
      56             : #include "Evolution/Systems/NewtonianEuler/Subcell/NeighborPackagedData.hpp"
      57             : #include "Evolution/Systems/NewtonianEuler/Subcell/PrimitiveGhostData.hpp"
      58             : #include "Evolution/Systems/NewtonianEuler/Subcell/PrimsAfterRollback.hpp"
      59             : #include "Evolution/Systems/NewtonianEuler/Subcell/ResizeAndComputePrimitives.hpp"
      60             : #include "Evolution/Systems/NewtonianEuler/Subcell/SetInitialRdmpData.hpp"
      61             : #include "Evolution/Systems/NewtonianEuler/Subcell/TciOnDgGrid.hpp"
      62             : #include "Evolution/Systems/NewtonianEuler/Subcell/TciOnFdGrid.hpp"
      63             : #include "Evolution/Systems/NewtonianEuler/Subcell/TimeDerivative.hpp"
      64             : #include "Evolution/Systems/NewtonianEuler/System.hpp"
      65             : #include "Evolution/Systems/NewtonianEuler/Tags.hpp"
      66             : #include "IO/Observer/Actions/RegisterEvents.hpp"
      67             : #include "IO/Observer/Helpers.hpp"
      68             : #include "IO/Observer/ObserverComponent.hpp"
      69             : #include "NumericalAlgorithms/DiscontinuousGalerkin/Formulation.hpp"
      70             : #include "NumericalAlgorithms/DiscontinuousGalerkin/Tags.hpp"
      71             : #include "NumericalAlgorithms/LinearOperators/Filters/Factory.hpp"
      72             : #include "Options/Protocols/FactoryCreation.hpp"
      73             : #include "Options/String.hpp"
      74             : #include "Parallel/Local.hpp"
      75             : #include "Parallel/Phase.hpp"
      76             : #include "Parallel/PhaseControl/CheckpointAndExitAfterWallclock.hpp"
      77             : #include "Parallel/PhaseControl/ExecutePhaseChange.hpp"
      78             : #include "Parallel/PhaseControl/Factory.hpp"
      79             : #include "Parallel/PhaseControl/VisitAndReturn.hpp"
      80             : #include "Parallel/PhaseDependentActionList.hpp"
      81             : #include "Parallel/Protocols/RegistrationMetavariables.hpp"
      82             : #include "ParallelAlgorithms/Actions/AddComputeTags.hpp"
      83             : #include "ParallelAlgorithms/Actions/InitializeItems.hpp"
      84             : #include "ParallelAlgorithms/Actions/MutateApply.hpp"
      85             : #include "ParallelAlgorithms/Actions/SpectralFilter.hpp"
      86             : #include "ParallelAlgorithms/Actions/TerminatePhase.hpp"
      87             : #include "ParallelAlgorithms/Events/ChangeFixedLtsRatio.hpp"
      88             : #include "ParallelAlgorithms/Events/Completion.hpp"
      89             : #include "ParallelAlgorithms/Events/Factory.hpp"
      90             : #include "ParallelAlgorithms/Events/Tags.hpp"
      91             : #include "ParallelAlgorithms/EventsAndDenseTriggers/DenseTrigger.hpp"
      92             : #include "ParallelAlgorithms/EventsAndDenseTriggers/DenseTriggers/Factory.hpp"
      93             : #include "ParallelAlgorithms/EventsAndTriggers/Event.hpp"
      94             : #include "ParallelAlgorithms/EventsAndTriggers/EventsAndTriggers.hpp"
      95             : #include "ParallelAlgorithms/EventsAndTriggers/LogicalTriggers.hpp"
      96             : #include "ParallelAlgorithms/EventsAndTriggers/Trigger.hpp"
      97             : #include "PointwiseFunctions/AnalyticData/AnalyticData.hpp"
      98             : #include "PointwiseFunctions/AnalyticData/NewtonianEuler/KhInstability.hpp"
      99             : #include "PointwiseFunctions/AnalyticSolutions/AnalyticSolution.hpp"
     100             : #include "PointwiseFunctions/AnalyticSolutions/NewtonianEuler/IsentropicVortex.hpp"
     101             : #include "PointwiseFunctions/AnalyticSolutions/NewtonianEuler/LaneEmdenStar.hpp"
     102             : #include "PointwiseFunctions/AnalyticSolutions/NewtonianEuler/RiemannProblem.hpp"
     103             : #include "PointwiseFunctions/AnalyticSolutions/NewtonianEuler/SmoothFlow.hpp"
     104             : #include "PointwiseFunctions/AnalyticSolutions/Tags.hpp"
     105             : #include "PointwiseFunctions/Hydro/EquationsOfState/Factory.hpp"
     106             : #include "PointwiseFunctions/Hydro/EquationsOfState/RegisterDerivedWithCharm.hpp"
     107             : #include "PointwiseFunctions/Hydro/Tags.hpp"
     108             : #include "Time/Actions/SelfStartActions.hpp"
     109             : #include "Time/AdvanceTime.hpp"
     110             : #include "Time/ChangeSlabSize/Action.hpp"
     111             : #include "Time/ChangeStepSize.hpp"
     112             : #include "Time/ChangeTimeStepperOrder.hpp"
     113             : #include "Time/CleanHistory.hpp"
     114             : #include "Time/RecordTimeStepperData.hpp"
     115             : #include "Time/StepChoosers/Factory.hpp"
     116             : #include "Time/StepChoosers/StepChooser.hpp"
     117             : #include "Time/Tags/Time.hpp"
     118             : #include "Time/Tags/TimeStepId.hpp"
     119             : #include "Time/TimeSequence.hpp"
     120             : #include "Time/TimeSteppers/Factory.hpp"
     121             : #include "Time/TimeSteppers/LtsTimeStepper.hpp"
     122             : #include "Time/TimeSteppers/TimeStepper.hpp"
     123             : #include "Time/Triggers/TimeTriggers.hpp"
     124             : #include "Time/UpdateU.hpp"
     125             : #include "Utilities/Functional.hpp"
     126             : #include "Utilities/ProtocolHelpers.hpp"
     127             : #include "Utilities/TMPL.hpp"
     128             : 
     129             : /// \cond
     130             : namespace PUP {
     131             : class er;
     132             : }  // namespace PUP
     133             : namespace Parallel {
     134             : template <typename Metavariables>
     135             : class CProxy_GlobalCache;
     136             : }  // namespace Parallel
     137             : /// \endcond
     138             : 
     139             : template <size_t Dim>
     140             : struct EvolutionMetavars {
     141             :   static constexpr size_t volume_dim = Dim;
     142             :   // The use_dg_subcell flag controls whether to use unlimited DG (false)
     143             :   // or a DG-FD hybrid scheme (true).
     144             :   static constexpr bool use_dg_subcell = true;
     145             : 
     146           0 :   using system = NewtonianEuler::System<Dim>;
     147             : 
     148           0 :   using temporal_id = Tags::TimeStepId;
     149           0 :   using TimeStepperBase = TimeStepper;
     150             : 
     151             :   static constexpr bool local_time_stepping =
     152             :       TimeStepperBase::local_time_stepping;
     153             :   static constexpr bool use_dg_element_collection = false;
     154             : 
     155           0 :   using initial_data_tag = evolution::initial_data::Tags::InitialData;
     156           0 :   using initial_data_list = NewtonianEuler::InitialData::initial_data_list<Dim>;
     157             : 
     158           0 :   using analytic_variables_tags =
     159             :       typename system::primitive_variables_tag::tags_list;
     160             : 
     161           0 :   using equation_of_state_tag = hydro::Tags::EquationOfState<false, 2>;
     162             : 
     163           0 :   using analytic_compute = evolution::Tags::AnalyticSolutionsCompute<
     164             :       volume_dim, analytic_variables_tags, use_dg_subcell, initial_data_list>;
     165           0 :   using error_compute = Tags::ErrorsCompute<analytic_variables_tags>;
     166           0 :   using error_tags = db::wrap_tags_in<Tags::Error, analytic_variables_tags>;
     167           0 :   using observe_fields = tmpl::push_back<
     168             :       tmpl::append<
     169             :           typename system::variables_tag::tags_list,
     170             :           typename system::primitive_variables_tag::tags_list, error_tags,
     171             :           tmpl::conditional_t<use_dg_subcell,
     172             :                               tmpl::list<evolution::dg::subcell::Tags::
     173             :                                              TciStatusCompute<volume_dim>>,
     174             :                               tmpl::list<>>>,
     175             :       tmpl::conditional_t<
     176             :           use_dg_subcell,
     177             :           evolution::dg::subcell::Tags::ObserverCoordinatesCompute<
     178             :               volume_dim, Frame::ElementLogical>,
     179             :           ::Events::Tags::ObserverCoordinatesCompute<volume_dim,
     180             :                                                      Frame::ElementLogical>>,
     181             :       tmpl::conditional_t<
     182             :           use_dg_subcell,
     183             :           evolution::dg::subcell::Tags::ObserverCoordinatesCompute<volume_dim,
     184             :                                                                    Frame::Grid>,
     185             :           domain::Tags::Coordinates<volume_dim, Frame::Grid>>,
     186             :       tmpl::conditional_t<
     187             :           use_dg_subcell,
     188             :           evolution::dg::subcell::Tags::ObserverCoordinatesCompute<
     189             :               volume_dim, Frame::Inertial>,
     190             :           domain::Tags::Coordinates<volume_dim, Frame::Inertial>>>;
     191           0 :   using non_tensor_compute_tags = tmpl::append<
     192             :       tmpl::conditional_t<
     193             :           use_dg_subcell,
     194             :           tmpl::list<
     195             :               evolution::dg::subcell::Tags::ObserverMeshCompute<volume_dim>,
     196             :               evolution::dg::subcell::Tags::ObserverInverseJacobianCompute<
     197             :                   volume_dim, Frame::ElementLogical, Frame::Inertial>,
     198             :               evolution::dg::subcell::Tags::
     199             :                   ObserverJacobianAndDetInvJacobianCompute<
     200             :                       volume_dim, Frame::ElementLogical, Frame::Inertial>>,
     201             :           tmpl::list<::Events::Tags::ObserverMeshCompute<volume_dim>,
     202             :                      ::Events::Tags::ObserverInverseJacobianCompute<
     203             :                          volume_dim, Frame::ElementLogical, Frame::Inertial>,
     204             :                      ::Events::Tags::ObserverJacobianCompute<
     205             :                          volume_dim, Frame::ElementLogical, Frame::Inertial>,
     206             :                      ::Events::Tags::ObserverDetInvJacobianCompute<
     207             :                          Frame::ElementLogical, Frame::Inertial>>>,
     208             :       tmpl::list<analytic_compute, error_compute>>;
     209             : 
     210             :   struct factory_creation
     211             :       : tt::ConformsTo<Options::protocols::FactoryCreation> {
     212           0 :     using factory_classes = tmpl::map<
     213             :         tmpl::pair<DenseTrigger, DenseTriggers::standard_dense_triggers>,
     214             :         tmpl::pair<DomainCreator<volume_dim>, domain_creators<volume_dim>>,
     215             :         tmpl::pair<NewtonianEuler::Sources::Source<Dim>,
     216             :                    NewtonianEuler::Sources::all_sources<Dim>>,
     217             :         tmpl::pair<evolution::initial_data::InitialData,
     218             :                    NewtonianEuler::InitialData::initial_data_list<Dim>>,
     219             :         tmpl::pair<
     220             :             Event,
     221             :             tmpl::flatten<tmpl::list<
     222             :                 Events::Completion,
     223             :                 dg::Events::field_observations<volume_dim, observe_fields,
     224             :                                                non_tensor_compute_tags>,
     225             :                 Events::time_events<system>,
     226             :                 tmpl::conditional_t<local_time_stepping,
     227             :                                     dg::Events::ChangeFixedLtsRatio<volume_dim>,
     228             :                                     tmpl::list<>>>>>,
     229             :         tmpl::pair<evolution::BoundaryCorrection,
     230             :                    NewtonianEuler::BoundaryCorrections::
     231             :                        standard_boundary_corrections<volume_dim>>,
     232             :         tmpl::pair<LtsTimeStepper, TimeSteppers::lts_time_steppers>,
     233             :         tmpl::pair<
     234             :             NewtonianEuler::BoundaryConditions::BoundaryCondition<volume_dim>,
     235             :             NewtonianEuler::BoundaryConditions::standard_boundary_conditions<
     236             :                 volume_dim>>,
     237             :         tmpl::pair<PhaseChange, PhaseControl::factory_creatable_classes>,
     238             :         tmpl::pair<StepChooser<StepChooserUse::LtsStep>,
     239             :                    StepChoosers::standard_step_choosers<system>>,
     240             :         tmpl::pair<StepChooser<StepChooserUse::Slab>,
     241             :                    tmpl::append<StepChoosers::standard_slab_choosers<system>,
     242             :                                 tmpl::conditional_t<
     243             :                                     local_time_stepping,
     244             :                                     tmpl::list<evolution::dg::StepChoosers::
     245             :                                                    FixedLtsRatio<volume_dim>>,
     246             :                                     tmpl::list<>>>>,
     247             :         tmpl::pair<TimeSequence<double>,
     248             :                    TimeSequences::all_time_sequences<double>>,
     249             :         tmpl::pair<TimeSequence<std::uint64_t>,
     250             :                    TimeSequences::all_time_sequences<std::uint64_t>>,
     251             :         tmpl::pair<TimeStepper, TimeSteppers::time_steppers>,
     252             :         tmpl::pair<Trigger, tmpl::append<Triggers::logical_triggers,
     253             :                                          Triggers::time_triggers>>,
     254             :         tmpl::pair<Filters::Filter<volume_dim,
     255             :                                    typename system::variables_tag::tags_list>,
     256             :                    Filters::all_filters<
     257             :                        volume_dim, typename system::variables_tag::tags_list>>>;
     258             :   };
     259             : 
     260           0 :   using observed_reduction_data_tags =
     261             :       observers::collect_reduction_data_tags<tmpl::flatten<tmpl::list<
     262             :           tmpl::at<typename factory_creation::factory_classes, Event>>>>;
     263             : 
     264           0 :   using dg_registration_list =
     265             :       tmpl::list<observers::Actions::RegisterEventsWithObservers>;
     266             : 
     267           0 :   using equal_rate_regions = tmpl::flatten<
     268             :       tmpl::list<evolution::dg::NonconformingEqualRateRegions<volume_dim>,
     269             :                  tmpl::conditional_t<
     270             :                      use_dg_subcell,
     271             :                      evolution::dg::subcell::SubcellEqualRateRegion<volume_dim>,
     272             :                      tmpl::list<>>>>;
     273             : 
     274           0 :   using initialization_actions = tmpl::flatten<tmpl::list<
     275             :       Initialization::Actions::InitializeItems<
     276             :           Initialization::TimeStepping<EvolutionMetavars, TimeStepperBase,
     277             :                                        false>,
     278             :           evolution::dg::Initialization::Domain<EvolutionMetavars>,
     279             :           Initialization::TimeStepperHistory<EvolutionMetavars>>,
     280             :       Initialization::Actions::ConservativeSystem<system>,
     281             : 
     282             :       tmpl::conditional_t<
     283             :           use_dg_subcell,
     284             :           tmpl::list<
     285             :               evolution::dg::subcell::Actions::SetSubcellGrid<volume_dim,
     286             :                                                               system, false>,
     287             :               Actions::UpdateConservatives,
     288             :               evolution::dg::subcell::Actions::SetAndCommunicateInitialRdmpData<
     289             :                   volume_dim,
     290             :                   NewtonianEuler::subcell::SetInitialRdmpData<volume_dim>>,
     291             :               evolution::dg::subcell::Actions::ComputeAndSendTciOnInitialGrid<
     292             :                   volume_dim, system,
     293             :                   NewtonianEuler::subcell::TciOnFdGrid<volume_dim>>,
     294             :               evolution::dg::subcell::Actions::SetInitialGridFromTciData<
     295             :                   volume_dim, system>,
     296             :               Actions::MutateApply<
     297             :                   NewtonianEuler::subcell::ResizeAndComputePrims<volume_dim>>,
     298             :               Actions::UpdateConservatives>,
     299             :           tmpl::list<evolution::Initialization::Actions::SetVariables<
     300             :                          domain::Tags::Coordinates<Dim, Frame::ElementLogical>>,
     301             :                      Actions::UpdateConservatives>>,
     302             :       Initialization::Actions::AddComputeTags<
     303             :           tmpl::list<NewtonianEuler::Tags::SoundSpeedSquaredCompute<DataVector>,
     304             :                      NewtonianEuler::Tags::SoundSpeedCompute<DataVector>>>,
     305             :       Initialization::Actions::AddComputeTags<
     306             :           StepChoosers::step_chooser_compute_tags<EvolutionMetavars>>,
     307             :       ::evolution::dg::Initialization::Mortars<volume_dim>,
     308             :       evolution::dg::Initialization::Actions::SetupEqualRateRegions<
     309             :           EvolutionMetavars, volume_dim, equal_rate_regions>,
     310             :       evolution::Actions::InitializeRunEventsAndDenseTriggers,
     311             :       Initialization::Actions::InitializeItems<
     312             :           evolution::dg::Initialization::SpectralFilters<
     313             :               volume_dim, typename system::variables_tag::tags_list>>,
     314             :       Parallel::Actions::TerminatePhase>>;
     315             : 
     316           0 :   using events_and_dense_triggers_postprocessors = tmpl::list<
     317             :       AlwaysReadyPostprocessor<typename system::primitive_from_conservative>>;
     318             : 
     319           0 :   using dg_step_actions = tmpl::flatten<tmpl::list<
     320             :       evolution::dg::Actions::ComputeTimeDerivative<
     321             :           volume_dim, system, AllStepChoosers, local_time_stepping,
     322             :           use_dg_element_collection>,
     323             :       evolution::dg::Actions::ApplyBoundaryCorrectionsToTimeDerivative<
     324             :           volume_dim, use_dg_element_collection>,
     325             :       Actions::MutateApply<RecordTimeStepperData<system>>,
     326             :       evolution::Actions::RunEventsAndDenseTriggers<tmpl::push_front<
     327             :           events_and_dense_triggers_postprocessors,
     328             :           evolution::dg::ApplyLtsDenseBoundaryCorrections<EvolutionMetavars>>>,
     329             :       Actions::MutateApply<UpdateU<system, local_time_stepping>>,
     330             :       evolution::dg::Actions::ApplyLtsBoundaryCorrections<
     331             :           volume_dim, use_dg_element_collection>,
     332             :       tmpl::conditional_t<
     333             :           local_time_stepping,
     334             :           tmpl::list<Actions::MutateApply<ChangeTimeStepperOrder<system>>>,
     335             :           tmpl::list<>>,
     336             :       tmpl::conditional_t<
     337             :           use_dg_subcell,
     338             :           // Note: The primitive variables are computed as part of the TCI.
     339             :           evolution::dg::subcell::Actions::TciAndRollback<
     340             :               NewtonianEuler::subcell::TciOnDgGrid<volume_dim>>,
     341             :           tmpl::list<
     342             :               // Conservative `UpdatePrimitives` expects system to
     343             :               // possess list of recovery schemes so we use
     344             :               // `MutateApply` instead.
     345             :               Actions::MutateApply<
     346             :                   typename system::primitive_from_conservative>>>,
     347             :       Actions::MutateApply<CleanHistory<system>>,
     348             :       Actions::MutateApply<evolution::dg::CleanMortarHistory<volume_dim>>,
     349             :       dg::Actions::SpectralFilter>>;
     350             : 
     351             :   struct SubcellOptions {
     352             :     static constexpr bool subcell_enabled = use_dg_subcell;
     353             :     static constexpr bool subcell_enabled_at_external_boundary = false;
     354             : 
     355             :     // We send `ghost_zone_size` cell-centered grid points for variable
     356             :     // reconstruction, of which we need `ghost_zone_size-1` for reconstruction
     357             :     // to the internal side of the element face, and `ghost_zone_size` for
     358             :     // reconstruction to the external side of the element face.
     359             :     template <typename DbTagsList>
     360           0 :     static constexpr size_t ghost_zone_size(
     361             :         const db::DataBox<DbTagsList>& box) {
     362             :       return db::get<NewtonianEuler::fd::Tags::Reconstructor<Dim>>(box)
     363             :           .ghost_zone_size();
     364             :     }
     365             : 
     366           0 :     using DgComputeSubcellNeighborPackagedData =
     367             :         NewtonianEuler::subcell::NeighborPackagedData;
     368             : 
     369           0 :     using GhostVariables =
     370             :         NewtonianEuler::subcell::PrimitiveGhostVariables<volume_dim>;
     371             :   };
     372             : 
     373           0 :   using dg_subcell_step_actions = tmpl::flatten<tmpl::list<
     374             :       evolution::dg::subcell::Actions::SelectNumericalMethod,
     375             : 
     376             :       Actions::Label<evolution::dg::subcell::Actions::Labels::BeginDg>,
     377             :       dg_step_actions,
     378             :       Actions::Goto<evolution::dg::subcell::Actions::Labels::EndOfSolvers>,
     379             : 
     380             :       Actions::Label<evolution::dg::subcell::Actions::Labels::BeginSubcell>,
     381             :       tmpl::conditional_t<local_time_stepping,
     382             :                           // This is just to adjust for FixedLtsRatio, so we
     383             :                           // can pass an empty list of StepChoosers.
     384             :                           Actions::MutateApply<ChangeStepSize<tmpl::list<>>>,
     385             :                           tmpl::list<>>,
     386             :       evolution::dg::subcell::Actions::SendDataForReconstruction<
     387             :           volume_dim,
     388             :           NewtonianEuler::subcell::PrimitiveGhostVariables<volume_dim>,
     389             :           use_dg_element_collection>,
     390             :       evolution::dg::subcell::Actions::ReceiveDataForReconstruction<volume_dim>,
     391             :       Actions::Label<
     392             :           evolution::dg::subcell::Actions::Labels::BeginSubcellAfterDgRollback>,
     393             :       Actions::MutateApply<
     394             :           NewtonianEuler::subcell::PrimsAfterRollback<volume_dim>>,
     395             :       evolution::dg::subcell::fd::Actions::TakeTimeStep<
     396             :           NewtonianEuler::subcell::TimeDerivative<volume_dim>>,
     397             :       Actions::MutateApply<RecordTimeStepperData<system>>,
     398             :       evolution::Actions::RunEventsAndDenseTriggers<
     399             :           events_and_dense_triggers_postprocessors>,
     400             :       Actions::MutateApply<UpdateU<system, local_time_stepping>>,
     401             :       Actions::MutateApply<CleanHistory<system>>,
     402             :       Actions::MutateApply<evolution::dg::CleanMortarHistory<volume_dim>>,
     403             :       Actions::MutateApply<typename system::primitive_from_conservative>,
     404             :       evolution::dg::subcell::Actions::TciAndSwitchToDg<
     405             :           NewtonianEuler::subcell::TciOnFdGrid<volume_dim>>,
     406             :       Actions::MutateApply<
     407             :           NewtonianEuler::subcell::ResizeAndComputePrims<volume_dim>>,
     408             : 
     409             :       Actions::Label<evolution::dg::subcell::Actions::Labels::EndOfSolvers>>>;
     410             : 
     411           0 :   using step_actions =
     412             :       tmpl::conditional_t<use_dg_subcell, dg_subcell_step_actions,
     413             :                           dg_step_actions>;
     414             : 
     415           0 :   using dg_element_array = DgElementArray<
     416             :       EvolutionMetavars,
     417             :       tmpl::list<
     418             :           Parallel::PhaseActions<Parallel::Phase::Initialization,
     419             :                                  initialization_actions>,
     420             : 
     421             :           Parallel::PhaseActions<
     422             :               Parallel::Phase::InitializeTimeStepperHistory,
     423             :               SelfStart::self_start_procedure<step_actions, system>>,
     424             : 
     425             :           Parallel::PhaseActions<Parallel::Phase::Register,
     426             :                                  tmpl::list<dg_registration_list,
     427             :                                             Parallel::Actions::TerminatePhase>>,
     428             : 
     429             :           Parallel::PhaseActions<Parallel::Phase::Restart,
     430             :                                  tmpl::list<dg_registration_list,
     431             :                                             Parallel::Actions::TerminatePhase>>,
     432             : 
     433             :           Parallel::PhaseActions<
     434             :               Parallel::Phase::WriteCheckpoint,
     435             :               tmpl::list<evolution::Actions::RunEventsAndTriggers<
     436             :                              Triggers::WhenToCheck::AtCheckpoints>,
     437             :                          Parallel::Actions::TerminatePhase>>,
     438             : 
     439             :           Parallel::PhaseActions<
     440             :               Parallel::Phase::Evolve,
     441             :               tmpl::flatten<tmpl::list<
     442             :                   std::conditional_t<local_time_stepping,
     443             :                                      evolution::Actions::RunEventsAndTriggers<
     444             :                                          Triggers::WhenToCheck::AtSteps>,
     445             :                                      tmpl::list<>>,
     446             :                   evolution::Actions::RunEventsAndTriggers<
     447             :                       Triggers::WhenToCheck::AtSlabs>,
     448             :                   Actions::ChangeSlabSize,
     449             :                   std::conditional_t<
     450             :                       local_time_stepping,
     451             :                       evolution::dg::Actions::ChangeFixedLtsRatio,
     452             :                       tmpl::list<>>,
     453             :                   step_actions, Actions::MutateApply<AdvanceTime<>>,
     454             :                   PhaseControl::Actions::ExecutePhaseChange>>>>>;
     455             : 
     456             :   struct registration
     457             :       : tt::ConformsTo<Parallel::protocols::RegistrationMetavariables> {
     458           0 :     using element_registrars =
     459             :         tmpl::map<tmpl::pair<dg_element_array, dg_registration_list>>;
     460             :   };
     461             : 
     462           0 :   using component_list =
     463             :       tmpl::list<observers::Observer<EvolutionMetavars>,
     464             :                  observers::ObserverWriter<EvolutionMetavars>,
     465             :                  dg_element_array>;
     466             : 
     467           0 :   using const_global_cache_tags = tmpl::push_back<
     468             :       tmpl::conditional_t<
     469             :           use_dg_subcell,
     470             :           tmpl::list<NewtonianEuler::fd::Tags::Reconstructor<volume_dim>>,
     471             :           tmpl::list<>>,
     472             :       initial_data_tag, equation_of_state_tag,
     473             :       NewtonianEuler::Tags::SourceTerm<Dim>>;
     474             : 
     475             :   static constexpr Options::String help{
     476             :       "Evolve the Newtonian Euler system in conservative form.\n\n"};
     477             : 
     478             :   static constexpr std::array<Parallel::Phase, 5> default_phase_order{
     479             :       {Parallel::Phase::Initialization,
     480             :        Parallel::Phase::InitializeTimeStepperHistory, Parallel::Phase::Register,
     481             :        Parallel::Phase::Evolve, Parallel::Phase::Exit}};
     482             : 
     483             :   // NOLINTNEXTLINE(google-runtime-references)
     484           0 :   void pup(PUP::er& /*p*/) {}
     485             : };

Generated by: LCOV version 1.14