SpECTRE Documentation Coverage Report
Current view: top level - Evolution/Executables/GeneralizedHarmonic - GeneralizedHarmonicBase.hpp Hit Total Coverage
Commit: 41764896b49d19fe02ce091a3d0f0299986eabbe Lines: 0 16 0.0 %
Date: 2024-12-20 19:18:20
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 <vector>
       8             : 
       9             : #include "ControlSystem/Actions/LimitTimeStep.hpp"
      10             : #include "DataStructures/DataBox/PrefixHelpers.hpp"
      11             : #include "DataStructures/DataBox/Tag.hpp"
      12             : #include "DataStructures/Tensor/EagerMath/RaiseOrLowerIndex.hpp"
      13             : #include "Domain/Creators/Factory1D.hpp"
      14             : #include "Domain/Creators/Factory2D.hpp"
      15             : #include "Domain/Creators/Factory3D.hpp"
      16             : #include "Domain/Tags.hpp"
      17             : #include "Domain/TagsCharacteristicSpeeds.hpp"
      18             : #include "Evolution/Actions/RunEventsAndDenseTriggers.hpp"
      19             : #include "Evolution/ComputeTags.hpp"
      20             : #include "Evolution/DiscontinuousGalerkin/Actions/ApplyBoundaryCorrections.hpp"
      21             : #include "Evolution/DiscontinuousGalerkin/Actions/ComputeTimeDerivative.hpp"
      22             : #include "Evolution/DiscontinuousGalerkin/DgElementArray.hpp"
      23             : #include "Evolution/DiscontinuousGalerkin/Initialization/Mortars.hpp"
      24             : #include "Evolution/Initialization/DgDomain.hpp"
      25             : #include "Evolution/Initialization/Evolution.hpp"
      26             : #include "Evolution/Initialization/NonconservativeSystem.hpp"
      27             : #include "Evolution/Systems/GeneralizedHarmonic/Actions/SetInitialData.hpp"
      28             : #include "Evolution/Systems/GeneralizedHarmonic/BoundaryConditions/Factory.hpp"
      29             : #include "Evolution/Systems/GeneralizedHarmonic/BoundaryCorrections/Factory.hpp"
      30             : #include "Evolution/Systems/GeneralizedHarmonic/Equations.hpp"
      31             : #include "Evolution/Systems/GeneralizedHarmonic/GaugeSourceFunctions/Factory.hpp"
      32             : #include "Evolution/Systems/GeneralizedHarmonic/GaugeSourceFunctions/Gauges.hpp"
      33             : #include "Evolution/Systems/GeneralizedHarmonic/GaugeSourceFunctions/SetPiAndPhiFromConstraints.hpp"
      34             : #include "Evolution/Systems/GeneralizedHarmonic/GaugeSourceFunctions/Tags/GaugeCondition.hpp"
      35             : #include "Evolution/Systems/GeneralizedHarmonic/Initialize.hpp"
      36             : #include "Evolution/Systems/GeneralizedHarmonic/System.hpp"
      37             : #include "Evolution/Systems/GeneralizedHarmonic/Tags.hpp"
      38             : #include "Evolution/Tags/Filter.hpp"
      39             : #include "Evolution/TypeTraits.hpp"
      40             : #include "IO/Importers/Actions/RegisterWithElementDataReader.hpp"
      41             : #include "IO/Importers/ElementDataReader.hpp"
      42             : #include "IO/Observer/Actions/ObserverRegistration.hpp"
      43             : #include "IO/Observer/Actions/RegisterEvents.hpp"
      44             : #include "IO/Observer/Helpers.hpp"
      45             : #include "IO/Observer/ObserverComponent.hpp"
      46             : #include "IO/Observer/Tags.hpp"
      47             : #include "NumericalAlgorithms/DiscontinuousGalerkin/Tags.hpp"
      48             : #include "NumericalAlgorithms/LinearOperators/ExponentialFilter.hpp"
      49             : #include "Options/Protocols/FactoryCreation.hpp"
      50             : #include "Options/String.hpp"
      51             : #include "Parallel/Algorithms/AlgorithmSingleton.hpp"
      52             : #include "Parallel/Local.hpp"
      53             : #include "Parallel/Phase.hpp"
      54             : #include "Parallel/PhaseControl/CheckpointAndExitAfterWallclock.hpp"
      55             : #include "Parallel/PhaseControl/ExecutePhaseChange.hpp"
      56             : #include "Parallel/PhaseControl/Factory.hpp"
      57             : #include "Parallel/PhaseControl/VisitAndReturn.hpp"
      58             : #include "Parallel/PhaseDependentActionList.hpp"
      59             : #include "Parallel/Reduction.hpp"
      60             : #include "ParallelAlgorithms/Actions/AddComputeTags.hpp"
      61             : #include "ParallelAlgorithms/Actions/FilterAction.hpp"
      62             : #include "ParallelAlgorithms/Actions/FunctionsOfTimeAreReady.hpp"
      63             : #include "ParallelAlgorithms/Actions/InitializeItems.hpp"
      64             : #include "ParallelAlgorithms/Actions/MemoryMonitor/ContributeMemoryData.hpp"
      65             : #include "ParallelAlgorithms/Actions/MutateApply.hpp"
      66             : #include "ParallelAlgorithms/Actions/TerminatePhase.hpp"
      67             : #include "ParallelAlgorithms/Amr/Actions/CollectDataFromChildren.hpp"
      68             : #include "ParallelAlgorithms/Amr/Actions/Component.hpp"
      69             : #include "ParallelAlgorithms/Amr/Actions/CreateChild.hpp"
      70             : #include "ParallelAlgorithms/Amr/Actions/Initialize.hpp"
      71             : #include "ParallelAlgorithms/Amr/Actions/SendAmrDiagnostics.hpp"
      72             : #include "ParallelAlgorithms/Amr/Criteria/Constraints.hpp"
      73             : #include "ParallelAlgorithms/Amr/Criteria/Criterion.hpp"
      74             : #include "ParallelAlgorithms/Amr/Criteria/DriveToTarget.hpp"
      75             : #include "ParallelAlgorithms/Amr/Criteria/Random.hpp"
      76             : #include "ParallelAlgorithms/Amr/Criteria/Tags/Criteria.hpp"
      77             : #include "ParallelAlgorithms/Amr/Criteria/TruncationError.hpp"
      78             : #include "ParallelAlgorithms/Amr/Projectors/CopyFromCreatorOrLeaveAsIs.hpp"
      79             : #include "ParallelAlgorithms/Amr/Projectors/DefaultInitialize.hpp"
      80             : #include "ParallelAlgorithms/Amr/Projectors/Tensors.hpp"
      81             : #include "ParallelAlgorithms/Amr/Projectors/Variables.hpp"
      82             : #include "ParallelAlgorithms/Amr/Protocols/AmrMetavariables.hpp"
      83             : #include "ParallelAlgorithms/ApparentHorizonFinder/Callbacks/FindApparentHorizon.hpp"
      84             : #include "ParallelAlgorithms/ApparentHorizonFinder/InterpolationTarget.hpp"
      85             : #include "ParallelAlgorithms/Events/Factory.hpp"
      86             : #include "ParallelAlgorithms/Events/MonitorMemory.hpp"
      87             : #include "ParallelAlgorithms/Events/ObserveTimeStep.hpp"
      88             : #include "ParallelAlgorithms/Events/ObserveTimeStepVolume.hpp"
      89             : #include "ParallelAlgorithms/Events/Tags.hpp"
      90             : #include "ParallelAlgorithms/EventsAndDenseTriggers/DenseTrigger.hpp"
      91             : #include "ParallelAlgorithms/EventsAndDenseTriggers/DenseTriggers/Factory.hpp"
      92             : #include "ParallelAlgorithms/EventsAndTriggers/Completion.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 "ParallelAlgorithms/Interpolation/Actions/CleanUpInterpolator.hpp"
      98             : #include "ParallelAlgorithms/Interpolation/Actions/InitializeInterpolationTarget.hpp"
      99             : #include "ParallelAlgorithms/Interpolation/Actions/InterpolationTargetReceiveVars.hpp"
     100             : #include "ParallelAlgorithms/Interpolation/Actions/InterpolatorReceivePoints.hpp"
     101             : #include "ParallelAlgorithms/Interpolation/Actions/InterpolatorReceiveVolumeData.hpp"
     102             : #include "ParallelAlgorithms/Interpolation/Actions/InterpolatorRegisterElement.hpp"
     103             : #include "ParallelAlgorithms/Interpolation/Actions/TryToInterpolate.hpp"
     104             : #include "ParallelAlgorithms/Interpolation/Callbacks/ObserveTimeSeriesOnSurface.hpp"
     105             : #include "ParallelAlgorithms/Interpolation/Events/Interpolate.hpp"
     106             : #include "ParallelAlgorithms/Interpolation/InterpolationTarget.hpp"
     107             : #include "ParallelAlgorithms/Interpolation/Interpolator.hpp"
     108             : #include "ParallelAlgorithms/Interpolation/Tags.hpp"
     109             : #include "PointwiseFunctions/AnalyticSolutions/GeneralRelativity/Factory.hpp"
     110             : #include "PointwiseFunctions/AnalyticSolutions/GeneralRelativity/KerrSchild.hpp"
     111             : #include "PointwiseFunctions/AnalyticSolutions/GeneralRelativity/SphericalKerrSchild.hpp"
     112             : #include "PointwiseFunctions/AnalyticSolutions/GeneralRelativity/WrappedGr.hpp"
     113             : #include "PointwiseFunctions/AnalyticSolutions/Tags.hpp"
     114             : #include "PointwiseFunctions/GeneralRelativity/Christoffel.hpp"
     115             : #include "PointwiseFunctions/GeneralRelativity/DetAndInverseSpatialMetric.hpp"
     116             : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/ConstraintGammas.hpp"
     117             : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/ExtrinsicCurvature.hpp"
     118             : #include "PointwiseFunctions/GeneralRelativity/Psi4Real.hpp"
     119             : #include "PointwiseFunctions/GeneralRelativity/Ricci.hpp"
     120             : #include "PointwiseFunctions/GeneralRelativity/SpacetimeNormalVector.hpp"
     121             : #include "PointwiseFunctions/GeneralRelativity/Surfaces/Tags.hpp"
     122             : #include "PointwiseFunctions/GeneralRelativity/Tags.hpp"
     123             : #include "PointwiseFunctions/GeneralRelativity/WeylElectric.hpp"
     124             : #include "PointwiseFunctions/GeneralRelativity/WeylTypeD1.hpp"
     125             : #include "PointwiseFunctions/InitialDataUtilities/InitialData.hpp"
     126             : #include "PointwiseFunctions/InitialDataUtilities/Tags/InitialData.hpp"
     127             : #include "PointwiseFunctions/MathFunctions/Factory.hpp"
     128             : #include "PointwiseFunctions/MathFunctions/MathFunction.hpp"
     129             : #include "Time/Actions/AdvanceTime.hpp"
     130             : #include "Time/Actions/CleanHistory.hpp"
     131             : #include "Time/Actions/RecordTimeStepperData.hpp"
     132             : #include "Time/Actions/UpdateU.hpp"
     133             : #include "Time/StepChoosers/Factory.hpp"
     134             : #include "Time/StepChoosers/StepChooser.hpp"
     135             : #include "Time/Tags/Time.hpp"
     136             : #include "Time/TimeSequence.hpp"
     137             : #include "Time/TimeSteppers/Factory.hpp"
     138             : #include "Time/TimeSteppers/LtsTimeStepper.hpp"
     139             : #include "Time/TimeSteppers/TimeStepper.hpp"
     140             : #include "Time/Triggers/TimeTriggers.hpp"
     141             : #include "Utilities/ErrorHandling/Error.hpp"
     142             : #include "Utilities/Functional.hpp"
     143             : #include "Utilities/ProtocolHelpers.hpp"
     144             : #include "Utilities/TMPL.hpp"
     145             : 
     146             : /// \cond
     147             : namespace Frame {
     148             : 
     149             : struct Inertial;
     150             : }  // namespace Frame
     151             : namespace PUP {
     152             : class er;
     153             : }  // namespace PUP
     154             : namespace Parallel {
     155             : template <typename Metavariables>
     156             : class CProxy_GlobalCache;
     157             : }  // namespace Parallel
     158             : /// \endcond
     159             : 
     160             : namespace detail {
     161             : 
     162             : template <size_t volume_dim>
     163             : struct ObserverTags {
     164             :   using system = gh::System<volume_dim>;
     165             : 
     166             :   using variables_tag = typename system::variables_tag;
     167             :   using analytic_solution_fields = typename variables_tag::tags_list;
     168             : 
     169             :   using initial_data_list = gh::Solutions::all_solutions<volume_dim>;
     170             : 
     171             :   using analytic_compute = evolution::Tags::AnalyticSolutionsCompute<
     172             :       volume_dim, analytic_solution_fields, false, initial_data_list>;
     173             :   using deriv_compute = ::Tags::DerivCompute<
     174             :       variables_tag, domain::Tags::Mesh<volume_dim>,
     175             :       domain::Tags::InverseJacobian<volume_dim, Frame::ElementLogical,
     176             :                                     Frame::Inertial>,
     177             :       typename system::gradient_variables>;
     178             :   using error_compute = Tags::ErrorsCompute<analytic_solution_fields>;
     179             :   using error_tags = db::wrap_tags_in<Tags::Error, analytic_solution_fields>;
     180             : 
     181             :   using observe_fields = tmpl::append<
     182             :       tmpl::push_back<
     183             :           analytic_solution_fields, gh::Tags::GaugeH<DataVector, volume_dim>,
     184             :           gh::Tags::SpacetimeDerivGaugeH<DataVector, volume_dim>,
     185             :           gr::Tags::SpatialMetric<DataVector, volume_dim>,
     186             :           gr::Tags::DetSpatialMetric<DataVector>,
     187             :           gr::Tags::InverseSpatialMetric<DataVector, volume_dim>,
     188             :           gr::Tags::Shift<DataVector, volume_dim>, gr::Tags::Lapse<DataVector>,
     189             :           gr::Tags::SqrtDetSpatialMetricCompute<DataVector, volume_dim,
     190             :                                                 Frame::Inertial>,
     191             :           gr::Tags::SpacetimeNormalOneFormCompute<DataVector, volume_dim,
     192             :                                                   Frame::Inertial>,
     193             :           gr::Tags::SpacetimeNormalVectorCompute<DataVector, volume_dim,
     194             :                                                  Frame::Inertial>,
     195             :           gr::Tags::InverseSpacetimeMetricCompute<DataVector, volume_dim,
     196             :                                                   Frame::Inertial>,
     197             : 
     198             :           gh::Tags::GaugeConstraintCompute<volume_dim, Frame::Inertial>,
     199             :           gh::Tags::TwoIndexConstraintCompute<volume_dim, Frame::Inertial>,
     200             :           gh::Tags::ThreeIndexConstraintCompute<volume_dim, Frame::Inertial>,
     201             :           gh::Tags::DerivSpatialMetricCompute<volume_dim, ::Frame::Inertial>,
     202             :           gr::Tags::SpatialChristoffelFirstKindCompute<DataVector, volume_dim,
     203             :                                                        ::Frame::Inertial>,
     204             :           gr::Tags::SpatialChristoffelSecondKindCompute<DataVector, volume_dim,
     205             :                                                         ::Frame::Inertial>,
     206             :           ::Tags::DerivTensorCompute<
     207             :               gr::Tags::SpatialChristoffelSecondKind<DataVector, volume_dim>,
     208             :               ::domain::Tags::InverseJacobian<volume_dim, Frame::ElementLogical,
     209             :                                               Frame::Inertial>,
     210             :               ::domain::Tags::Mesh<volume_dim>>,
     211             :           gr::Tags::SpatialRicciCompute<DataVector, volume_dim,
     212             :                                         ::Frame::Inertial>,
     213             :           gr::Tags::SpatialRicciScalarCompute<DataVector, volume_dim,
     214             :                                               ::Frame::Inertial>,
     215             :           // following tags added to observe constraints
     216             :           ::Tags::PointwiseL2NormCompute<
     217             :               gh::Tags::GaugeConstraint<DataVector, volume_dim>>,
     218             :           ::Tags::PointwiseL2NormCompute<
     219             :               gh::Tags::TwoIndexConstraint<DataVector, volume_dim>>,
     220             :           ::Tags::PointwiseL2NormCompute<
     221             :               gh::Tags::ThreeIndexConstraint<DataVector, volume_dim>>,
     222             :           ::domain::Tags::Coordinates<volume_dim, Frame::Grid>,
     223             :           ::domain::Tags::Coordinates<volume_dim, Frame::Inertial>>,
     224             :       error_tags,
     225             :       // The 4-index constraint is only implemented in 3d
     226             :       tmpl::conditional_t<
     227             :           volume_dim == 3,
     228             :           tmpl::list<
     229             :               gh::Tags::FourIndexConstraintCompute<3, Frame::Inertial>,
     230             :               gh::Tags::FConstraintCompute<3, Frame::Inertial>,
     231             :               ::Tags::PointwiseL2NormCompute<
     232             :                   gh::Tags::FConstraint<DataVector, 3>>,
     233             :               ::Tags::PointwiseL2NormCompute<
     234             :                   gh::Tags::FourIndexConstraint<DataVector, 3>>,
     235             :               gh::Tags::ConstraintEnergyCompute<3, Frame::Inertial>,
     236             :               gh::Tags::ExtrinsicCurvatureCompute<3, Frame::Inertial>,
     237             :               ::Tags::DerivTensorCompute<
     238             :                   gr::Tags::ExtrinsicCurvature<DataVector, 3>,
     239             :                   ::domain::Tags::InverseJacobian<
     240             :                       volume_dim, Frame::ElementLogical, Frame::Inertial>,
     241             :                   ::domain::Tags::Mesh<volume_dim>>,
     242             :               gr::Tags::WeylElectricCompute<DataVector, 3, Frame::Inertial>,
     243             :               gr::Tags::WeylElectricScalarCompute<DataVector, 3,
     244             :                                                   Frame::Inertial>,
     245             :               gr::Tags::WeylTypeD1Compute<DataVector, 3, Frame::Inertial>,
     246             :               gr::Tags::WeylTypeD1ScalarCompute<DataVector, 3, Frame::Inertial>,
     247             :               gr::Tags::Psi4RealCompute<Frame::Inertial>>,
     248             :           tmpl::list<>>>;
     249             :   using non_tensor_compute_tags = tmpl::list<
     250             :       ::Events::Tags::ObserverMeshCompute<volume_dim>,
     251             :       ::Events::Tags::ObserverCoordinatesCompute<volume_dim, Frame::Inertial>,
     252             :       ::Events::Tags::ObserverInverseJacobianCompute<
     253             :           volume_dim, Frame::ElementLogical, Frame::Inertial>,
     254             :       ::Events::Tags::ObserverJacobianCompute<volume_dim, Frame::ElementLogical,
     255             :                                               Frame::Inertial>,
     256             :       ::Events::Tags::ObserverDetInvJacobianCompute<Frame::ElementLogical,
     257             :                                                     Frame::Inertial>,
     258             :       ::Events::Tags::ObserverMeshVelocityCompute<volume_dim, Frame::Inertial>,
     259             :       analytic_compute, error_compute,
     260             :       gh::gauges::Tags::GaugeAndDerivativeCompute<volume_dim>>;
     261             : 
     262             :   using field_observations =
     263             :       dg::Events::field_observations<volume_dim, observe_fields,
     264             :                                      non_tensor_compute_tags>;
     265             : };
     266             : 
     267             : template <size_t volume_dim, bool LocalTimeStepping>
     268             : struct FactoryCreation : tt::ConformsTo<Options::protocols::FactoryCreation> {
     269             :   using system = gh::System<volume_dim>;
     270             : 
     271             :   using factory_classes = tmpl::map<
     272             :       tmpl::pair<
     273             :           amr::Criterion,
     274             :           tmpl::list<
     275             :               amr::Criteria::DriveToTarget<volume_dim>,
     276             :               amr::Criteria::Constraints<
     277             :                   volume_dim, tmpl::list<gh::Tags::ThreeIndexConstraintCompute<
     278             :                                   volume_dim, Frame::Inertial>>>,
     279             :               amr::Criteria::TruncationError<
     280             :                   volume_dim, typename system::variables_tag::tags_list>>>,
     281             :       tmpl::pair<DenseTrigger, DenseTriggers::standard_dense_triggers>,
     282             :       tmpl::pair<DomainCreator<volume_dim>, domain_creators<volume_dim>>,
     283             :       tmpl::pair<
     284             :           Event,
     285             :           tmpl::flatten<tmpl::list<
     286             :               Events::Completion, Events::MonitorMemory<volume_dim>,
     287             :               typename detail::ObserverTags<volume_dim>::field_observations,
     288             :               Events::time_events<system>,
     289             :               dg::Events::ObserveTimeStepVolume<volume_dim>>>>,
     290             :       tmpl::pair<
     291             :           gh::BoundaryConditions::BoundaryCondition<volume_dim>,
     292             :           gh::BoundaryConditions::standard_boundary_conditions<volume_dim>>,
     293             :       tmpl::pair<gh::gauges::GaugeCondition, gh::gauges::all_gauges>,
     294             :       tmpl::pair<
     295             :           evolution::initial_data::InitialData,
     296             :           tmpl::append<gh::Solutions::all_solutions<volume_dim>,
     297             :                        tmpl::conditional_t<volume_dim == 3,
     298             :                                            tmpl::list<gh::NumericInitialData>,
     299             :                                            tmpl::list<>>>>,
     300             :       tmpl::pair<LtsTimeStepper, TimeSteppers::lts_time_steppers>,
     301             :       tmpl::pair<MathFunction<1, Frame::Inertial>,
     302             :                  MathFunctions::all_math_functions<1, Frame::Inertial>>,
     303             :       tmpl::pair<PhaseChange, PhaseControl::factory_creatable_classes>,
     304             :       tmpl::pair<StepChooser<StepChooserUse::LtsStep>,
     305             :                  StepChoosers::standard_step_choosers<system>>,
     306             :       tmpl::pair<
     307             :           StepChooser<StepChooserUse::Slab>,
     308             :           StepChoosers::standard_slab_choosers<system, LocalTimeStepping>>,
     309             :       tmpl::pair<TimeSequence<double>,
     310             :                  TimeSequences::all_time_sequences<double>>,
     311             :       tmpl::pair<TimeSequence<std::uint64_t>,
     312             :                  TimeSequences::all_time_sequences<std::uint64_t>>,
     313             :       tmpl::pair<TimeStepper, TimeSteppers::time_steppers>,
     314             :       tmpl::pair<Trigger, tmpl::append<Triggers::logical_triggers,
     315             :                                        Triggers::time_triggers>>>;
     316             : };
     317             : }  // namespace detail
     318             : 
     319             : template <size_t VolumeDim, bool LocalTimeStepping>
     320           0 : struct GeneralizedHarmonicTemplateBase {
     321           0 :   static constexpr size_t volume_dim = VolumeDim;
     322           0 :   using system = gh::System<volume_dim>;
     323           0 :   using TimeStepperBase =
     324             :       tmpl::conditional_t<LocalTimeStepping, LtsTimeStepper, TimeStepper>;
     325             : 
     326           0 :   static constexpr bool local_time_stepping =
     327             :       TimeStepperBase::local_time_stepping;
     328           0 :   static constexpr bool use_dg_element_collection = false;
     329             : 
     330             :   // NOLINTNEXTLINE(google-runtime-references)
     331           0 :   void pup(PUP::er& /*p*/) {}
     332             : 
     333           0 :   using factory_creation =
     334             :       detail::FactoryCreation<volume_dim, local_time_stepping>;
     335             : 
     336           0 :   using observed_reduction_data_tags =
     337             :       observers::collect_reduction_data_tags<tmpl::push_back<
     338             :           tmpl::at<typename factory_creation::factory_classes, Event>>>;
     339             : 
     340           0 :   using initialize_initial_data_dependent_quantities_actions = tmpl::list<
     341             :       Actions::MutateApply<gh::gauges::SetPiAndPhiFromConstraints<volume_dim>>,
     342             :       Parallel::Actions::TerminatePhase>;
     343             : 
     344             :   // A tmpl::list of tags to be added to the GlobalCache by the
     345             :   // metavariables
     346           0 :   using const_global_cache_tags =
     347             :       tmpl::list<gh::gauges::Tags::GaugeCondition,
     348             :                  evolution::initial_data::Tags::InitialData,
     349             :                  gh::ConstraintDamping::Tags::DampingFunctionGamma0<
     350             :                      volume_dim, Frame::Grid>,
     351             :                  gh::ConstraintDamping::Tags::DampingFunctionGamma1<
     352             :                      volume_dim, Frame::Grid>,
     353             :                  gh::ConstraintDamping::Tags::DampingFunctionGamma2<
     354             :                      volume_dim, Frame::Grid>>;
     355             : 
     356           0 :   using dg_registration_list =
     357             :       tmpl::list<observers::Actions::RegisterEventsWithObservers>;
     358             : 
     359             :   // Register needs to be before InitializeTimeStepperHistory so that CCE is
     360             :   // properly registered when the self-start happens
     361           0 :   static constexpr auto default_phase_order =
     362             :       std::array{Parallel::Phase::Initialization,
     363             :                  Parallel::Phase::RegisterWithElementDataReader,
     364             :                  Parallel::Phase::ImportInitialData,
     365             :                  Parallel::Phase::InitializeInitialDataDependentQuantities,
     366             :                  Parallel::Phase::Register,
     367             :                  Parallel::Phase::InitializeTimeStepperHistory,
     368             :                  Parallel::Phase::CheckDomain,
     369             :                  Parallel::Phase::Evolve,
     370             :                  Parallel::Phase::Exit};
     371             : 
     372             :   template <typename ControlSystems>
     373           0 :   using step_actions = tmpl::list<
     374             :       evolution::dg::Actions::ComputeTimeDerivative<
     375             :           volume_dim, system, AllStepChoosers, local_time_stepping,
     376             :           use_dg_element_collection>,
     377             :       tmpl::conditional_t<
     378             :           local_time_stepping,
     379             :           tmpl::list<evolution::Actions::RunEventsAndDenseTriggers<tmpl::list<
     380             :                          ::domain::CheckFunctionsOfTimeAreReadyPostprocessor<
     381             :                              volume_dim>,
     382             :                          evolution::dg::ApplyBoundaryCorrections<
     383             :                              local_time_stepping, system, volume_dim, true>>>,
     384             :                      evolution::dg::Actions::ApplyLtsBoundaryCorrections<
     385             :                          system, volume_dim, false, use_dg_element_collection>>,
     386             :           tmpl::list<
     387             :               evolution::dg::Actions::ApplyBoundaryCorrectionsToTimeDerivative<
     388             :                   system, volume_dim, false, use_dg_element_collection>,
     389             :               Actions::RecordTimeStepperData<system>,
     390             :               evolution::Actions::RunEventsAndDenseTriggers<tmpl::list<>>,
     391             :               control_system::Actions::LimitTimeStep<ControlSystems>,
     392             :               Actions::UpdateU<system>>>,
     393             :       Actions::CleanHistory<system, local_time_stepping>,
     394             :       dg::Actions::Filter<
     395             :           Filters::Exponential<0>,
     396             :           tmpl::list<gr::Tags::SpacetimeMetric<DataVector, volume_dim>,
     397             :                      gh::Tags::Pi<DataVector, volume_dim>,
     398             :                      gh::Tags::Phi<DataVector, volume_dim>>>>;
     399             : 
     400             :   template <typename DerivedMetavars, bool UseControlSystems>
     401           0 :   using initialization_actions = tmpl::list<
     402             :       Initialization::Actions::InitializeItems<
     403             :           Initialization::TimeStepping<DerivedMetavars, TimeStepperBase>,
     404             :           evolution::dg::Initialization::Domain<volume_dim, UseControlSystems>,
     405             :           ::amr::Initialization::Initialize<volume_dim>,
     406             :           Initialization::TimeStepperHistory<DerivedMetavars>>,
     407             :       Initialization::Actions::NonconservativeSystem<system>,
     408             :       Initialization::Actions::AddComputeTags<::Tags::DerivCompute<
     409             :           typename system::variables_tag, domain::Tags::Mesh<volume_dim>,
     410             :           domain::Tags::InverseJacobian<volume_dim, Frame::ElementLogical,
     411             :                                         Frame::Inertial>,
     412             :           typename system::gradient_variables>>,
     413             :       gh::Actions::InitializeGhAnd3Plus1Variables<volume_dim>,
     414             :       Initialization::Actions::AddComputeTags<
     415             :           tmpl::push_back<StepChoosers::step_chooser_compute_tags<
     416             :               GeneralizedHarmonicTemplateBase, local_time_stepping>>>,
     417             :       ::evolution::dg::Initialization::Mortars<volume_dim, system>,
     418             :       evolution::Actions::InitializeRunEventsAndDenseTriggers,
     419             :       Parallel::Actions::TerminatePhase>;
     420             : };

Generated by: LCOV version 1.14