Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <optional>
       7             : #include <tuple>
       8             : #include <utility>
       9             : 
      10             : #include "DataStructures/DataBox/DataBox.hpp"
      11             : #include "DataStructures/DataVector.hpp"
      12             : #include "DataStructures/Tensor/EagerMath/Determinant.hpp"
      13             : #include "DataStructures/Tensor/Tensor.hpp"
      14             : #include "Domain/CoordinateMaps/CoordinateMap.hpp"
      15             : #include "Domain/CoordinateMaps/Tags.hpp"
      16             : #include "Domain/ElementMap.hpp"
      17             : #include "Domain/Tags.hpp"
      18             : #include "Evolution/DgSubcell/Tags/Coordinates.hpp"
      19             : #include "Evolution/DgSubcell/Tags/Jacobians.hpp"
      20             : #include "Evolution/DiscontinuousGalerkin/MortarData.hpp"
      21             : #include "Evolution/DiscontinuousGalerkin/MortarDataHolder.hpp"
      22             : #include "Evolution/DiscontinuousGalerkin/MortarTags.hpp"
      23             : #include "Parallel/AlgorithmExecution.hpp"
      24             : #include "Utilities/ErrorHandling/Assert.hpp"
      25             : 
      26             : /// \cond
      27             : namespace Parallel {
      28             : template <typename Metavariables>
      29             : class GlobalCache;
      30             : }  // namespace Parallel
      31             : namespace tuples {
      32             : template <typename...>
      33             : class TaggedTuple;
      34             : }  // namespace tuples
      35             : /// \endcond
      36             : 
      37             : namespace evolution::dg::subcell::fd::Actions {
      38             : /*!
      39             :  * \brief Take a finite-difference time step on the subcell grid.
      40             :  *
      41             :  * The template parameter `TimeDerivative` must have a `static apply` function
      42             :  * that takes the `DataBox` by `gsl::not_null` as the first argument, the
      43             :  * cell-centered inverse Jacobian from the logical to the grid frame as the
      44             :  * second argument, and its determinant as the third argument.
      45             :  *
      46             :  * GlobalCache: nothing
      47             :  *
      48             :  * DataBox:
      49             :  * - Uses:
      50             :  *   - `subcell::fd::Tags::InverseJacobianLogicalToGrid<Dim>`
      51             :  *   - `subcell::fd::Tags::DetInverseJacobianLogicalToGrid`
      52             :  *   - `domain::Tags::ElementMap<Dim, Frame::Grid>`
      53             :  *   - `domain::CoordinateMaps::Tags::CoordinateMap<Dim, Grid, Inertial>`
      54             :  *   - `subcell::Tags::Coordinates<Dim, Frame::ElementLogical>`
      55             :  *   - Anything that `Metavariables::SubcellOptions::TimeDerivative` uses
      56             :  * - Adds: nothing
      57             :  * - Removes: nothing
      58             :  * - Modifies:
      59             :  *   - Anything that `Metavariables::SubcellOptions::TimeDerivative` modifies
      60             :  */
      61             : template <typename TimeDerivative>
      62           1 : struct TakeTimeStep {
      63             :   template <typename DbTags, typename... InboxTags, typename Metavariables,
      64             :             typename ArrayIndex, typename ActionList,
      65             :             typename ParallelComponent, size_t Dim = Metavariables::volume_dim>
      66           0 :   static Parallel::iterable_action_return_t apply(
      67             :       db::DataBox<DbTags>& box, tuples::TaggedTuple<InboxTags...>& /*inboxes*/,
      68             :       const Parallel::GlobalCache<Metavariables>& /*cache*/,
      69             :       const ArrayIndex& /*array_index*/, const ActionList /*meta*/,
      70             :       const ParallelComponent* const /*meta*/) {
      71             :     TimeDerivative::apply(make_not_null(&box));
      72             : 
      73             :     db::mutate<evolution::dg::Tags::MortarData<Dim>>(
      74             :         [](const auto mortar_data_ptr) {
      75             :           for (auto& data : *mortar_data_ptr) {
      76             :             data.second = evolution::dg::MortarDataHolder<Dim>{};
      77             :           }
      78             :         },
      79             :         make_not_null(&box));
      80             :     return {Parallel::AlgorithmExecution::Continue, std::nullopt};
      81             :   }
      82             : };
      83             : }  // namespace evolution::dg::subcell::fd::Actions