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             : #include <iosfwd>
       8             : #include <optional>
       9             : 
      10             : #include "DataStructures/DataVector.hpp"
      11             : #include "Evolution/DiscontinuousGalerkin/InterpolatedBoundaryData.hpp"
      12             : #include "NumericalAlgorithms/Spectral/Mesh.hpp"
      13             : #include "Time/TimeStepId.hpp"
      14             : 
      15             : namespace evolution::dg {
      16             : /*!
      17             :  * \brief The data communicated between neighber elements.
      18             :  *
      19             :  * The stored data consists of the following:
      20             :  *
      21             :  * 1. the volume mesh of the element.
      22             :  * 2. the volume mesh corresponding to the ghost cell data. This allows eliding
      23             :  *    projection when all neighboring elements are doing DG.
      24             :  * 3. the mortar mesh of the data on the mortar
      25             :  * 4. the variables at the ghost zone cells for finite difference/volume
      26             :  *    reconstruction
      27             :  * 5. the data on the mortar needed for computing the boundary corrections (e.g.
      28             :  *    fluxes, characteristic speeds, conserved variables)
      29             :  * 6. the TimeStepId beyond which the boundary terms are no longer valid, when
      30             :  *    using local time stepping.
      31             :  * 7. the troubled cell indicator status used for determining halos around
      32             :  *    troubled cells.
      33             :  * 8. the integration order of the time-stepper
      34             :  * 9. the InterpolatedBoundaryData sent by a non-conforming Element that
      35             :  *    interpolates its data to a subset of the points of the Element receiving
      36             :  *    this BoundaryData
      37             :  */
      38             : template <size_t Dim>
      39           1 : struct BoundaryData {
      40             :   // NOLINTNEXTLINE(google-runtime-references)
      41           0 :   void pup(PUP::er& p);
      42             : 
      43           0 :   Mesh<Dim> volume_mesh{};
      44           0 :   std::optional<Mesh<Dim>> volume_mesh_ghost_cell_data{};
      45           0 :   std::optional<Mesh<Dim - 1>> boundary_correction_mesh{};
      46           0 :   std::optional<DataVector> ghost_cell_data{};
      47           0 :   std::optional<DataVector> boundary_correction_data{};
      48           0 :   ::TimeStepId validity_range{};
      49           0 :   int tci_status{};
      50           0 :   size_t integration_order{std::numeric_limits<size_t>::max()};
      51           0 :   std::optional<InterpolatedBoundaryData<Dim>> interpolated_boundary_data{};
      52             : };
      53             : 
      54             : template <size_t Dim>
      55           0 : bool operator==(const BoundaryData<Dim>& lhs, const BoundaryData<Dim>& rhs);
      56             : template <size_t Dim>
      57           0 : bool operator!=(const BoundaryData<Dim>& lhs, const BoundaryData<Dim>& rhs);
      58             : template <size_t Dim>
      59           0 : std::ostream& operator<<(std::ostream& os, const BoundaryData<Dim>& value);
      60             : }  // namespace evolution::dg