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 <optional>
       8             : #include <tuple>
       9             : #include <utility>
      10             : 
      11             : #include "DataStructures/DataBox/DataBox.hpp"
      12             : #include "Evolution/Systems/Cce/GaugeTransformBoundaryData.hpp"
      13             : #include "Evolution/Systems/Cce/OptionTags.hpp"
      14             : #include "Parallel/AlgorithmExecution.hpp"
      15             : #include "Parallel/GlobalCache.hpp"
      16             : #include "Utilities/Gsl.hpp"
      17             : #include "Utilities/TMPL.hpp"
      18             : 
      19             : namespace Cce {
      20             : namespace Actions {
      21             : 
      22             : /*!
      23             :  * \ingroup ActionsGroup
      24             :  * \brief Updates all of the gauge quantities associated with the additional
      25             :  * regularity-preserving gauge transformation on the boundaries for a new set of
      26             :  * Cauchy and partially flat Bondi-like coordinates.
      27             :  *
      28             :  * \details This action is to be called after `Tags::CauchyCartesianCoords`
      29             :  * and `Tags::PartiallyFlatCartesianCoords` have been updated, typically via a
      30             :  * time step of a set of coordinate evolution equations. It prepares the
      31             :  * gauge quantities in the \ref DataBoxGroup for calls to the individual
      32             :  * `GaugeAdjustedBoundaryValue` specializations.
      33             :  *
      34             :  * Internally, this dispatches to `GaugeUpdateAngularFromCartesian`,
      35             :  * `GaugeUpdateJacobianFromCoordinates`, `GaugeUpdateInterpolator`, and
      36             :  * `GaugeUpdateOmega` to perform the computations. Refer to the documentation
      37             :  * for those mutators for mathematical details.
      38             :  */
      39             : template <bool EvolveCcm>
      40           1 : struct UpdateGauge {
      41           0 :   using const_global_cache_tags = tmpl::list<Tags::LMax>;
      42             : 
      43             :   template <typename DbTags, typename... InboxTags, typename Metavariables,
      44             :             typename ArrayIndex, typename ActionList,
      45             :             typename ParallelComponent>
      46           0 :   static Parallel::iterable_action_return_t apply(
      47             :       db::DataBox<DbTags>& box,
      48             :       const tuples::TaggedTuple<InboxTags...>& /*inboxes*/,
      49             :       const Parallel::GlobalCache<Metavariables>& /*cache*/,
      50             :       const ArrayIndex& /*array_index*/, const ActionList /*meta*/,
      51             :       const ParallelComponent* const /*meta*/) {
      52             :     db::mutate_apply<GaugeUpdateAngularFromCartesian<
      53             :         Tags::CauchyAngularCoords, Tags::CauchyCartesianCoords>>(
      54             :         make_not_null(&box));
      55             :     db::mutate_apply<GaugeUpdateJacobianFromCoordinates<
      56             :         Tags::PartiallyFlatGaugeC, Tags::PartiallyFlatGaugeD,
      57             :         Tags::CauchyAngularCoords, Tags::CauchyCartesianCoords>>(
      58             :         make_not_null(&box));
      59             :     db::mutate_apply<GaugeUpdateInterpolator<Tags::CauchyAngularCoords>>(
      60             :         make_not_null(&box));
      61             :     db::mutate_apply<
      62             :         GaugeUpdateOmega<Tags::PartiallyFlatGaugeC, Tags::PartiallyFlatGaugeD,
      63             :                          Tags::PartiallyFlatGaugeOmega>>(make_not_null(&box));
      64             : 
      65             :     if constexpr (EvolveCcm) {
      66             :       db::mutate_apply<
      67             :           GaugeUpdateAngularFromCartesian<Tags::PartiallyFlatAngularCoords,
      68             :                                           Tags::PartiallyFlatCartesianCoords>>(
      69             :           make_not_null(&box));
      70             :       db::mutate_apply<GaugeUpdateJacobianFromCoordinates<
      71             :           Tags::CauchyGaugeC, Tags::CauchyGaugeD,
      72             :           Tags::PartiallyFlatAngularCoords,
      73             :           Tags::PartiallyFlatCartesianCoords>>(make_not_null(&box));
      74             :       db::mutate_apply<
      75             :           GaugeUpdateInterpolator<Tags::PartiallyFlatAngularCoords>>(
      76             :           make_not_null(&box));
      77             :       db::mutate_apply<GaugeUpdateOmega<Tags::CauchyGaugeC, Tags::CauchyGaugeD,
      78             :                                         Tags::CauchyGaugeOmega>>(
      79             :           make_not_null(&box));
      80             :     }
      81             :     return {Parallel::AlgorithmExecution::Continue, std::nullopt};
      82             :   }
      83             : };
      84             : }  // namespace Actions
      85             : }  // namespace Cce