Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include "DataStructures/DataBox/Protocols/Mutator.hpp"
       7              : #include "DataStructures/Tensor/Tensor.hpp"
       8              : #include "Domain/Structure/DirectionalIdMap.hpp"
       9              : #include "Evolution/DgSubcell/GhostData.hpp"
      10              : #include "Evolution/DgSubcell/Mesh.hpp"
      11              : #include "Evolution/DgSubcell/Tags/GhostDataForReconstruction.hpp"
      12              : #include "Evolution/DgSubcell/Tags/Mesh.hpp"
      13              : #include "Evolution/Systems/Ccz4/FiniteDifference/System.hpp"
      14              : #include "Evolution/Systems/Ccz4/FiniteDifference/Tags.hpp"
      15              : #include "Evolution/Systems/Ccz4/Tags.hpp"
      16              : #include "Utilities/Gsl.hpp"
      17              : #include "Utilities/ProtocolHelpers.hpp"
      18              : #include "Utilities/TMPL.hpp"
      19              : 
      20            1 : namespace Ccz4::fd {
      21              : /*!
      22              :  * \brief Apply the Kreiss-Oliger filter to the evolved variables
      23              :  *
      24              :  */
      25            1 : struct ApplyFilter : tt::ConformsTo<db::protocols::Mutator> {
      26            0 :   using return_tags = tmpl::list<System::variables_tag>;
      27            0 :   using argument_tags =
      28              :       tmpl::list<evolution::dg::subcell::Tags::Mesh<3>,
      29              :                  Tags::EvolveLapseAndShift, Tags::KreissOligerEpsilon,
      30              :                  evolution::dg::subcell::Tags::GhostDataForReconstruction<3>>;
      31              : 
      32            0 :   static void apply(
      33              :       gsl::not_null<typename System::variables_tag::type*>
      34              :           evolved_vars_ptr,
      35              :       const Mesh<3>& subcell_mesh, bool evolve_lapse_and_shift,
      36              :       double kreiss_oliger_epsilon,
      37              :       const DirectionalIdMap<3, evolution::dg::subcell::GhostData>& ghost_data);
      38              : };
      39              : }  // namespace Ccz4::fd