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 "Domain/Structure/Element.hpp"
       9             : #include "NumericalAlgorithms/DiscontinuousGalerkin/MortarHelpers.hpp"
      10             : #include "Utilities/TaggedTuple.hpp"
      11             : 
      12             : namespace dg {
      13             : 
      14             : /*!
      15             :  * \ingroup DiscontinuousGalerkinGroup
      16             :  * \brief Determines if data on all mortars has been received for the `InboxTag`
      17             :  * at time `temporal_id`.
      18             :  *
      19             :  * The `InboxTag` must hold a container indexed by the `TemporalIdType`,
      20             :  * representing received neighbor data at a given time. Each element in the
      21             :  * container, i.e. the received neighbor data at a given time, must be another
      22             :  * container indexed by the `dg::MortarId<Dim>`. The value it holds is not
      23             :  * relevant for this function. Here's an example for such a type:
      24             :  *
      25             :  * \snippet Test_HasReceivedFromAllMortars.cpp inbox_example
      26             :  */
      27             : template <typename InboxTag, size_t Dim, typename TemporalIdType,
      28             :           typename... InboxTags>
      29           1 : bool has_received_from_all_mortars(
      30             :     const TemporalIdType& temporal_id, const Element<Dim>& element,
      31             :     const tuples::TaggedTuple<InboxTags...>& inboxes) {
      32             :   if (UNLIKELY(element.number_of_neighbors() == 0)) {
      33             :     return true;
      34             :   }
      35             :   const auto& inbox = tuples::get<InboxTag>(inboxes);
      36             :   const auto temporal_received = inbox.find(temporal_id);
      37             :   if (temporal_received == inbox.end()) {
      38             :     return false;
      39             :   }
      40             :   const auto& received_neighbor_data = temporal_received->second;
      41             :   for (const auto& direction_and_neighbors : element.neighbors()) {
      42             :     const auto& direction = direction_and_neighbors.first;
      43             :     for (const auto& neighbor : direction_and_neighbors.second) {
      44             :       const auto neighbor_received =
      45             :           received_neighbor_data.find(MortarId<Dim>{direction, neighbor});
      46             :       if (neighbor_received == received_neighbor_data.end()) {
      47             :         return false;
      48             :       }
      49             :     }
      50             :   }
      51             :   return true;
      52             : }
      53             : 
      54             : }  // namespace dg