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 <tuple>
       8             : 
       9             : #include "DataStructures/Tensor/TypeAliases.hpp"
      10             : #include "Domain/Tags.hpp"
      11             : #include "Evolution/DgSubcell/RdmpTciData.hpp"
      12             : #include "Evolution/DgSubcell/Tags/DataForRdmpTci.hpp"
      13             : #include "Evolution/DgSubcell/Tags/Mesh.hpp"
      14             : #include "Evolution/DgSubcell/Tags/SubcellOptions.hpp"
      15             : #include "Evolution/Systems/ScalarAdvection/Subcell/TciOptions.hpp"
      16             : #include "Evolution/Systems/ScalarAdvection/Tags.hpp"
      17             : #include "Utilities/TMPL.hpp"
      18             : 
      19             : /// \cond
      20             : class DataVector;
      21             : template <size_t Dim>
      22             : class Mesh;
      23             : /// \endcond
      24             : 
      25             : namespace ScalarAdvection::subcell {
      26             : /*!
      27             :  * \brief Troubled-cell indicator applied to the finite difference subcell
      28             :  * solution to check if the corresponding DG solution is admissible.
      29             :  *
      30             :  * Applies 1) the RDMP TCI to \f$U\f$ and 2) the Persson TCI to \f$U\f$ if the
      31             :  * \f$\max(|U|)\f$ on the DG grid is greater than `tci_options.u_cutoff`.
      32             :  */
      33             : template <size_t Dim>
      34           1 : struct TciOnFdGrid {
      35           0 :   using return_tags = tmpl::list<>;
      36           0 :   using argument_tags =
      37             :       tmpl::list<ScalarAdvection::Tags::U, ::domain::Tags::Mesh<Dim>,
      38             :                  evolution::dg::subcell::Tags::Mesh<Dim>,
      39             :                  evolution::dg::subcell::Tags::DataForRdmpTci,
      40             :                  evolution::dg::subcell::Tags::SubcellOptions<Dim>,
      41             :                  Tags::TciOptions>;
      42             : 
      43           0 :   static std::tuple<bool, evolution::dg::subcell::RdmpTciData> apply(
      44             :       const Scalar<DataVector>& subcell_u, const Mesh<Dim>& dg_mesh,
      45             :       const Mesh<Dim>& subcell_mesh,
      46             :       const evolution::dg::subcell::RdmpTciData& past_rdmp_tci_data,
      47             :       const evolution::dg::subcell::SubcellOptions& subcell_options,
      48             :       const TciOptions& tci_options, double persson_exponent,
      49             :       bool need_rdmp_data_only);
      50             : };
      51             : }  // namespace ScalarAdvection::subcell