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             : 
       8             : #include "DataStructures/DataBox/DataBox.hpp"
       9             : #include "Evolution/Systems/GeneralizedHarmonic/Bbh/CompletionCriteria.hpp"
      10             : #include "Evolution/Systems/GeneralizedHarmonic/Bbh/CompletionSingleton.hpp"
      11             : #include "NumericalAlgorithms/SphericalHarmonics/Tags.hpp"
      12             : #include "Parallel/GlobalCache.hpp"
      13             : #include "Parallel/Invoke.hpp"
      14             : #include "ParallelAlgorithms/ApparentHorizonFinder/FastFlow.hpp"
      15             : #include "ParallelAlgorithms/ApparentHorizonFinder/Protocols/Callback.hpp"
      16             : #include "ParallelAlgorithms/ApparentHorizonFinder/Tags.hpp"
      17             : #include "Utilities/ErrorHandling/Assert.hpp"
      18             : #include "Utilities/ProtocolHelpers.hpp"
      19             : #include "Utilities/TMPL.hpp"
      20             : 
      21           0 : namespace gh::bbh::callbacks {
      22             : /*!
      23             :  * \brief Apparent-horizon callback that updates binary-black-hole completion
      24             :  * criteria from successful common-horizon finds.
      25             :  *
      26             :  * \details This callback forwards each successful AhC find (`temporal_id`,
      27             :  * `LMax`) to `gh::bbh::CompletionSingleton`, where completion criteria are
      28             :  * evaluated in a temporally robust, centralized way.
      29             :  */
      30             : template <typename HorizonMetavars>
      31           1 : struct UpdateCompletionCriteria : tt::ConformsTo<ah::protocols::Callback> {
      32           0 :   using const_global_cache_tags = tmpl::list<>;
      33           0 :   using mutable_global_cache_tags = tmpl::list<>;
      34             : 
      35             :   template <typename DbTags, typename Metavariables>
      36           0 :   static void apply(const db::DataBox<DbTags>& box,
      37             :                     Parallel::GlobalCache<Metavariables>& cache,
      38             :                     const FastFlow::Status /*status*/) {
      39             :     const auto& current_time = db::get<ah::Tags::CurrentTime>(box);
      40             :     ASSERT(current_time.has_value(),
      41             :            "AhC completion callback requires a current temporal id.");
      42             :     const size_t l_max =
      43             :         db::get<ylm::Tags::Strahlkorper<typename HorizonMetavars::frame>>(box)
      44             :             .l_max();
      45             :     Parallel::simple_action<gh::bbh::Actions::RecordCommonHorizonSuccess>(
      46             :         Parallel::get_parallel_component<
      47             :             gh::bbh::CompletionSingleton<Metavariables>>(cache),
      48             :         *current_time, l_max);
      49             :   }
      50             : };
      51             : }  // namespace gh::bbh::callbacks