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

Generated by: LCOV version 1.14