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              : 
       8              : #include "DataStructures/DataBox/Prefixes.hpp"
       9              : #include "DataStructures/DataVector.hpp"
      10              : #include "DataStructures/Tensor/TypeAliases.hpp"
      11              : #include "Evolution/DiscontinuousGalerkin/TimeDerivativeDecisions.hpp"
      12              : #include "Evolution/Systems/ScalarAdvection/Tags.hpp"
      13              : #include "Utilities/Gsl.hpp"
      14              : #include "Utilities/TMPL.hpp"
      15              : 
      16              : namespace ScalarAdvection {
      17              : /*!
      18              :  * \brief Computes the time derivative terms needed for the ScalarAdvection
      19              :  * system, which are just the fluxes.
      20              :  */
      21              : template <size_t Dim>
      22            1 : struct TimeDerivativeTerms {
      23            0 :   using temporary_tags = tmpl::list<Tags::VelocityField<Dim>>;
      24            0 :   using argument_tags = tmpl::list<Tags::U, Tags::VelocityField<Dim>>;
      25              : 
      26            0 :   static evolution::dg::TimeDerivativeDecisions<Dim> apply(
      27              :       // Time derivatives returned by reference. No source terms or
      28              :       // nonconservative products, so not used. All the tags in the
      29              :       // variables_tag in the system struct.
      30              :       gsl::not_null<Scalar<DataVector>*> /*non_flux_terms_dt_vars*/,
      31              : 
      32              :       // Fluxes returned by reference. Listed in the system struct as
      33              :       // flux_variables.
      34              :       gsl::not_null<tnsr::I<DataVector, Dim>*> flux,
      35              : 
      36              :       // Temporary tags
      37              :       gsl::not_null<tnsr::I<DataVector, Dim>*> temp_velocity_field,
      38              : 
      39              :       // Arguments listed in argument_tags above
      40              :       const Scalar<DataVector>& u,
      41              :       const tnsr::I<DataVector, Dim>& velocity_field);
      42              : };
      43              : }  // namespace ScalarAdvection