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/DataVector.hpp"
       7              : #include "Evolution/Systems/ScalarAdvection/Tags.hpp"
       8              : #include "Utilities/Gsl.hpp"
       9              : #include "Utilities/TMPL.hpp"
      10              : 
      11              : namespace ScalarAdvection {
      12            1 : namespace Tags {
      13              : /*!
      14              :  * \brief Compute the largest characteristic speed of the ScalarAdvection system
      15              :  *
      16              :  * \f{align*}
      17              :  * \lambda_\text{maxabs} = \sqrt{v^iv_i}
      18              :  * \f}
      19              :  *
      20              :  * where \f$v^i\f$ is the velocity field.
      21              :  */
      22              : template <size_t Dim>
      23            1 : struct LargestCharacteristicSpeedCompute : LargestCharacteristicSpeed,
      24              :                                            db::ComputeTag {
      25            0 :   using argument_tags = tmpl::list<Tags::VelocityField<Dim>>;
      26            0 :   using return_type = double;
      27            0 :   using base = LargestCharacteristicSpeed;
      28            0 :   static void function(const gsl::not_null<double*> speed,
      29              :                        const tnsr::I<DataVector, Dim>& velocity_field);
      30              : };
      31              : }  // namespace Tags
      32              : }  // namespace ScalarAdvection