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 <iterator>
       8             : #include <memory>
       9             : #include <unordered_map>
      10             : #include <utility>
      11             : 
      12             : #include "DataStructures/TaggedTuple.hpp"
      13             : #include "Domain/Structure/Element.hpp"
      14             : #include "Domain/Structure/ElementId.hpp"
      15             : #include "NumericalAlgorithms/LinearOperators/Filters/Filter.hpp"
      16             : #include "NumericalAlgorithms/LinearOperators/Filters/Tag.hpp"
      17             : #include "NumericalAlgorithms/Spectral/Mesh.hpp"
      18             : #include "ParallelAlgorithms/Amr/Protocols/Projector.hpp"
      19             : #include "Utilities/ErrorHandling/Error.hpp"
      20             : #include "Utilities/Gsl.hpp"
      21             : #include "Utilities/TMPL.hpp"
      22             : 
      23             : namespace evolution::dg::Initialization {
      24             : 
      25             : /// \brief AMR projector for `Filters::Tags::SpectralFilter<Dim, TagList>`.
      26             : ///
      27             : /// \details
      28             : /// - For p-refinement: leaves the filter unchanged.
      29             : /// - For h-refinement (splitting): clones the parent's filter for the child.
      30             : /// - For h-coarsening (joining): clones the first child's filter for the
      31             : ///   parent (all siblings share the same block and mesh basis/quadrature, so
      32             : ///   they carry identical filters).
      33             : template <size_t Dim, typename TagList>
      34           1 : struct ProjectSpectralFilters : tt::ConformsTo<amr::protocols::Projector> {
      35           0 :   using return_tags = tmpl::list<Filters::Tags::SpectralFilter<Dim, TagList>>;
      36           0 :   using argument_tags = tmpl::list<>;
      37             : 
      38             :   // p-refinement: leave the filter unchanged.
      39           0 :   static void apply(
      40             :       const gsl::not_null<
      41             :           std::unique_ptr<Filters::Filter<Dim, TagList>>*> /*filter*/,
      42             :       const std::pair<Mesh<Dim>, Element<Dim>>& /*old_mesh_and_element*/) {}
      43             : 
      44             :   // Splitting: clone the parent's filter for the new child.
      45             :   template <typename... ParentTags>
      46           0 :   static void apply(
      47             :       const gsl::not_null<std::unique_ptr<Filters::Filter<Dim, TagList>>*>
      48             :           filter,
      49             :       const tuples::TaggedTuple<ParentTags...>& parent_items) {
      50             :     *filter =
      51             :         tuples::get<Filters::Tags::SpectralFilter<Dim, TagList>>(parent_items)
      52             :             ->get_clone();
      53             :   }
      54             : 
      55             :   // Joining: clone the first child's filter for the parent.
      56             :   template <typename... ChildrenTags>
      57           0 :   static void apply(
      58             :       const gsl::not_null<std::unique_ptr<Filters::Filter<Dim, TagList>>*>
      59             :           filter,
      60             :       const std::unordered_map<ElementId<Dim>,
      61             :                                tuples::TaggedTuple<ChildrenTags...>>&
      62             :           children_items) {
      63             :     const auto first_child_items = children_items.begin();
      64             :     const auto& first_filter =
      65             :         tuples::get<Filters::Tags::SpectralFilter<Dim, TagList>>(
      66             :             first_child_items->second);
      67             :     for (auto it = std::next(first_child_items); it != children_items.end();
      68             :          ++it) {
      69             :       const auto& other_filter =
      70             :           tuples::get<Filters::Tags::SpectralFilter<Dim, TagList>>(it->second);
      71             :       if (not first_filter->is_equal(*other_filter)) {
      72             :         ERROR("Children do not agree on all items!");
      73             :       }
      74             :     }
      75             :     *filter = first_filter->get_clone();
      76             :   }
      77             : };
      78             : 
      79             : }  // namespace evolution::dg::Initialization