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 "Domain/Tags.hpp"
      10             : #include "Evolution/Systems/ScalarAdvection/Tags.hpp"
      11             : #include "Utilities/Gsl.hpp"
      12             : 
      13             : /// \cond
      14             : class DataVector;
      15             : namespace gsl {
      16             : template <typename T>
      17             : class not_null;
      18             : }  // namespace gsl
      19             : /// \endcond
      20             : 
      21             : namespace ScalarAdvection::Tags {
      22             : /*!
      23             :  * \brief Compute the advection velocity field \f$v\f$ of the ScalarAdvection
      24             :  * system
      25             :  *
      26             :  * - For 1D problem, \f$v(x) = 1.0\f$
      27             :  * - For 2D problem, \f$v(x,y) = (0.5-y,-0.5+x)\f$
      28             :  *
      29             :  */
      30             : template <size_t Dim>
      31           1 : struct VelocityFieldCompute : VelocityField<Dim>, db::ComputeTag {
      32           0 :   using argument_tags =
      33             :       tmpl::list<domain::Tags::Coordinates<Dim, Frame::Inertial>>;
      34           0 :   using return_type = tnsr::I<DataVector, Dim>;
      35           0 :   using base = VelocityField<Dim>;
      36             : 
      37           0 :   static void function(
      38             :       const gsl::not_null<tnsr::I<DataVector, Dim>*> velocity_field,
      39             :       const tnsr::I<DataVector, Dim, Frame::Inertial>& inertial_coords);
      40             : };
      41             : }  // namespace ScalarAdvection::Tags