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 <tuple>
       8             : #include <utility>
       9             : 
      10             : #include "DataStructures/DataBox/DataBox.hpp"
      11             : #include "Evolution/Systems/Cce/NewmanPenrose.hpp"
      12             : #include "Evolution/Systems/Cce/PreSwshDerivatives.hpp"
      13             : 
      14             : namespace Cce {
      15             : namespace Actions {
      16             : 
      17             : /*!
      18             :  * \ingroup ActionsGroup
      19             :  * \brief Calculate \f$\Psi_0\f$ and its radial derivative
      20             :  * \f$\partial_\underline{\lambda}\Psi_0\f$ at the inner boundary of
      21             :  * the CCE grid.
      22             :  *
      23             :  * \details This action is to be called after `Tags::BondiJ` has been updated.
      24             :  * Note that `Tags::BondiJ` is in the asymptotically inertial frame, while the
      25             :  * returned \f$\Psi_0\f$ and its radial derivative are in the Cauchy frame.
      26             :  * \f$\underline \lambda\f$ is an affine parameter along the null
      27             :  * rays generated by \f$l\f$, see Eq. (19a) of \cite Moxon2020gha.
      28             :  */
      29           1 : struct CalculatePsi0AndDerivAtInnerBoundary {
      30           0 :   using const_global_cache_tags = tmpl::list<Tags::LMax>;
      31           0 :   using mutators =
      32             :       tmpl::list<TransformBondiJToCauchyCoords,
      33             :                  PreSwshDerivatives<Tags::Dy<Tags::BondiJCauchyView>>,
      34             :                  PreSwshDerivatives<Tags::Dy<Tags::Dy<Tags::BondiJCauchyView>>>,
      35             :                  VolumeWeyl<Tags::Psi0Match>,
      36             :                  PreSwshDerivatives<Tags::Dy<Tags::Psi0Match>>,
      37             :                  InnerBoundaryWeyl>;
      38             : 
      39             :   template <typename DbTags, typename... InboxTags, typename Metavariables,
      40             :             typename ArrayIndex, typename ActionList,
      41             :             typename ParallelComponent>
      42           0 :   static Parallel::iterable_action_return_t apply(
      43             :       db::DataBox<DbTags>& box,
      44             :       const tuples::TaggedTuple<InboxTags...>& /*inboxes*/,
      45             :       const Parallel::GlobalCache<Metavariables>& /*cache*/,
      46             :       const ArrayIndex& /*array_index*/, const ActionList /*meta*/,
      47             :       const ParallelComponent* const /*meta*/) {
      48             :     tmpl::for_each<mutators>([&box](auto mutator_v) {
      49             :       using mutator = typename decltype(mutator_v)::type;
      50             :       db::mutate_apply<mutator>(make_not_null(&box));
      51             :     });
      52             :     return {Parallel::AlgorithmExecution::Continue, std::nullopt};
      53             :   }
      54             : };
      55             : }  // namespace Actions
      56             : }  // namespace Cce