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             : #include <optional>
       8             : 
       9             : #include "DataStructures/Tensor/Tensor.hpp"
      10             : #include "Evolution/Systems/GeneralizedHarmonic/Constraints.hpp"
      11             : #include "Evolution/Systems/GeneralizedHarmonic/System.hpp"
      12             : #include "Evolution/Systems/GeneralizedHarmonic/Tags.hpp"
      13             : #include "Evolution/Systems/ScalarTensor/Tags.hpp"
      14             : #include "Parallel/AlgorithmExecution.hpp"
      15             : #include "Parallel/GlobalCache.hpp"
      16             : #include "PointwiseFunctions/GeneralRelativity/Christoffel.hpp"
      17             : #include "PointwiseFunctions/GeneralRelativity/DerivativesOfSpacetimeMetric.hpp"
      18             : #include "PointwiseFunctions/GeneralRelativity/DetAndInverseSpatialMetric.hpp"
      19             : #include "PointwiseFunctions/GeneralRelativity/ExtrinsicCurvature.hpp"
      20             : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/ConstraintDampingTags.hpp"
      21             : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/ConstraintGammas.hpp"
      22             : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/DerivSpatialMetric.hpp"
      23             : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/ExtrinsicCurvature.hpp"
      24             : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/SpatialDerivOfLapse.hpp"
      25             : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/SpatialDerivOfShift.hpp"
      26             : #include "PointwiseFunctions/GeneralRelativity/InverseSpacetimeMetric.hpp"
      27             : #include "PointwiseFunctions/GeneralRelativity/Lapse.hpp"
      28             : #include "PointwiseFunctions/GeneralRelativity/Ricci.hpp"
      29             : #include "PointwiseFunctions/GeneralRelativity/Shift.hpp"
      30             : #include "PointwiseFunctions/GeneralRelativity/SpacetimeMetric.hpp"
      31             : #include "PointwiseFunctions/GeneralRelativity/SpacetimeNormalOneForm.hpp"
      32             : #include "PointwiseFunctions/GeneralRelativity/SpacetimeNormalVector.hpp"
      33             : #include "PointwiseFunctions/GeneralRelativity/SpatialMetric.hpp"
      34             : #include "PointwiseFunctions/GeneralRelativity/Tags.hpp"
      35             : #include "PointwiseFunctions/GeneralRelativity/WeylElectric.hpp"
      36             : #include "PointwiseFunctions/GeneralRelativity/WeylMagnetic.hpp"
      37             : #include "PointwiseFunctions/ScalarTensor/ConstraintDampingTags.hpp"
      38             : #include "PointwiseFunctions/ScalarTensor/ConstraintGammas.hpp"
      39             : #include "PointwiseFunctions/ScalarTensor/ScalarGaussBonnet/ScalarSource.hpp"
      40             : #include "PointwiseFunctions/ScalarTensor/ScalarGaussBonnet/Tags.hpp"
      41             : #include "PointwiseFunctions/ScalarTensor/SourceTags.hpp"
      42             : #include "Utilities/TMPL.hpp"
      43             : 
      44             : /// \cond
      45             : namespace Parallel {
      46             : template <typename Metavariables>
      47             : class GlobalCache;
      48             : }  // namespace Parallel
      49             : /// \endcond
      50             : 
      51             : namespace ScalarTensor {
      52           0 : namespace Initialization {
      53             : /// \brief List of basic compute tags to initialize the the ScalarTensor without
      54             : /// scalar sources
      55             : template <size_t Dim, typename Fr = Frame::Inertial>
      56           1 : using scalar_tensor_basic_compute_tags = tmpl::list<
      57             :     // Needed to compute the characteristic speeds for the AH finder
      58             :     gr::Tags::SpatialMetricCompute<DataVector, Dim, Fr>,
      59             :     gr::Tags::DetAndInverseSpatialMetricCompute<DataVector, Dim, Fr>,
      60             :     gr::Tags::ShiftCompute<DataVector, Dim, Fr>,
      61             :     gr::Tags::LapseCompute<DataVector, Dim, Fr>,
      62             : 
      63             :     gr::Tags::SpacetimeNormalVectorCompute<DataVector, Dim, Fr>,
      64             :     gh::Tags::DerivLapseCompute<Dim, Fr>,
      65             : 
      66             :     gr::Tags::InverseSpacetimeMetricCompute<DataVector, Dim, Fr>,
      67             :     gh::Tags::DerivShiftCompute<Dim, Fr>,
      68             : 
      69             :     gh::Tags::DerivSpatialMetricCompute<Dim, Fr>,
      70             : 
      71             :     // Compute tags for Trace of Christoffel and Extrinsic curvature
      72             :     gr::Tags::SpatialChristoffelFirstKindCompute<DataVector, Dim, Fr>,
      73             :     gr::Tags::SpatialChristoffelSecondKindCompute<DataVector, Dim, Fr>,
      74             :     gr::Tags::TraceSpatialChristoffelSecondKindCompute<DataVector, Dim, Fr>,
      75             :     gh::Tags::ExtrinsicCurvatureCompute<Dim, Fr>,
      76             :     gh::Tags::TraceExtrinsicCurvatureCompute<Dim, Fr>,
      77             : 
      78             :     // Compute constraint damping parameters.
      79             :     gh::Tags::ConstraintGamma0Compute<Dim, Frame::Grid>,
      80             :     gh::Tags::ConstraintGamma1Compute<Dim, Frame::Grid>,
      81             :     gh::Tags::ConstraintGamma2Compute<Dim, Frame::Grid>,
      82             : 
      83             :     ScalarTensor::Tags::ConstraintGamma1Compute<Dim, Frame::Grid>,
      84             :     ScalarTensor::Tags::ConstraintGamma2Compute<Dim, Frame::Grid>,
      85             : 
      86             :     ScalarTensor::Tags::ScalarSourceCompute>;
      87             : 
      88             : /// \brief List of compute tags to the coupling to curvature
      89             : template <size_t Dim, typename Fr = Frame::Inertial>
      90           1 : using sgb_extra_compute_tags = tmpl::list<
      91             :     ::Tags::DerivTensorCompute<
      92             :         gr::Tags::ExtrinsicCurvature<DataVector, Dim, Fr>,
      93             :         ::domain::Tags::InverseJacobian<Dim, ::Frame::ElementLogical,
      94             :                                         ::Frame::Inertial>,
      95             :         ::domain::Tags::Mesh<Dim>>,
      96             :     gr::Tags::CovariantDerivativeOfExtrinsicCurvatureCompute<Dim, Fr>,
      97             :     ::Tags::DerivTensorCompute<
      98             :         gr::Tags::SpatialChristoffelSecondKind<DataVector, Dim, Fr>,
      99             :         ::domain::Tags::InverseJacobian<Dim, Frame::ElementLogical,
     100             :                                         Frame::Inertial>,
     101             :         ::domain::Tags::Mesh<Dim>>,
     102             :     gr::Tags::SpatialRicciCompute<DataVector, Dim, Fr>,
     103             :     gr::Tags::SpatialRicciScalarCompute<DataVector, Dim, Fr>,
     104             :     gr::Tags::WeylElectricCompute<DataVector, Dim, Fr>,
     105             :     gr::Tags::WeylElectricScalarCompute<DataVector, Dim, Fr>,
     106             :     gr::Tags::SqrtDetSpatialMetricCompute<DataVector, Dim, Fr>,
     107             :     gr::Tags::WeylMagneticCompute<DataVector, Dim, Fr>,
     108             :     gr::Tags::WeylMagneticScalarCompute<DataVector, Dim, Fr>>;
     109             : 
     110             : /// \brief List of compute tags to be initialized in the ScalarTensor system
     111             : ///
     112             : /// \details The compute tags required include those specified in
     113             : /// ::gh::Actions::InitializeGhAnd3Plus1Variables as well as the tags required
     114             : /// to compute spacetime quantities appearing in the scalar evolution equations.
     115             : /// Namely, we include the compute tags associated to the trace of the extrinsic
     116             : /// curvature and the trace of the spatial Christoffel symbol, as well as the
     117             : /// compute tag required to calculate the source term of the scalar equation.
     118             : template <size_t Dim, typename Fr = Frame::Inertial>
     119           1 : using scalar_tensor_3plus1_compute_tags =
     120             :     tmpl::append<scalar_tensor_basic_compute_tags<Dim, Fr>,
     121             :                  sgb_extra_compute_tags<Dim, Fr>>;
     122             : }  // namespace Initialization
     123             : 
     124             : namespace Actions {
     125           0 : struct InitializeGhAnd3Plus1Variables {
     126           0 :   static constexpr size_t volume_dim = 3;
     127           0 :   using frame = Frame::Inertial;
     128           0 :   using compute_tags = db::AddComputeTags<
     129             :       Initialization::scalar_tensor_3plus1_compute_tags<volume_dim, frame>>;
     130             : 
     131           0 :   using const_global_cache_tags = tmpl::list<
     132             :       gh::Tags::DampingFunctionGamma0<volume_dim, Frame::Grid>,
     133             :       gh::Tags::DampingFunctionGamma1<volume_dim, Frame::Grid>,
     134             :       gh::Tags::DampingFunctionGamma2<volume_dim, Frame::Grid>,
     135             :       ScalarTensor::Tags::DampingFunctionGamma1<volume_dim, Frame::Grid>,
     136             :       ScalarTensor::Tags::DampingFunctionGamma2<volume_dim, Frame::Grid>,
     137             :       ScalarTensor::Tags::RampUpParameters,
     138             :       ScalarTensor::Tags::CouplingParameters, ScalarTensor::Tags::ScalarMass>;
     139             : 
     140             :   template <typename DbTagsList, typename... InboxTags, typename Metavariables,
     141             :             typename ArrayIndex, typename ActionList,
     142             :             typename ParallelComponent>
     143           0 :   static Parallel::iterable_action_return_t apply(
     144             :       db::DataBox<DbTagsList>& /*box*/,
     145             :       const tuples::TaggedTuple<InboxTags...>& /*inboxes*/,
     146             :       const Parallel::GlobalCache<Metavariables>& /*cache*/,
     147             :       const ArrayIndex& /*array_index*/, const ActionList /*meta*/,
     148             :       const ParallelComponent* const /*meta*/) {
     149             :     return {Parallel::AlgorithmExecution::Continue, std::nullopt};
     150             :   }
     151             : };
     152             : }  // namespace Actions
     153             : 
     154             : }  // namespace ScalarTensor