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              : 
       8              : #include "DataStructures/DataBox/DataBox.hpp"
       9              : #include "DataStructures/Tensor/EagerMath/Norms.hpp"
      10              : #include "DataStructures/Variables.hpp"
      11              : #include "Evolution/Systems/ScalarWave/Constraints.hpp"
      12              : #include "Evolution/Systems/ScalarWave/System.hpp"
      13              : #include "Evolution/Systems/ScalarWave/Tags.hpp"
      14              : #include "NumericalAlgorithms/LinearOperators/PartialDerivatives.hpp"
      15              : #include "Parallel/AlgorithmExecution.hpp"
      16              : #include "Parallel/GlobalCache.hpp"
      17              : #include "ParallelAlgorithms/Initialization/MutateAssign.hpp"
      18              : #include "Utilities/ErrorHandling/Assert.hpp"
      19              : 
      20              : namespace ScalarWave {
      21            0 : namespace Actions {
      22              : /// \ingroup InitializationGroup
      23              : /// \brief Initialize items related to constraints of the ScalarWave system
      24              : ///
      25              : /// We add both constraints and the constraint damping parameter to the
      26              : /// evolution databox.
      27              : ///
      28              : /// DataBox changes:
      29              : /// - Adds:
      30              : ///   * `ScalarWave::Tags::ConstraintGamma2`
      31              : /// - Removes: nothing
      32              : /// - Modifies: nothing
      33              : template <size_t Dim>
      34            1 : struct InitializeConstraints {
      35            0 :   using simple_tags = tmpl::list<ScalarWave::Tags::ConstraintGamma2>;
      36              : 
      37            0 :   using compute_tags = tmpl::list<>;
      38              : 
      39              :   template <typename DbTagsList, typename... InboxTags, typename Metavariables,
      40              :             typename ArrayIndex, typename ActionList,
      41              :             typename ParallelComponent>
      42            0 :   static Parallel::iterable_action_return_t apply(
      43              :       db::DataBox<DbTagsList>& box,
      44              :       const tuples::TaggedTuple<InboxTags...>& /*inboxes*/,
      45              :       const Parallel::GlobalCache<Metavariables>& /*cache*/,
      46              :       const ArrayIndex& /*array_index*/, const ActionList /*meta*/,
      47              :       const ParallelComponent* const /*meta*/) {
      48              :     const auto& mesh = db::get<domain::Tags::Mesh<Dim>>(box);
      49              :     Scalar<DataVector> gamma_2{mesh.number_of_grid_points(), 0.};
      50              : 
      51              :     Initialization::mutate_assign<simple_tags>(make_not_null(&box),
      52              :                                                std::move(gamma_2));
      53              : 
      54              :     return {Parallel::AlgorithmExecution::Continue, std::nullopt};
      55              :   }
      56              : };
      57              : }  // namespace Actions
      58              : }  // namespace ScalarWave