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              : #include <pup_stl.h>
       9              : 
      10              : #include "DataStructures/Tensor/Tensor.hpp"
      11              : #include "Utilities/Serialization/CharmPupable.hpp"
      12              : 
      13              : /// \cond
      14              : class DataVector;
      15              : namespace gsl {
      16              : template <class T>
      17              : class not_null;
      18              : }  // namespace gsl
      19              : namespace EquationsOfState {
      20              : template <bool IsRelativistic, size_t ThermodynamicDim>
      21              : class EquationOfState;
      22              : }  // namespace EquationsOfState
      23              : /// \endcond
      24              : 
      25              : namespace NewtonianEuler {
      26              : /*!
      27              :  * Holds classes implementing sources for the Newtonian Euler system.
      28              :  */
      29              : namespace Sources {
      30              : /*!
      31              :  * \brief Source terms base class.
      32              :  */
      33              : template <size_t Dim>
      34            1 : class Source : public PUP::able {
      35              :  protected:
      36            0 :   Source() = default;
      37              : 
      38              :  public:
      39            0 :   ~Source() override = default;
      40              : 
      41              :   /// \cond
      42              :   explicit Source(CkMigrateMessage* msg) : PUP::able(msg) {}
      43              :   WRAPPED_PUPable_abstract(Source);
      44              :   /// \endcond
      45              : 
      46            0 :   virtual auto get_clone() const -> std::unique_ptr<Source> = 0;
      47              : 
      48            0 :   virtual void operator()(
      49              :       gsl::not_null<Scalar<DataVector>*> source_mass_density_cons,
      50              :       gsl::not_null<tnsr::I<DataVector, Dim>*> source_momentum_density,
      51              :       gsl::not_null<Scalar<DataVector>*> source_energy_density,
      52              :       const Scalar<DataVector>& mass_density_cons,
      53              :       const tnsr::I<DataVector, Dim>& momentum_density,
      54              :       const Scalar<DataVector>& energy_density,
      55              :       const tnsr::I<DataVector, Dim>& velocity,
      56              :       const Scalar<DataVector>& pressure,
      57              :       const Scalar<DataVector>& specific_internal_energy,
      58              :       const EquationsOfState::EquationOfState<false, 2>& eos,
      59              :       const tnsr::I<DataVector, Dim>& coords, double time) const = 0;
      60              : };
      61              : }  // namespace Sources
      62              : }  // namespace NewtonianEuler