Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <array>
       7             : #include <cstddef>
       8             : #include <optional>
       9             : 
      10             : #include "DataStructures/DataVector.hpp"
      11             : #include "DataStructures/Tensor/EagerMath/DotProduct.hpp"
      12             : #include "DataStructures/Tensor/EagerMath/Magnitude.hpp"  // For Tags::Normalized
      13             : #include "DataStructures/Tensor/Tensor.hpp"
      14             : #include "Domain/InterfaceHelpers.hpp"
      15             : #include "Domain/Tags.hpp"
      16             : #include "Domain/TagsTimeDependent.hpp"
      17             : #include "Utilities/Gsl.hpp"
      18             : #include "Utilities/TMPL.hpp"
      19             : 
      20             : namespace domain {
      21             : namespace Tags {
      22             : /// Compute the characteristic speeds on the moving mesh given the
      23             : /// characteristic speeds if the mesh were stationary.
      24             : ///
      25             : /// \note Assumes that `typename CharSpeedsComputeTag::return_type` is a
      26             : /// `std::array<DataVector, NumberOfCharSpeeds>`
      27             : template <typename CharSpeedsComputeTag, size_t Dim>
      28           1 : struct CharSpeedCompute : CharSpeedsComputeTag::base, db::ComputeTag {
      29           0 :   using base = typename CharSpeedsComputeTag::base;
      30           0 :   using return_type = typename CharSpeedsComputeTag::return_type;
      31             : 
      32             :   template <typename... Ts, typename T, size_t NumberOfCharSpeeds>
      33           0 :   static void function(
      34             :       const gsl::not_null<std::array<T, NumberOfCharSpeeds>*> result,
      35             :       const std::optional<tnsr::I<DataVector, Dim, Frame::Inertial>>&
      36             :           grid_velocity,
      37             :       const tnsr::i<DataVector, Dim, Frame::Inertial>& unit_normal_covector,
      38             :       const Ts&... ts) {
      39             :     // Note that while the CharSpeedsComputeTag almost certainly also needs the
      40             :     // unit normal covector for computing the original characteristic speeds, we
      41             :     // don't know which of the `ts` it is, and thus we need the unit normal
      42             :     // covector to be passed explicitly.
      43             :     CharSpeedsComputeTag::function(result, ts...);
      44             :     if (grid_velocity.has_value()) {
      45             :       const Scalar<DataVector> normal_dot_velocity =
      46             :           dot_product(*grid_velocity, unit_normal_covector);
      47             :       for (size_t i = 0; i < result->size(); ++i) {
      48             :         gsl::at(*result, i) -= get(normal_dot_velocity);
      49             :       }
      50             :     }
      51             :   }
      52             : 
      53           0 :   using argument_tags =
      54             :       tmpl::push_front<typename CharSpeedsComputeTag::argument_tags,
      55             :                        MeshVelocity<Dim, Frame::Inertial>,
      56             :                        ::Tags::Normalized<UnnormalizedFaceNormal<Dim>>>;
      57           0 :   using volume_tags = get_volume_tags<CharSpeedsComputeTag>;
      58             : };
      59             : }  // namespace Tags
      60             : }  // namespace domain