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 "Utilities/Gsl.hpp"
      10             : 
      11             : namespace gh {
      12             : /// @{
      13             : /*!
      14             :  * \ingroup GeneralRelativityGroup
      15             :  * \brief Computes spacetime derivatives of the norm of the shift vector.
      16             :  *
      17             :  * \details The same is computed as:
      18             :  * \f{align*}
      19             :  * \partial_a (\beta^i \beta_i) =
      20             :  *     (\beta_i \partial_0 \beta^i + \beta^i \partial_0 \beta_i,
      21             :  *      \beta_i \partial_j \beta^i + \beta^i \partial_j \beta_i)
      22             :  * \f}
      23             :  */
      24             : template <typename DataType, size_t SpatialDim, typename Frame>
      25           1 : void spacetime_deriv_of_norm_of_shift(
      26             :     gsl::not_null<tnsr::a<DataType, SpatialDim, Frame>*> d4_norm_of_shift,
      27             :     const Scalar<DataType>& lapse,
      28             :     const tnsr::I<DataType, SpatialDim, Frame>& shift,
      29             :     const tnsr::ii<DataType, SpatialDim, Frame>& spatial_metric,
      30             :     const tnsr::II<DataType, SpatialDim, Frame>& inverse_spatial_metric,
      31             :     const tnsr::AA<DataType, SpatialDim, Frame>& inverse_spacetime_metric,
      32             :     const tnsr::A<DataType, SpatialDim, Frame>& spacetime_unit_normal,
      33             :     const tnsr::iaa<DataType, SpatialDim, Frame>& phi,
      34             :     const tnsr::aa<DataType, SpatialDim, Frame>& pi);
      35             : 
      36             : template <typename DataType, size_t SpatialDim, typename Frame>
      37           1 : tnsr::a<DataType, SpatialDim, Frame> spacetime_deriv_of_norm_of_shift(
      38             :     const Scalar<DataType>& lapse,
      39             :     const tnsr::I<DataType, SpatialDim, Frame>& shift,
      40             :     const tnsr::ii<DataType, SpatialDim, Frame>& spatial_metric,
      41             :     const tnsr::II<DataType, SpatialDim, Frame>& inverse_spatial_metric,
      42             :     const tnsr::AA<DataType, SpatialDim, Frame>& inverse_spacetime_metric,
      43             :     const tnsr::A<DataType, SpatialDim, Frame>& spacetime_unit_normal,
      44             :     const tnsr::iaa<DataType, SpatialDim, Frame>& phi,
      45             :     const tnsr::aa<DataType, SpatialDim, Frame>& pi);
      46             : /// @}
      47             : }  // namespace gh