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             : 
       9             : #include "DataStructures/DataBox/Tag.hpp"
      10             : #include "Evolution/DgSubcell/SubcellOptions.hpp"
      11             : #include "Evolution/DgSubcell/Tags/SubcellOptions.hpp"
      12             : #include "Evolution/DgSubcell/Tags/SubcellSolver.hpp"
      13             : #include "Evolution/Systems/GrMhd/ValenciaDivClean/FiniteDifference/Reconstructor.hpp"
      14             : #include "Options/String.hpp"
      15             : #include "Utilities/ErrorHandling/Error.hpp"
      16             : #include "Utilities/TMPL.hpp"
      17             : 
      18             : namespace grmhd::ValenciaDivClean::fd {
      19             : /// Option tags for reconstruction
      20           1 : namespace OptionTags {
      21             : /// \brief Option tag for the reconstructor
      22           1 : struct Reconstructor {
      23           0 :   using type = std::unique_ptr<fd::Reconstructor>;
      24             : 
      25           0 :   static constexpr Options::String help = {"The reconstruction scheme to use."};
      26           0 :   using group = evolution::dg::subcell::OptionTags::SubcellSolverGroup;
      27             : };
      28             : }  // namespace OptionTags
      29             : 
      30             : /// %Tags for reconstruction
      31           1 : namespace Tags {
      32             : /// \brief Tag for the reconstructor
      33           1 : struct Reconstructor : db::SimpleTag {
      34           0 :   using type = std::unique_ptr<fd::Reconstructor>;
      35           0 :   using option_tags =
      36             :       tmpl::list<OptionTags::Reconstructor,
      37             :                  ::evolution::dg::subcell::OptionTags::SubcellOptions>;
      38             : 
      39           0 :   static constexpr bool pass_metavariables = false;
      40           0 :   static type create_from_options(
      41             :       const type& reconstructor,
      42             :       const ::evolution::dg::subcell::SubcellOptions& subcell_options) {
      43             :     if (static_cast<int>(subcell_options.finite_difference_derivative_order()) <
      44             :             0 and
      45             :         not reconstructor->supports_adaptive_order()) {
      46             :       ERROR_NO_TRACE(
      47             :           "Cannot use adaptive finite difference derivative order with "
      48             :           "specified reconstructor.");
      49             :     }
      50             :     if ((static_cast<int>(
      51             :              subcell_options.finite_difference_derivative_order()) /
      52             :              2 -
      53             :          1) > static_cast<int>(reconstructor->ghost_zone_size())) {
      54             :       ERROR_NO_TRACE(
      55             :           "The derivative order chosen ("
      56             :           << subcell_options.finite_difference_derivative_order()
      57             :           << ") requires more ghost zones ("
      58             :           << (static_cast<int>(
      59             :                   subcell_options.finite_difference_derivative_order()) /
      60             :                   2 -
      61             :               1)
      62             :           << ") than the reconstructor sends, "
      63             :           << reconstructor->ghost_zone_size());
      64             :     }
      65             :     return reconstructor->get_clone();
      66             :   }
      67             : };
      68             : }  // namespace Tags
      69             : }  // namespace NewtonianEuler::fd