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/DataBox/Tag.hpp"
       9              : #include "DataStructures/Tensor/TypeAliases.hpp"
      10              : #include "Utilities/TMPL.hpp"
      11              : 
      12              : /// \cond
      13              : class ComplexDataVector;
      14              : class DataVector;
      15              : /// \endcond
      16              : 
      17              : namespace Cce {
      18            0 : namespace Frame {
      19              : /// The frame for the spherical metric in which the radial coordinate is an
      20              : /// affine parameter along outward-pointing null geodesics.
      21            1 : struct RadialNull {};
      22              : }  // namespace Frame
      23              : 
      24              : // tensor aliases for brevity
      25            0 : using SphericaliCartesianJ = Tensor<
      26              :     DataVector, tmpl::integral_list<std::int32_t, 2, 1>,
      27              :     index_list<SpatialIndex<3, UpLo::Lo, ::Frame::Spherical<::Frame::Inertial>>,
      28              :                SpatialIndex<3, UpLo::Lo, ::Frame::Inertial>>>;
      29              : 
      30            0 : using CartesianiSphericalJ =
      31              :     Tensor<DataVector, tmpl::integral_list<std::int32_t, 2, 1>,
      32              :            index_list<SpatialIndex<3, UpLo::Lo, ::Frame::Inertial>,
      33              :                       SpatialIndex<3, UpLo::Up,
      34              :                                    ::Frame::Spherical<::Frame::Inertial>>>>;
      35              : 
      36            0 : using AngulariCartesianA = Tensor<
      37              :   DataVector, tmpl::integral_list<std::int32_t, 2, 1>,
      38              :   index_list<SpatialIndex<2, UpLo::Lo, ::Frame::Spherical<::Frame::Inertial>>,
      39              :              SpacetimeIndex<3, UpLo::Lo, ::Frame::Inertial>>>;
      40              : 
      41            0 : using SphericaliCartesianjj = Tensor<
      42              :     DataVector, tmpl::integral_list<std::int32_t, 2, 1, 1>,
      43              :     index_list<SpatialIndex<3, UpLo::Lo, ::Frame::Spherical<::Frame::Inertial>>,
      44              :                SpatialIndex<3, UpLo::Lo, ::Frame::Inertial>,
      45              :                SpatialIndex<3, UpLo::Lo, ::Frame::Inertial>>>;
      46              : 
      47              : namespace Tags {
      48              : namespace detail {
      49              : // this provides a set of tags for the purposes of allocating once in the entire
      50              : // Boundary data computation; these tags are currently not used outside
      51              : // intermediate steps of the procedure in `BoundaryData.hpp`
      52              : 
      53              : struct CosPhi : db::SimpleTag {
      54              :   using type = Scalar<DataVector>;
      55              : };
      56              : 
      57              : struct CosTheta : db::SimpleTag {
      58              :   using type = Scalar<DataVector>;
      59              : };
      60              : 
      61              : struct SinPhi : db::SimpleTag {
      62              :   using type = Scalar<DataVector>;
      63              : };
      64              : 
      65              : struct SinTheta : db::SimpleTag {
      66              :   using type = Scalar<DataVector>;
      67              : };
      68              : 
      69              : struct CartesianCoordinates : db::SimpleTag {
      70              :   using type = tnsr::I<DataVector, 3>;
      71              : };
      72              : 
      73              : struct CartesianToSphericalJacobian : db::SimpleTag {
      74              :   using type = SphericaliCartesianJ;
      75              : };
      76              : 
      77              : struct InverseCartesianToSphericalJacobian : db::SimpleTag {
      78              :   using type = CartesianiSphericalJ;
      79              : };
      80              : 
      81              : struct WorldtubeNormal : db::SimpleTag {
      82              :   using type = tnsr::I<DataVector, 3>;
      83              : };
      84              : 
      85              : struct UpDyad : db::SimpleTag {
      86              :   using type = tnsr::I<ComplexDataVector, 2, Frame::RadialNull>;
      87              : };
      88              : 
      89              : struct DownDyad : db::SimpleTag {
      90              :   using type = tnsr::i<ComplexDataVector, 2, Frame::RadialNull>;
      91              : };
      92              : 
      93              : struct RealBondiR : db::SimpleTag {
      94              :   using type = Scalar<DataVector>;
      95              : };
      96              : 
      97              : struct AngularDNullL : db::SimpleTag {
      98              :   using type = AngulariCartesianA;
      99              : };
     100              : 
     101              : struct NullL : db::SimpleTag {
     102              :   using type = tnsr::A<DataVector, 3>;
     103              : };
     104              : 
     105              : template <typename Tag>
     106              : struct DLambda : db::SimpleTag {
     107              :   using type = typename Tag::type;
     108              : };
     109              : 
     110              : }  // namespace detail
     111              : }  // namespace Tags
     112              : }  // namespace Cce