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/DataBox/Tag.hpp"
       9              : #include "DataStructures/DataVector.hpp"
      10              : #include "DataStructures/Tensor/Tensor.hpp"
      11              : #include "Domain/Tags.hpp"
      12              : #include "Utilities/Gsl.hpp"
      13              : 
      14              : namespace domain {
      15              : 
      16              : /*!
      17              :  * \brief The "flat logical metric" $\sum_k \frac{\partial x^k}{\partial
      18              :  * \xi^\hat{i}} \frac{\partial x^k}{\partial \xi^\hat{j}}$, which is the flat
      19              :  * spatial metric in element-logical coordinates.
      20              :  *
      21              :  * We define the "flat logical metric" to be the matrix of inner products of the
      22              :  * $\hat{i}$-th logical coordinate basis vector $\partial_{\xi^\hat{i}}$ with
      23              :  * the $\hat{j}$-th logical coordinate basis vector $\partial_{\xi^\hat{j}}$,
      24              :  * where the inner product is taken assuming a flat spatial metric. When
      25              :  * expressed in the ("inertial") $x$-coordinate system, each basis vector is a
      26              :  * column of the Jacobian.
      27              :  */
      28              : template <size_t Dim>
      29            1 : void flat_logical_metric(
      30              :     const gsl::not_null<tnsr::ii<DataVector, Dim, Frame::ElementLogical>*>
      31              :         result,
      32              :     const Jacobian<DataVector, Dim, Frame::ElementLogical, Frame::Inertial>&
      33              :         jacobian);
      34              : 
      35              : namespace Tags {
      36              : /// The flat metric in element-logical coordinates
      37              : ///
      38              : /// \see domain::flat_logical_metric
      39              : template <size_t Dim>
      40            1 : struct FlatLogicalMetric : db::SimpleTag {
      41            0 :   using type = tnsr::ii<DataVector, Dim, Frame::ElementLogical>;
      42              : };
      43              : 
      44              : /// Compute the flat metric in element-logical coordinates from the inverse
      45              : /// Jacobian.
      46              : ///
      47              : /// \see domain::flat_logical_metric
      48              : template <size_t Dim>
      49            1 : struct FlatLogicalMetricCompute : FlatLogicalMetric<Dim>, db::ComputeTag {
      50            0 :   using base = FlatLogicalMetric<Dim>;
      51            0 :   using return_type = typename base::type;
      52            0 :   using argument_tags =
      53              :       tmpl::list<domain::Tags::InverseJacobian<Dim, Frame::ElementLogical,
      54              :                                                Frame::Inertial>>;
      55            0 :   static void function(
      56              :       const gsl::not_null<tnsr::ii<DataVector, Dim, Frame::ElementLogical>*>
      57              :           result,
      58              :       const ::InverseJacobian<DataVector, Dim, Frame::ElementLogical,
      59              :                               Frame::Inertial>& inv_jacobian);
      60              : };
      61              : 
      62              : }  // namespace Tags
      63              : }  // namespace domain