Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <boost/functional/hash.hpp>
       7             : #include <cstddef>
       8             : #include <iomanip>
       9             : #include <map>
      10             : #include <sstream>
      11             : #include <string>
      12             : #include <unordered_map>
      13             : #include <utility>
      14             : 
      15             : #include "Domain/Structure/DirectionalId.hpp"
      16             : #include "Evolution/DgSubcell/InitialTciData.hpp"
      17             : #include "Utilities/Gsl.hpp"
      18             : 
      19             : namespace evolution::dg::subcell::Tags {
      20             : /*!
      21             :  * \brief Inbox tag for communicating the RDMP and TCI status/decision during
      22             :  * initialization.
      23             :  *
      24             :  * \note The inner map uses `std::unordered_map` rather than
      25             :  * `DirectionalIdMap` because elements at non-conforming interfaces (e.g.
      26             :  *  inner wedges bordering a spherical shell) can have more than
      27             :  * `maximum_number_of_neighbors(Dim)` neighbors sending messages. This inbox
      28             :  * is only populated during initialization and is erased after use, so the
      29             :  * heap-allocation cost is negligible.
      30             :  */
      31             : template <size_t Dim>
      32           1 : struct InitialTciData {
      33           0 :   using temporal_id = int;
      34           0 :   using type =
      35             :       std::map<temporal_id,
      36             :                std::unordered_map<DirectionalId<Dim>,
      37             :                                   evolution::dg::subcell::InitialTciData,
      38             :                                   boost::hash<DirectionalId<Dim>>>>;
      39             : 
      40             :   template <typename ReceiveDataType>
      41           0 :   static bool insert_into_inbox(const gsl::not_null<type*> inbox,
      42             :                                 const temporal_id& time_step_id,
      43             :                                 ReceiveDataType&& data) {
      44             :     auto& current_inbox = (*inbox)[time_step_id];
      45             : 
      46             :     const auto& direction_and_element_id = data.first;
      47             : 
      48             :     ASSERT(current_inbox.find(direction_and_element_id) == current_inbox.end(),
      49             :            "Received data from direction "
      50             :                << direction_and_element_id.direction() << " and element ID "
      51             :                << direction_and_element_id.id() << " more than once");
      52             :     current_inbox.emplace(std::forward<ReceiveDataType>(data));
      53             :     return true;
      54             :   }
      55             : 
      56           0 :   static std::string output_inbox(const type& inbox,
      57             :                                   const size_t padding_size) {
      58             :     std::stringstream ss{};
      59             :     const std::string pad(padding_size, ' ');
      60             : 
      61             :     ss << std::scientific << std::setprecision(16);
      62             :     ss << pad << "InitialTciDataInbox:\n";
      63             :     for (const auto& [action_number, hash_map] : inbox) {
      64             :       ss << pad << " Action number: " << action_number << "\n";
      65             :       for (const auto& [key, initial_tci_data] : hash_map) {
      66             :         using ::operator<<;
      67             :         ss << pad << "  Key: " << key
      68             :            << ", TCI: " << initial_tci_data.tci_status << "\n";
      69             :       }
      70             :     }
      71             : 
      72             :     return ss.str();
      73             :   }
      74             : };
      75             : }  // namespace evolution::dg::subcell::Tags