Line data Source code
1 0 : // Distributed under the MIT License. 2 : // See LICENSE.txt for details. 3 : 4 : #pragma once 5 : 6 : #include <array> 7 : #include <cstddef> 8 : #include <memory> 9 : #include <utility> 10 : 11 : #include "DataStructures/DataBox/Prefixes.hpp" 12 : #include "DataStructures/Tensor/TypeAliases.hpp" 13 : #include "DataStructures/VariablesTag.hpp" 14 : #include "Domain/Structure/DirectionalIdMap.hpp" 15 : #include "Domain/Structure/Element.hpp" 16 : #include "Domain/Tags.hpp" 17 : #include "Evolution/DgSubcell/Tags/GhostDataForReconstruction.hpp" 18 : #include "Evolution/DgSubcell/Tags/Mesh.hpp" 19 : #include "Evolution/Systems/Burgers/FiniteDifference/Reconstructor.hpp" 20 : #include "Evolution/Systems/Burgers/Tags.hpp" 21 : #include "Options/String.hpp" 22 : #include "Utilities/TMPL.hpp" 23 : 24 : /// \cond 25 : class DataVector; 26 : template <size_t Dim> 27 : class Direction; 28 : template <size_t Dim> 29 : class Element; 30 : template <size_t Dim> 31 : class ElementId; 32 : template <size_t Dim> 33 : class Mesh; 34 : template <typename TagsList> 35 : class Variables; 36 : namespace gsl { 37 : template <typename> 38 : class not_null; 39 : } // namespace gsl 40 : namespace PUP { 41 : class er; 42 : } // namespace PUP 43 : namespace evolution::dg::subcell { 44 : class GhostData; 45 : } // namespace evolution::dg::subcell 46 : /// \endcond 47 : 48 : namespace Burgers::fd { 49 : /*! 50 : * \brief Monotonised central reconstruction. See 51 : * ::fd::reconstruction::monotonised_central() for details. 52 : */ 53 1 : class MonotonisedCentral : public Reconstructor { 54 : private: 55 0 : using face_vars_tags = tmpl::list< 56 : Burgers::Tags::U, 57 : ::Tags::Flux<Burgers::Tags::U, tmpl::size_t<1>, Frame::Inertial>>; 58 0 : using volume_vars_tags = tmpl::list<Burgers::Tags::U>; 59 : 60 : public: 61 0 : using options = tmpl::list<>; 62 0 : static constexpr Options::String help{ 63 : "Monotonised central reconstruction scheme."}; 64 : 65 0 : MonotonisedCentral() = default; 66 0 : MonotonisedCentral(MonotonisedCentral&&) = default; 67 0 : MonotonisedCentral& operator=(MonotonisedCentral&&) = default; 68 0 : MonotonisedCentral(const MonotonisedCentral&) = default; 69 0 : MonotonisedCentral& operator=(const MonotonisedCentral&) = default; 70 0 : ~MonotonisedCentral() override = default; 71 : 72 0 : void pup(PUP::er& p) override; 73 : 74 : /// \cond 75 : explicit MonotonisedCentral(CkMigrateMessage* msg); 76 : WRAPPED_PUPable_decl_base_template(Reconstructor, MonotonisedCentral); 77 : /// \endcond 78 : 79 0 : auto get_clone() const -> std::unique_ptr<Reconstructor> override; 80 : 81 0 : size_t ghost_zone_size() const override { return 2; } 82 : 83 0 : using reconstruction_argument_tags = 84 : tmpl::list<::Tags::Variables<volume_vars_tags>, domain::Tags::Element<1>, 85 : evolution::dg::subcell::Tags::GhostDataForReconstruction<1>, 86 : evolution::dg::subcell::Tags::Mesh<1>>; 87 : 88 0 : void reconstruct( 89 : gsl::not_null<std::array<Variables<face_vars_tags>, 1>*> 90 : vars_on_lower_face, 91 : gsl::not_null<std::array<Variables<face_vars_tags>, 1>*> 92 : vars_on_upper_face, 93 : const Variables<tmpl::list<Burgers::Tags::U>>& volume_vars, 94 : const Element<1>& element, 95 : const DirectionalIdMap<1, evolution::dg::subcell::GhostData>& ghost_data, 96 : const Mesh<1>& subcell_mesh) const; 97 : 98 0 : void reconstruct_fd_neighbor( 99 : gsl::not_null<Variables<face_vars_tags>*> vars_on_face, 100 : const Variables<volume_vars_tags>& volume_vars, const Element<1>& element, 101 : const DirectionalIdMap<1, evolution::dg::subcell::GhostData>& ghost_data, 102 : const Mesh<1>& subcell_mesh, 103 : const Direction<1> direction_to_reconstruct) const; 104 : }; 105 : 106 0 : bool operator==(const MonotonisedCentral& /*lhs*/, 107 : const MonotonisedCentral& /*rhs*/); 108 : 109 0 : bool operator!=(const MonotonisedCentral& lhs, const MonotonisedCentral& rhs); 110 : } // namespace Burgers::fd