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 sum of the Ricci scalar and twice the divergence of the
      15             :  * upper spatial Z4 constraint
      16             :  *
      17             :  * \details Computes the expression as:
      18             :  *
      19             :  * \f{align}
      20             :  *     R + 2 \nabla_k Z^k &=
      21             :  *         \phi^2 \tilde{\gamma}^{ij} (R_{ij} + \nabla_i Z_j + \nabla_j Z_i)
      22             :  * \f}
      23             :  *
      24             :  * where \f$R\f$ is the spatial Ricci scalar, \f$Z^i\f$ is the upper spatial Z4
      25             :  * constraint defined by `Ccz4::Tags::Z4ConstraintUp`, \f$phi^2\f$ is the square
      26             :  * of the conformal factor defined by `Ccz4::Tags::ConformalFactorSquared`,
      27             :  * \f$\tilde{\gamma}^{ij}\f$ is the inverse conformal spatial metric defined by
      28             :  * `Ccz4::Tags::InverseConformalMetric`, \f$R_{ij}\f$ is the spatial Ricci
      29             :  * tensor defined by `Ccz4::Tags::Ricci`, and \f$\nabla_j Z_i\f$ is the gradient
      30             :  * of the spatial Z4 constraint defined by `Ccz4::Tags::GradZ4Constraint`.
      31             :  */
      32             : template <typename DataType, size_t Dim, typename Frame>
      33           1 : void ricci_scalar_plus_divergence_z4_constraint(
      34             :     const gsl::not_null<Scalar<DataType>*> result,
      35             :     const Scalar<DataType>& conformal_factor_squared,
      36             :     const tnsr::II<DataType, Dim, Frame>& inverse_conformal_spatial_metric,
      37             :     const tnsr::ii<DataType, Dim, Frame>& spatial_ricci_tensor,
      38             :     const tnsr::ij<DataType, Dim, Frame>& grad_spatial_z4_constraint);
      39             : 
      40             : template <typename DataType, size_t Dim, typename Frame>
      41           1 : Scalar<DataType> ricci_scalar_plus_divergence_z4_constraint(
      42             :     const Scalar<DataType>& conformal_factor_squared,
      43             :     const tnsr::II<DataType, Dim, Frame>& inverse_conformal_spatial_metric,
      44             :     const tnsr::ii<DataType, Dim, Frame>& spatial_ricci_tensor,
      45             :     const tnsr::ij<DataType, Dim, Frame>& grad_spatial_z4_constraint);
      46             : /// @}
      47             : }  // namespace Ccz4