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/DataVector.hpp"
       9              : #include "DataStructures/Tensor/Tensor.hpp"
      10              : #include "Domain/Tags.hpp"
      11              : #include "Evolution/Systems/SecondOrderScalarWave/Tags.hpp"
      12              : #include "NumericalAlgorithms/LinearOperators/PartialDerivatives.hpp"
      13              : #include "NumericalAlgorithms/Spectral/Mesh.hpp"
      14              : #include "Utilities/Gsl.hpp"
      15              : #include "Utilities/TMPL.hpp"
      16              : 
      17              : namespace SecondOrderScalarWave {
      18              : /*!
      19              :  * \brief Computes the auxiliary variable \f$\Phi_i = \partial_i \Psi\f$
      20              :  *
      21              :  * \details In the LDG scheme, the auxiliary variable \f$\Phi_i\f$ is not
      22              :  * evolved in time but is recomputed from \f$\Psi\f$ at each substep.
      23              :  * After this mutator runs, the boundary correction from
      24              :  * `ApplyAuxiliaryBoundaryCorrectionsToVariables` will correct \f$\Phi_i\f$
      25              :  * at element interfaces.
      26              :  */
      27              : template <size_t Dim>
      28            1 : struct UpdateAuxiliaryVariables {
      29            0 :   using return_tags = tmpl::list<Tags::Phi<Dim>>;
      30            0 :   using argument_tags =
      31              :       tmpl::list<Tags::Psi, domain::Tags::Mesh<Dim>,
      32              :                  domain::Tags::InverseJacobian<Dim, Frame::ElementLogical,
      33              :                                                Frame::Inertial>>;
      34              : 
      35            0 :   static void apply(
      36              :       const gsl::not_null<tnsr::i<DataVector, Dim, Frame::Inertial>*> phi,
      37              :       const Scalar<DataVector>& psi, const Mesh<Dim>& mesh,
      38              :       const InverseJacobian<DataVector, Dim, Frame::ElementLogical,
      39              :                             Frame::Inertial>& inverse_jacobian) {
      40              :     partial_derivative(phi, psi, mesh, inverse_jacobian);
      41              :   }
      42              : };
      43              : }  // namespace SecondOrderScalarWave