Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include "DataStructures/Tensor/TypeAliases.hpp"
       7              : #include "DataStructures/Variables.hpp"
       8              : #include "Evolution/DiscontinuousGalerkin/TimeDerivativeDecisions.hpp"
       9              : #include "Evolution/Systems/ScalarWave/Tags.hpp"
      10              : #include "NumericalAlgorithms/LinearOperators/PartialDerivatives.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 ScalarWave {
      23              : /*!
      24              :  * \brief Compute the time derivatives for scalar wave system
      25              :  */
      26              : template <size_t Dim>
      27            1 : struct TimeDerivative {
      28            0 :   using temporary_tags = tmpl::list<Tags::ConstraintGamma2>;
      29            0 :   using argument_tags =
      30              :       tmpl::list<Tags::Pi, Tags::Phi<Dim>, Tags::ConstraintGamma2>;
      31              : 
      32            0 :   static evolution::dg::TimeDerivativeDecisions<Dim> apply(
      33              :       // Time derivatives returned by reference. All the tags in the
      34              :       // variables_tag in the system struct.
      35              :       gsl::not_null<Scalar<DataVector>*> dt_psi,
      36              :       gsl::not_null<Scalar<DataVector>*> dt_pi,
      37              :       gsl::not_null<tnsr::i<DataVector, Dim, Frame::Inertial>*> dt_phi,
      38              : 
      39              :       gsl::not_null<Scalar<DataVector>*> result_gamma2,
      40              : 
      41              :       // Partial derivative arguments. Listed in the system struct as
      42              :       // gradient_variables.
      43              :       const tnsr::i<DataVector, Dim, Frame::Inertial>& d_psi,
      44              :       const tnsr::i<DataVector, Dim, Frame::Inertial>& d_pi,
      45              :       const tnsr::ij<DataVector, Dim, Frame::Inertial>& d_phi,
      46              : 
      47              :       // Terms list in argument_tags above
      48              :       const Scalar<DataVector>& pi,
      49              :       const tnsr::i<DataVector, Dim, Frame::Inertial>& phi,
      50              :       const Scalar<DataVector>& gamma2);
      51              : };
      52              : }  // namespace ScalarWave