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              : /// \cond
       9              : class DataVector;
      10              : /// \endcond
      11              : 
      12              : #include "DataStructures/DataBox/Tag.hpp"
      13              : #include "DataStructures/Tensor/EagerMath/DeterminantAndInverse.hpp"
      14              : #include "DataStructures/Tensor/EagerMath/RaiseOrLowerIndex.hpp"
      15              : #include "DataStructures/Tensor/TypeAliases.hpp"
      16              : #include "Evolution/Systems/GeneralizedHarmonic/Tags.hpp"
      17              : #include "PointwiseFunctions/GeneralRelativity/Christoffel.hpp"
      18              : #include "PointwiseFunctions/GeneralRelativity/ExtrinsicCurvature.hpp"
      19              : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/DerivSpatialMetric.hpp"
      20              : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/ExtrinsicCurvature.hpp"
      21              : #include "PointwiseFunctions/GeneralRelativity/Lapse.hpp"
      22              : #include "PointwiseFunctions/GeneralRelativity/Shift.hpp"
      23              : #include "PointwiseFunctions/GeneralRelativity/SpacetimeNormalVector.hpp"
      24              : #include "PointwiseFunctions/GeneralRelativity/SpatialMetric.hpp"
      25              : #include "PointwiseFunctions/GeneralRelativity/Tags.hpp"
      26              : #include "Utilities/Gsl.hpp"
      27              : #include "Utilities/SetNumberOfGridPoints.hpp"
      28              : #include "Utilities/TMPL.hpp"
      29              : 
      30            1 : namespace gr::surfaces::Tags {
      31              : /// @{
      32              : /// These ComputeItems are different from those used in
      33              : /// GeneralizedHarmonic evolution because these live only on the
      34              : /// intrp::Actions::ApparentHorizon DataBox, not in the volume
      35              : /// DataBox.  And these ComputeItems can do fewer allocations than the
      36              : /// volume ones, because (for example) Lapse, SpaceTimeNormalVector,
      37              : /// etc.  can be inlined instead of being allocated as a separate
      38              : /// ComputeItem.
      39              : template <size_t Dim, typename Frame>
      40            1 : struct InverseSpatialMetricCompute
      41              :     : gr::Tags::InverseSpatialMetric<DataVector, Dim, Frame>,
      42              :       db::ComputeTag {
      43            0 :   using return_type = tnsr::II<DataVector, Dim, Frame>;
      44            0 :   static void function(
      45              :       const gsl::not_null<tnsr::II<DataVector, Dim, Frame>*> result,
      46              :       const tnsr::aa<DataVector, Dim, Frame>& psi) {
      47              :     *result = determinant_and_inverse(gr::spatial_metric(psi)).second;
      48              :   };
      49            0 :   using argument_tags =
      50              :       tmpl::list<gr::Tags::SpacetimeMetric<DataVector, Dim, Frame>>;
      51            0 :   using base = gr::Tags::InverseSpatialMetric<DataVector, Dim, Frame>;
      52              : };
      53              : template <size_t Dim, typename Frame>
      54            0 : struct ExtrinsicCurvatureCompute
      55              :     : gr::Tags::ExtrinsicCurvature<DataVector, Dim, Frame>,
      56              :       db::ComputeTag {
      57            0 :   using return_type = tnsr::ii<DataVector, Dim, Frame>;
      58            0 :   static void function(
      59              :       const gsl::not_null<tnsr::ii<DataVector, Dim, Frame>*> result,
      60              :       const tnsr::aa<DataVector, Dim, Frame>& psi,
      61              :       const tnsr::aa<DataVector, Dim, Frame>& pi,
      62              :       const tnsr::iaa<DataVector, Dim, Frame>& phi,
      63              :       const tnsr::II<DataVector, Dim, Frame>& inv_g) {
      64              :     const auto shift = gr::shift(psi, inv_g);
      65              :     set_number_of_grid_points(result, psi);
      66              :     gh::extrinsic_curvature(
      67              :         result, gr::spacetime_normal_vector(gr::lapse(shift, psi), shift), pi,
      68              :         phi);
      69              :   }
      70            0 :   using argument_tags =
      71              :       tmpl::list<gr::Tags::SpacetimeMetric<DataVector, Dim, Frame>,
      72              :                  gh::Tags::Pi<DataVector, Dim, Frame>,
      73              :                  gh::Tags::Phi<DataVector, Dim, Frame>,
      74              :                  gr::Tags::InverseSpatialMetric<DataVector, Dim, Frame>>;
      75            0 :   using base = gr::Tags::ExtrinsicCurvature<DataVector, Dim, Frame>;
      76              : };
      77              : template <size_t Dim, typename Frame>
      78            0 : struct SpatialChristoffelSecondKindCompute
      79              :     : ::gr::Tags::SpatialChristoffelSecondKind<DataVector, Dim, Frame>,
      80              :       db::ComputeTag {
      81            0 :   using return_type = tnsr::Ijj<DataVector, Dim, Frame>;
      82            0 :   static void function(
      83              :       const gsl::not_null<tnsr::Ijj<DataVector, Dim, Frame>*> result,
      84              :       const tnsr::iaa<DataVector, Dim, Frame>& phi,
      85              :       const tnsr::II<DataVector, Dim, Frame>& inv_g) {
      86              :     set_number_of_grid_points(result, phi);
      87              :     raise_or_lower_first_index(
      88              :         result, gr::christoffel_first_kind(gh::deriv_spatial_metric(phi)),
      89              :         inv_g);
      90              :   }
      91            0 :   using argument_tags =
      92              :       tmpl::list<gh::Tags::Phi<DataVector, Dim, Frame>,
      93              :                  gr::Tags::InverseSpatialMetric<DataVector, Dim, Frame>>;
      94            0 :   using base = ::gr::Tags::SpatialChristoffelSecondKind<DataVector, Dim, Frame>;
      95              : };
      96              : /// @}
      97              : }  // namespace gr::surfaces::Tags