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              : 
       8              : #include "DataStructures/LinkedMessageId.hpp"
       9              : #include "IO/Logging/Verbosity.hpp"
      10              : #include "Parallel/GlobalCache.hpp"
      11              : #include "Parallel/ParallelComponentHelpers.hpp"
      12              : #include "ParallelAlgorithms/ApparentHorizonFinder/FastFlow.hpp"
      13              : #include "ParallelAlgorithms/ApparentHorizonFinder/OptionTags.hpp"
      14              : #include "ParallelAlgorithms/ApparentHorizonFinder/Tags.hpp"
      15              : #include "Utilities/Gsl.hpp"
      16              : #include "Utilities/TMPL.hpp"
      17              : 
      18              : namespace ah {
      19              : /*!
      20              :  * \brief Initialize items related to the horizon finder
      21              :  *
      22              :  * GlobalCache:
      23              :  * - Uses:
      24              :  *   - `ah::Tags::ApparentHorizonOptions`
      25              :  *   - `ah::Tags::BlocksForHorizonFind`
      26              :  *   - `ah::Tags::LMax`
      27              :  *
      28              :  * DataBox:
      29              :  * - Uses: Nothing
      30              :  * - Adds:
      31              :  *   - `ah::Tags::Verbosity`
      32              :  *   - `ah::Tags::FastFlow`
      33              :  *   - `ah::Tags::CurrentTime`
      34              :  *   - `ah::Tags::PendingTimes`
      35              :  *   - `ah::Tags::CompletedTimes`
      36              :  *   - `ah::Tags::Storage`
      37              :  *   - `ah::Tags::PreviousSurfaces`
      38              :  *   - `ah::Tags::Strahlkorper`
      39              :  *   - `ah::Tags::TimeDerivStrahlkorper`
      40              :  *   - `ah::Tags::Dependency`
      41              :  *   - `::Tags::Variables<ah::vars_to_interpolate_to_target>`
      42              :  *   - `ah::Tags::CurrentResolutionL`
      43              :  * - Modifies:
      44              :  *   - `ah::Tags::Verbosity`
      45              :  *   - `ah::Tags::FastFlow`
      46              :  *   - `ah::Tags::CurrentTime`
      47              :  */
      48              : template <typename HorizonMetavars>
      49            1 : struct Initialize {
      50            0 :   using Fr = typename HorizonMetavars::frame;
      51              : 
      52            0 :   using simple_tags_from_options = tmpl::list<>;
      53              : 
      54            0 :   using simple_tags = tmpl::append<
      55              :       tmpl::list<Tags::CurrentResolutionL, Tags::Verbosity, Tags::FastFlow,
      56              :                  Tags::CurrentTime, Tags::PendingTimes, Tags::CompletedTimes,
      57              :                  Tags::Storage<Fr>, Tags::BlockSearchOrder,
      58              :                  Tags::PreviousSurfaces<Fr>, ylm::Tags::Strahlkorper<Fr>,
      59              :                  ylm::Tags::TimeDerivStrahlkorper<Fr>, ah::Tags::Dependency,
      60              :                  ::Tags::Variables<ah::vars_to_interpolate_to_target<3, Fr>>>,
      61              :       tmpl::conditional_t<
      62              :           std::is_same_v<Fr, Frame::Inertial>, tmpl::list<>,
      63              :           tmpl::list<ylm::Tags::CartesianCoords<Frame::Inertial>>>>;
      64              : 
      65            0 :   using const_global_cache_tags =
      66              :       tmpl::remove_duplicates<tmpl::flatten<tmpl::append<
      67              :           tmpl::list<Tags::ApparentHorizonOptions<HorizonMetavars>,
      68              :                      Tags::BlocksForHorizonFind, Tags::LMax>,
      69              :           Parallel::get_const_global_cache_tags_from_actions<tmpl::flatten<
      70              :               tmpl::list<typename HorizonMetavars::horizon_find_callbacks,
      71              :                          typename HorizonMetavars::
      72              :                              horizon_find_failure_callbacks>>>>>>;
      73              : 
      74            0 :   using mutable_global_cache_tags =
      75              :       tmpl::list<Tags::PreviousSurface<HorizonMetavars>>;
      76              : 
      77            0 :   using compute_tags = ah::compute_items_on_target<3, Fr>;
      78              : 
      79            0 :   using return_tags = tmpl::list<Tags::CurrentResolutionL, Tags::Verbosity,
      80              :                                  Tags::FastFlow, Tags::CurrentTime>;
      81              : 
      82            0 :   using argument_tags =
      83              :       tmpl::list<Tags::ApparentHorizonOptions<HorizonMetavars>>;
      84              : 
      85            0 :   static void apply(
      86              :       const gsl::not_null<std::optional<size_t>*> current_resolution_l,
      87              :       const gsl::not_null<::Verbosity*> verbosity,
      88              :       const gsl::not_null<::FastFlow*> fast_flow,
      89              :       const gsl::not_null<std::optional<LinkedMessageId<double>>*> current_time,
      90              :       const HorizonOptions<Fr>& options) {
      91              :     (*verbosity) = options.verbosity;
      92              :     (*fast_flow) = options.fast_flow;
      93              :     current_time->reset();
      94              :     current_resolution_l->reset();
      95              :   }
      96              : };
      97              : }  // namespace ah