Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <optional>
       7             : #include <tuple>
       8             : 
       9             : #include "DataStructures/DataBox/DataBox.hpp"
      10             : #include "DataStructures/Variables.hpp"
      11             : #include "Evolution/Systems/Cce/ReceiveTags.hpp"
      12             : #include "Parallel/AlgorithmExecution.hpp"
      13             : #include "Parallel/GlobalCache.hpp"
      14             : #include "Parallel/Invoke.hpp"
      15             : #include "Time/TimeStepId.hpp"
      16             : #include "Utilities/Gsl.hpp"
      17             : #include "Utilities/TMPL.hpp"
      18             : 
      19             : /// \cond
      20             : namespace Tags {
      21             : struct TimeStepId;
      22             : }  // namespace Tags
      23             : /// \endcond
      24             : 
      25             : namespace Cce {
      26             : namespace Actions {
      27             : 
      28             : /*!
      29             :  * \ingroup ActionsGroup
      30             :  * \brief Takes the boundary data needed to perform the CCE linear solves as
      31             :  * arguments and puts them in the \ref DataBoxGroup, updating the
      32             :  * `Cce::Tags::BoundaryTime` accordingly.
      33             :  *
      34             :  * \details The boundary data is computed by a separate component, and packaged
      35             :  * into a `Variables<TagList>` which is sent in the argument
      36             :  * of the simple action invocation. The `TimeStepId` is also provided to confirm
      37             :  * the time associated with the passed boundary data.
      38             :  *
      39             :  * \ref DataBoxGroup changes:
      40             :  * - Adds: nothing
      41             :  * - Removes: nothing
      42             :  * - Modifies:
      43             :  *   - All tags in `TagList`
      44             :  */
      45             : template <typename Metavariables, typename TagList>
      46           1 : struct ReceiveWorldtubeData {
      47           0 :   using inbox_tags = tmpl::list<Cce::ReceiveTags::BoundaryData<TagList>>;
      48             : 
      49             :   template <typename DbTags, typename... InboxTags, typename ArrayIndex,
      50             :             typename ActionList, typename ParallelComponent>
      51           0 :   static Parallel::iterable_action_return_t apply(
      52             :       db::DataBox<DbTags>& box, tuples::TaggedTuple<InboxTags...>& inboxes,
      53             :       const Parallel::GlobalCache<Metavariables>& /*cache*/,
      54             :       const ArrayIndex& /*array_index*/, const ActionList /*meta*/,
      55             :       const ParallelComponent* const /*meta*/) {
      56             :     auto& inbox = tuples::get<Cce::ReceiveTags::BoundaryData<TagList>>(inboxes);
      57             :     if (inbox.count(db::get<::Tags::TimeStepId>(box)) != 1) {
      58             :       return {Parallel::AlgorithmExecution::Pause,
      59             :               tmpl::index_of<ActionList, ReceiveWorldtubeData>::value};
      60             :     }
      61             : 
      62             :     tmpl::for_each<TagList>(
      63             :         [&inbox, &box](auto tag_v) {
      64             :           using tag = typename decltype(tag_v)::type;
      65             :           db::mutate<tag>(
      66             :               [&inbox](const gsl::not_null<typename tag::type*> destination,
      67             :                        const TimeStepId& time) {
      68             :                 *destination = get<tag>(inbox[time]);
      69             :               },
      70             :               make_not_null(&box), db::get<::Tags::TimeStepId>(box));
      71             :         });
      72             :     inbox.erase(db::get<::Tags::TimeStepId>(box));
      73             :     return {Parallel::AlgorithmExecution::Continue,
      74             :             tmpl::index_of<ActionList, ReceiveWorldtubeData>::value + 1};
      75             :   }
      76             : };
      77             : }  // namespace Actions
      78             : }  // namespace Cce