Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include "DataStructures/Tensor/TypeAliases.hpp"
       7              : #include "Evolution/Systems/GeneralizedHarmonic/GaugeSourceFunctions/Gauges.hpp"
       8              : #include "Utilities/Gsl.hpp"
       9              : 
      10              : /// \cond
      11              : class DataVector;
      12              : template <size_t Dim>
      13              : class Mesh;
      14              : /// \endcond
      15              : 
      16              : namespace gh::gauges {
      17              : /*!
      18              :  * \brief Dispatch to the derived gauge condition.
      19              :  *
      20              :  * Which of the arguments to this function are used will depend on the gauge
      21              :  * condition, but since that is a runtime choice we need support for all gauge
      22              :  * conditions.
      23              :  */
      24              : template <class AllSolutionsForChristoffelAnalytic, size_t Dim>
      25            1 : void dispatch(
      26              :     gsl::not_null<tnsr::a<DataVector, Dim, Frame::Inertial>*> gauge_h,
      27              :     gsl::not_null<tnsr::ab<DataVector, Dim, Frame::Inertial>*> d4_gauge_h,
      28              :     const Scalar<DataVector>& lapse,
      29              :     const tnsr::I<DataVector, Dim, Frame::Inertial>& shift,
      30              :     const Scalar<DataVector>& sqrt_det_spatial_metric,
      31              :     const tnsr::II<DataVector, Dim, Frame::Inertial>& inverse_spatial_metric,
      32              :     const tnsr::abb<DataVector, Dim, Frame::Inertial>& d4_spacetime_metric,
      33              :     const Scalar<DataVector>& half_pi_two_normals,
      34              :     const tnsr::i<DataVector, Dim, Frame::Inertial>& half_phi_two_normals,
      35              :     const tnsr::aa<DataVector, Dim, Frame::Inertial>& spacetime_metric,
      36              :     const tnsr::iaa<DataVector, Dim, Frame::Inertial>& phi,
      37              :     const Mesh<Dim>& mesh, double time,
      38              :     const tnsr::I<DataVector, Dim, Frame::Inertial>& inertial_coords,
      39              :     const InverseJacobian<DataVector, Dim, Frame::ElementLogical,
      40              :                           Frame::Inertial>& inverse_jacobian,
      41              :     const GaugeCondition& gauge_condition);
      42              : }  // namespace gh::gauges