Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <string>
       7              : #include <tuple>
       8              : 
       9              : #include "DataStructures/DataBox/DataBox.hpp"
      10              : #include "DataStructures/DataVector.hpp"
      11              : #include "DataStructures/LinkedMessageId.hpp"
      12              : #include "DataStructures/Tensor/Tensor.hpp"
      13              : #include "DataStructures/Variables.hpp"
      14              : #include "Evolution/Systems/Cce/Actions/BoundaryComputeAndSendToEvolution.hpp"
      15              : #include "Evolution/Systems/Cce/InterfaceManagers/GhInterfaceManager.hpp"
      16              : #include "Evolution/Systems/Cce/OptionTags.hpp"
      17              : #include "Evolution/Systems/Cce/WorldtubeDataManager.hpp"
      18              : #include "Parallel/GlobalCache.hpp"
      19              : #include "Parallel/Invoke.hpp"
      20              : #include "Time/TimeStepId.hpp"
      21              : #include "Utilities/Gsl.hpp"
      22              : #include "Utilities/Requires.hpp"
      23              : #include "Utilities/TMPL.hpp"
      24              : #include "Utilities/TypeTraits.hpp"
      25              : 
      26              : namespace Cce {
      27              : namespace Actions {
      28              : 
      29              : /*!
      30              :  * \ingroup ActionsGroup
      31              :  * \brief Stores the boundary data from the GH evolution in the
      32              :  * `Cce::InterfaceManagers::GhInterfaceManager`, and sends to the
      33              :  * `EvolutionComponent` (template argument) if the data fulfills a prior
      34              :  * request.
      35              :  *
      36              :  * \details If the new data fulfills a prior request submitted to the
      37              :  * `Cce::InterfaceManagers::GhInterfaceManager`, this will dispatch the result
      38              :  * to `Cce::Actions::SendToEvolution<GhWorldtubeBoundary<Metavariables>,
      39              :  * EvolutionComponent>` for sending the processed boundary data to
      40              :  * the `EvolutionComponent`.
      41              :  *
      42              :  * \ref DataBoxGroup changes:
      43              :  * - Adds: nothing
      44              :  * - Removes: nothing
      45              :  * - Modifies:
      46              :  *   - `Tags::GhInterfaceManager`
      47              :  */
      48              : template <typename EvolutionComponent, bool DuringSelfStart>
      49            1 : struct ReceiveGhWorldtubeData {
      50              :   template <typename ParallelComponent, typename... DbTags,
      51              :             typename Metavariables, typename ArrayIndex>
      52            0 :   static void apply(db::DataBox<tmpl::list<DbTags...>>& box,
      53              :                     Parallel::GlobalCache<Metavariables>& cache,
      54              :                     const ArrayIndex& /*array_index*/,
      55              :                     const tmpl::conditional_t<DuringSelfStart, TimeStepId,
      56              :                                               LinkedMessageId<double>>
      57              :                         time,
      58              :                     const tnsr::aa<DataVector, 3>& spacetime_metric,
      59              :                     const tnsr::iaa<DataVector, 3>& phi,
      60              :                     const tnsr::aa<DataVector, 3>& pi) {
      61              :     auto insert_gh_data_to_interface_manager =
      62              :         [&spacetime_metric, &phi, &pi, &time,
      63              :          &cache](const auto interface_manager) {
      64              :           interface_manager->insert_gh_data(time, spacetime_metric, phi, pi);
      65              :           const auto gh_data =
      66              :               interface_manager->retrieve_and_remove_first_ready_gh_data();
      67              :           if (static_cast<bool>(gh_data)) {
      68              :             Parallel::simple_action<Actions::SendToEvolution<
      69              :                 GhWorldtubeBoundary<Metavariables>, EvolutionComponent>>(
      70              :                 Parallel::get_parallel_component<
      71              :                     GhWorldtubeBoundary<Metavariables>>(cache),
      72              :                 get<0>(*gh_data), get<1>(*gh_data));
      73              :           }
      74              :         };
      75              :     if constexpr (DuringSelfStart) {
      76              :       db::mutate<Tags::SelfStartGhInterfaceManager>(
      77              :           insert_gh_data_to_interface_manager, make_not_null(&box));
      78              :     } else {
      79              :       db::mutate<Tags::GhInterfaceManager>(insert_gh_data_to_interface_manager,
      80              :                                            make_not_null(&box));
      81              :     }
      82              :   }
      83              : };
      84              : }  // namespace Actions
      85              : }  // namespace Cce