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/ComplexDataVector.hpp"
       9              : #include "DataStructures/ComplexModalVector.hpp"
      10              : #include "DataStructures/SpinWeighted.hpp"
      11              : #include "DataStructures/Tensor/TypeAliases.hpp"
      12              : #include "Evolution/Systems/Cce/BoundaryDataTags.hpp"
      13              : 
      14              : /// \cond
      15              : class DataVector;
      16              : namespace gsl {
      17              : template <class T>
      18              : class not_null;
      19              : }  // namespace gsl
      20              : /// \endcond
      21              : 
      22              : namespace Cce {
      23              : /*
      24              :  * \brief Compute \f$\gamma_{i j}\f$, \f$\gamma^{i j}\f$,
      25              :  * \f$\partial_i \gamma_{j k}\f$, and
      26              :  * \f$\partial_t g_{i j}\f$ from input libsharp-compatible modal spatial
      27              :  * metric quantities.
      28              :  *
      29              :  * \details This function will apply a correction factor associated with a SpEC
      30              :  * bug.
      31              :  */
      32            0 : void cartesian_spatial_metric_and_derivatives_from_unnormalized_spec_modes(
      33              :     gsl::not_null<tnsr::ii<DataVector, 3>*> cartesian_spatial_metric,
      34              :     gsl::not_null<tnsr::II<DataVector, 3>*> inverse_cartesian_spatial_metric,
      35              :     gsl::not_null<tnsr::ijj<DataVector, 3>*> d_cartesian_spatial_metric,
      36              :     gsl::not_null<tnsr::ii<DataVector, 3>*> dt_cartesian_spatial_metric,
      37              :     gsl::not_null<Scalar<SpinWeighted<ComplexModalVector, 0>>*>
      38              :         interpolation_modal_buffer,
      39              :     gsl::not_null<Scalar<SpinWeighted<ComplexDataVector, 0>>*>
      40              :         interpolation_buffer,
      41              :     gsl::not_null<Scalar<SpinWeighted<ComplexDataVector, 1>>*> eth_buffer,
      42              :     gsl::not_null<Scalar<DataVector>*> radial_correction_factor,
      43              :     const tnsr::ii<ComplexModalVector, 3>& spatial_metric_coefficients,
      44              :     const tnsr::ii<ComplexModalVector, 3>& dr_spatial_metric_coefficients,
      45              :     const tnsr::ii<ComplexModalVector, 3>& dt_spatial_metric_coefficients,
      46              :     const CartesianiSphericalJ& inverse_cartesian_to_spherical_jacobian,
      47              :     const tnsr::I<DataVector, 3>& unit_cartesian_coords, size_t l_max);
      48              : 
      49              : /*!
      50              :  * \brief Compute \f$\beta^{i}\f$, \f$\partial_i \beta^{j}\f$, and
      51              :  * \f$\partial_t \beta^i\f$ from input libsharp-compatible modal spatial
      52              :  * metric quantities.
      53              :  *
      54              :  * \details This function will apply a correction factor associated with a SpEC
      55              :  * bug.
      56              :  */
      57            1 : void cartesian_shift_and_derivatives_from_unnormalized_spec_modes(
      58              :     gsl::not_null<tnsr::I<DataVector, 3>*> cartesian_shift,
      59              :     gsl::not_null<tnsr::iJ<DataVector, 3>*> d_cartesian_shift,
      60              :     gsl::not_null<tnsr::I<DataVector, 3>*> dt_cartesian_shift,
      61              :     gsl::not_null<Scalar<SpinWeighted<ComplexModalVector, 0>>*>
      62              :         interpolation_modal_buffer,
      63              :     gsl::not_null<Scalar<SpinWeighted<ComplexDataVector, 0>>*>
      64              :         interpolation_buffer,
      65              :     gsl::not_null<Scalar<SpinWeighted<ComplexDataVector, 1>>*> eth_buffer,
      66              :     const tnsr::I<ComplexModalVector, 3>& shift_coefficients,
      67              :     const tnsr::I<ComplexModalVector, 3>& dr_shift_coefficients,
      68              :     const tnsr::I<ComplexModalVector, 3>& dt_shift_coefficients,
      69              :     const CartesianiSphericalJ& inverse_cartesian_to_spherical_jacobian,
      70              :     const Scalar<DataVector>& radial_derivative_correction_factor,
      71              :     size_t l_max);
      72              : 
      73              : /*!
      74              :  * \brief Compute \f$\alpha\f$, \f$\partial_i \alpha\f$, and
      75              :  * \f$\partial_t \beta^i\f$ from input libsharp-compatible modal spatial
      76              :  * metric quantities.
      77              :  *
      78              :  * \details This function will apply a correction factor associated with a SpEC
      79              :  * bug.
      80              :  */
      81            1 : void cartesian_lapse_and_derivatives_from_unnormalized_spec_modes(
      82              :     gsl::not_null<Scalar<DataVector>*> cartesian_lapse,
      83              :     gsl::not_null<tnsr::i<DataVector, 3>*> d_cartesian_lapse,
      84              :     gsl::not_null<Scalar<DataVector>*> dt_cartesian_lapse,
      85              :     gsl::not_null<Scalar<SpinWeighted<ComplexModalVector, 0>>*>
      86              :         interpolation_modal_buffer,
      87              :     gsl::not_null<Scalar<SpinWeighted<ComplexDataVector, 0>>*>
      88              :         interpolation_buffer,
      89              :     gsl::not_null<Scalar<SpinWeighted<ComplexDataVector, 1>>*> eth_buffer,
      90              :     const Scalar<ComplexModalVector>& lapse_coefficients,
      91              :     const Scalar<ComplexModalVector>& dr_lapse_coefficients,
      92              :     const Scalar<ComplexModalVector>& dt_lapse_coefficients,
      93              :     const CartesianiSphericalJ& inverse_cartesian_to_spherical_jacobian,
      94              :     const Scalar<DataVector>& radial_derivative_correction_factor,
      95              :     size_t l_max);
      96              : }  // namespace Cce