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 Ccz4 {
      12             : /// @{
      13             : /*!
      14             :  * \brief Computes the spatial part of the Z4 constraint
      15             :  *
      16             :  * \details Computes the constraint as:
      17             :  *
      18             :  * \f{align}
      19             :  *     Z_i &= \frac{1}{2} \tilde{\gamma}_{ij} \left(
      20             :  *         \hat{\Gamma}^j - \tilde{\Gamma}^j\right)
      21             :  * \f}
      22             :  *
      23             :  * where \f$\tilde{\gamma}_{ij}\f$ is the conformal spatial metric defined by
      24             :  * `Ccz4::Tags::ConformalMetric` and
      25             :  * \f$\left(\hat{\Gamma}^i - \tilde{\Gamma}^i\right)\f$ is the CCZ4 temporary
      26             :  * expression defined by
      27             :  * `Ccz4::Tags::GammaHatMinusContractedConformalChristoffel`.
      28             :  */
      29             : template <typename DataType, size_t Dim, typename Frame>
      30           1 : void spatial_z4_constraint(
      31             :     const gsl::not_null<tnsr::i<DataType, Dim, Frame>*> result,
      32             :     const tnsr::ii<DataType, Dim, Frame>& conformal_spatial_metric,
      33             :     const tnsr::I<DataType, Dim, Frame>&
      34             :         gamma_hat_minus_contracted_conformal_christoffel);
      35             : 
      36             : template <typename DataType, size_t Dim, typename Frame>
      37           1 : tnsr::i<DataType, Dim, Frame> spatial_z4_constraint(
      38             :     const tnsr::ii<DataType, Dim, Frame>& conformal_spatial_metric,
      39             :     const tnsr::I<DataType, Dim, Frame>&
      40             :         gamma_hat_minus_contracted_conformal_christoffel);
      41             : /// @}
      42             : 
      43             : /// @{
      44             : /*!
      45             :  * \brief Computes the spatial part of the upper Z4 constraint
      46             :  *
      47             :  * \details Computes the constraint as:
      48             :  *
      49             :  * \f{align}
      50             :  *     Z^i &= \frac{1}{2} \phi^2 \left(
      51             :  *         \hat{\Gamma}^i - \tilde{\Gamma}^i\right)
      52             :  * \f}
      53             :  *
      54             :  * where \f$\phi^2 \f$ is the square of the conformal factor defined by
      55             :  * `Ccz4::Tags::ConformalFactorSquared` and
      56             :  * \f$\left(\hat{\Gamma}^i - \tilde{\Gamma}^i\right)\f$ is the CCZ4 temporary
      57             :  * expression defined by
      58             :  * `Ccz4::Tags::GammaHatMinusContractedConformalChristoffel`.
      59             :  */
      60             : template <typename DataType, size_t Dim, typename Frame>
      61           1 : void upper_spatial_z4_constraint(
      62             :     const gsl::not_null<tnsr::I<DataType, Dim, Frame>*> result,
      63             :     const Scalar<DataType>& half_conformal_factor_squared,
      64             :     const tnsr::I<DataType, Dim, Frame>&
      65             :         gamma_hat_minus_contracted_conformal_christoffel);
      66             : 
      67             : template <typename DataType, size_t Dim, typename Frame>
      68           1 : tnsr::I<DataType, Dim, Frame> upper_spatial_z4_constraint(
      69             :     const Scalar<DataType>& half_conformal_factor_squared,
      70             :     const tnsr::I<DataType, Dim, Frame>&
      71             :         gamma_hat_minus_contracted_conformal_christoffel);
      72             : /// @}
      73             : }  // namespace Ccz4