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/Tensor/TypeAliases.hpp"
       7             : 
       8             : /// \cond
       9             : class DataVector;
      10             : namespace gsl {
      11             : template <typename>
      12             : struct not_null;
      13             : }  // namespace gsl
      14             : /// \endcond
      15             : 
      16             : namespace gr::surfaces {
      17             : /// @{
      18             : /// \ingroup SurfacesGroup
      19             : /// \brief Expansion of a `Strahlkorper`. Should be zero on apparent horizons.
      20             : ///
      21             : /// \details Implements Eq. (5) in \cite Baumgarte1996hh.  The input argument
      22             : /// `grad_normal` is the quantity returned by
      23             : /// `gr::surfaces::grad_unit_normal_one_form`, and `inverse_surface_metric`
      24             : /// is the quantity returned by `gr::surfaces::inverse_surface_metric`.
      25             : template <typename Frame>
      26           1 : void expansion(gsl::not_null<Scalar<DataVector>*> result,
      27             :                const tnsr::ii<DataVector, 3, Frame>& grad_normal,
      28             :                const tnsr::II<DataVector, 3, Frame>& inverse_surface_metric,
      29             :                const tnsr::ii<DataVector, 3, Frame>& extrinsic_curvature);
      30             : 
      31             : template <typename Frame>
      32           1 : Scalar<DataVector> expansion(
      33             :     const tnsr::ii<DataVector, 3, Frame>& grad_normal,
      34             :     const tnsr::II<DataVector, 3, Frame>& inverse_surface_metric,
      35             :     const tnsr::ii<DataVector, 3, Frame>& extrinsic_curvature);
      36             : /// @}
      37             : }  // namespace gr::surfaces