Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <memory>
       7              : #include <pup.h>
       8              : 
       9              : #include "Utilities/Serialization/CharmPupable.hpp"
      10              : #include "Utilities/TMPL.hpp"
      11              : 
      12              : namespace NewtonianEuler::fd {
      13              : /// \cond
      14              : template <size_t Dim>
      15              : class AoWeno53Prim;
      16              : template <size_t Dim>
      17              : class MonotonisedCentralPrim;
      18              : /// \endcond
      19              : 
      20              : /*!
      21              :  * \brief The base class from which all reconstruction schemes must inherit
      22              :  *
      23              :  * Currently we have hard-coded reconstructing \f$\rho, p, v^i\f$. However, the
      24              :  * DG-subcell solver is coded generally enough that an efficient implementation
      25              :  * of reconstructing the conserved or characteristic variables is also possible.
      26              :  * It is not yet clear how much info about what is being reconstructed is needed
      27              :  * at compile time and so we currently append `Prim` to the end of the
      28              :  * reconstruction schemes to clarify that they are reconstructing the primitive
      29              :  * variables. Ideally the choice of what variables to reconstruct can be made by
      30              :  * a runtime argument to the individual reconstruction schemes.
      31              :  */
      32              : template <size_t Dim>
      33            1 : class Reconstructor : public PUP::able {
      34              :  public:
      35            0 :   Reconstructor() = default;
      36            0 :   Reconstructor(const Reconstructor&) = default;
      37            0 :   Reconstructor& operator=(const Reconstructor&) = default;
      38            0 :   Reconstructor(Reconstructor&&) = default;
      39            0 :   Reconstructor& operator=(Reconstructor&&) = default;
      40            0 :   ~Reconstructor() override = default;
      41              : 
      42              :   /// \cond
      43              :   explicit Reconstructor(CkMigrateMessage* msg);
      44              :   WRAPPED_PUPable_abstract(Reconstructor<Dim>);  // NOLINT
      45              :   /// \endcond
      46              : 
      47            0 :   using creatable_classes =
      48              :       tmpl::list<AoWeno53Prim<Dim>, MonotonisedCentralPrim<Dim>>;
      49              : 
      50            0 :   virtual std::unique_ptr<Reconstructor<Dim>> get_clone() const = 0;
      51              : 
      52            0 :   virtual size_t ghost_zone_size() const = 0;
      53              : 
      54            0 :   void pup(PUP::er& p) override;
      55              : };
      56              : }  // namespace NewtonianEuler::fd