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/Tensor/TypeAliases.hpp"
       9              : #include "Evolution/DiscontinuousGalerkin/TimeDerivativeDecisions.hpp"
      10              : #include "Evolution/Systems/SecondOrderScalarWave/Tags.hpp"
      11              : #include "Utilities/TMPL.hpp"
      12              : 
      13              : /// \cond
      14              : namespace gsl {
      15              : template <typename T>
      16              : class not_null;
      17              : }  // namespace gsl
      18              : 
      19              : class DataVector;
      20              : /// \endcond
      21              : 
      22              : namespace SecondOrderScalarWave {
      23              : /*!
      24              :  * \brief Compute the time derivatives for the second-order scalar wave system
      25              :  */
      26              : template <size_t Dim>
      27            1 : struct TimeDerivative {
      28            0 :   using temporary_tags = tmpl::list<>;
      29              : 
      30            0 :   using argument_tags = tmpl::list<Tags::Pi>;
      31              : 
      32            0 :   static evolution::dg::TimeDerivativeDecisions<Dim> apply(
      33              :       gsl::not_null<Scalar<DataVector>*> dt_psi,
      34              :       gsl::not_null<Scalar<DataVector>*> dt_pi,
      35              : 
      36              :       // Partial derivative arguments. Listed in the system struct as
      37              :       // gradient_variables. Only the derivative of Phi enters the time
      38              :       // derivatives; the evolved variables' derivatives are in
      39              :       // gradient_variables solely for the framework's moving-mesh term.
      40              :       const tnsr::i<DataVector, Dim, Frame::Inertial>& /*d_psi*/,
      41              :       const tnsr::i<DataVector, Dim, Frame::Inertial>& /*d_pi*/,
      42              :       const tnsr::ij<DataVector, Dim, Frame::Inertial>& d_phi,
      43              : 
      44              :       // Terms list in argument_tags above
      45              :       const Scalar<DataVector>& pi);
      46              : };
      47              : }  // namespace SecondOrderScalarWave