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/Tag.hpp"
       9             : #include "DataStructures/Tensor/TypeAliases.hpp"
      10             : #include "Domain/Structure/Direction.hpp"
      11             : #include "Domain/Structure/DirectionalId.hpp"
      12             : #include "Domain/Structure/DirectionalIdMap.hpp"
      13             : #include "Domain/Structure/ElementId.hpp"
      14             : #include "Domain/Tags.hpp"
      15             : #include "Evolution/DgSubcell/ActiveGrid.hpp"
      16             : #include "Evolution/DgSubcell/Tags/ActiveGrid.hpp"
      17             : #include "Evolution/DgSubcell/Tags/Mesh.hpp"
      18             : #include "Utilities/Gsl.hpp"
      19             : 
      20             : /// \cond
      21             : class DataVector;
      22             : template <size_t Dim>
      23             : class Mesh;
      24             : /// \endcond
      25             : 
      26             : namespace evolution::dg::subcell::Tags {
      27             : /// Stores the status of the troubled cell indicator in the element as an `int`.
      28             : ///
      29             : /// A non-zero value indicates the TCI decided the element is troubled.
      30           1 : struct TciDecision : db::SimpleTag {
      31           0 :   using type = int;
      32             : };
      33             : 
      34             : /// The TCI decision of neighboring elements.
      35             : template <size_t Dim>
      36           1 : struct NeighborTciDecisions : db::SimpleTag {
      37           0 :   using type = DirectionalIdMap<Dim, int>;
      38             : };
      39             : 
      40             : /// Stores the status of the troubled cell indicator in the element
      41             : /// (TciDecision) as a `Scalar<DataVector>` so it can be observed.
      42           1 : struct TciStatus : db::SimpleTag {
      43           0 :   using type = Scalar<DataVector>;
      44             : };
      45             : 
      46             : /// Compute tag to get a `TciStatus` from a `TciDecision`.
      47             : template <size_t Dim>
      48           1 : struct TciStatusCompute : db::ComputeTag, TciStatus {
      49           0 :   using base = TciStatus;
      50           0 :   using return_type = typename base::type;
      51           0 :   using argument_tags = tmpl::list<Tags::TciDecision, Tags::ActiveGrid,
      52             :                                    Tags::Mesh<Dim>, ::domain::Tags::Mesh<Dim>>;
      53           0 :   static void function(gsl::not_null<return_type*> result, int tci_decision,
      54             :                        subcell::ActiveGrid active_grid,
      55             :                        const ::Mesh<Dim>& subcell_mesh,
      56             :                        const ::Mesh<Dim>& dg_mesh);
      57             : };
      58             : }  // namespace evolution::dg::subcell::Tags