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 <memory>
       8             : #include <string>
       9             : #include <unordered_map>
      10             : #include <utility>
      11             : #include <vector>
      12             : 
      13             : #include "Domain/Creators/OptionTags.hpp"
      14             : #include "Domain/Structure/Direction.hpp"
      15             : #include "Utilities/TMPL.hpp"
      16             : 
      17             : /// \cond
      18             : template <size_t VolumeDim>
      19             : class DomainCreator;
      20             : template <size_t VolumeDim>
      21             : class ElementId;
      22             : namespace PUP {
      23             : class er;
      24             : }  // namespace PUP
      25             : /// \endcond
      26             : 
      27             : namespace evolution::dg {
      28             : /// Equal-rate region generator for regions on nonconforming block boundaries.
      29             : ///
      30             : /// For each block that borders multiple other blocks on a single
      31             : /// side, defines a region "NonconformingN", where "N" is the block
      32             : /// number, consisting of all elements on the nonconforming boundary
      33             : /// in that block and the neighboring blocks.  If the block has
      34             : /// nonconforming boundaries on both sides and is more than one
      35             : /// element across, defines two regions "NonconformingN+" and
      36             : /// "NonconformingN-".
      37             : template <size_t Dim>
      38           1 : class NonconformingEqualRateRegions {
      39             :  public:
      40           0 :   NonconformingEqualRateRegions() = default;
      41             : 
      42           0 :   using creation_tags = tmpl::list<domain::OptionTags::DomainCreator<Dim>>;
      43             : 
      44           0 :   explicit NonconformingEqualRateRegions(
      45             :       const std::unique_ptr<DomainCreator<Dim>>& domain_creator);
      46             : 
      47           0 :   std::unordered_map<std::string, size_t> regions() const;
      48             : 
      49           0 :   bool is_in_region(size_t region, const ElementId<Dim>& element_id) const;
      50             : 
      51           0 :   void pup(PUP::er& p);
      52             : 
      53             :  private:
      54           0 :   std::vector<std::string> region_names_{};
      55           0 :   std::vector<std::vector<std::pair<size_t, Direction<Dim>>>> regions_{};
      56             : };
      57             : }  // namespace evolution::dg