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/Tag.hpp"
       9              : #include "DataStructures/DataVector.hpp"
      10              : #include "DataStructures/Tensor/EagerMath/Magnitude.hpp"
      11              : #include "DataStructures/Tensor/Tensor.hpp"
      12              : #include "PointwiseFunctions/GeneralRelativity/Tags.hpp"
      13              : #include "Utilities/Gsl.hpp"
      14              : #include "Utilities/TMPL.hpp"
      15              : 
      16            1 : namespace grmhd::GhValenciaDivClean::Tags {
      17            0 : struct LargestCharacteristicSpeed : db::SimpleTag {
      18            0 :   using type = double;
      19              : };
      20              : 
      21              : /*!
      22              :  * \brief Computes the largest magnitude of the characteristic speeds.
      23              :  *
      24              :  * \warning Assumes \f$-1\le\gamma_1\le0\f$.
      25              :  */
      26              : template <typename Frame = Frame::Inertial>
      27            1 : struct ComputeLargestCharacteristicSpeed : db::ComputeTag,
      28              :                                            LargestCharacteristicSpeed {
      29            0 :   using argument_tags =
      30              :       tmpl::list<gr::Tags::Lapse<DataVector>,
      31              :                  gr::Tags::Shift<DataVector, 3, Frame>,
      32              :                  gr::Tags::SpatialMetric<DataVector, 3, Frame>>;
      33            0 :   using return_type = double;
      34            0 :   using base = LargestCharacteristicSpeed;
      35            0 :   static void function(const gsl::not_null<double*> speed,
      36              :                        const Scalar<DataVector>& lapse,
      37              :                        const tnsr::I<DataVector, 3, Frame>& shift,
      38              :                        const tnsr::ii<DataVector, 3, Frame>& spatial_metric) {
      39              :     const auto shift_magnitude = magnitude(shift, spatial_metric);
      40              :     *speed = max(get(shift_magnitude) + get(lapse));
      41              :   }
      42              : };
      43              : }  // namespace grmhd::GhValenciaDivClean::Tags