Line data    Source code

       1            1 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : /// \file
       5              : /// Defines tags related to domain quantities
       6              : 
       7              : #pragma once
       8              : 
       9              : #include "Domain/Tags.hpp"
      10              : #include "PointwiseFunctions/GeneralRelativity/Tags.hpp"
      11              : #include "PointwiseFunctions/Hydro/Tags.hpp"
      12              : 
      13              : namespace Frame {
      14              : struct Fluid;
      15              : struct Inertial;
      16              : }  // namespace Frame
      17              : 
      18              : namespace Particles::MonteCarlo {
      19              : 
      20              : /// Inverse Jacobian of the map from inertial coordinate to an orthonormal frame
      21              : /// comoving with the fluid. That frame uses the 4-velocity as its time axis,
      22              : /// and constructs the other members of the tetrads using Gram-Schmidt's
      23              : /// algorithm.
      24            1 : struct InverseJacobianInertialToFluidCompute
      25              :     : domain::Tags::InverseJacobian<4, Frame::Inertial, Frame::Fluid>,
      26              :       db::ComputeTag {
      27            0 :   using base = domain::Tags::InverseJacobian<4, typename Frame::Inertial,
      28              :                                              typename Frame::Fluid>;
      29            0 :   using return_type = typename base::type;
      30            0 :   using argument_tags =
      31              :       tmpl::list<hydro::Tags::SpatialVelocity<DataVector, 3, Frame::Inertial>,
      32              :                  hydro::Tags::LorentzFactor<DataVector>,
      33              :                  gr::Tags::Lapse<DataVector>,
      34              :                  gr::Tags::Shift<DataVector, 3, Frame::Inertial>,
      35              :                  gr::Tags::SpatialMetric<DataVector, 3> >;
      36              : 
      37            0 :   static void function(gsl::not_null<return_type*> inv_jacobian,
      38              :                        const tnsr::I<DataVector, 3>& spatial_velocity,
      39              :                        const Scalar<DataVector>& lorentz_factor,
      40              :                        const Scalar<DataVector>& lapse,
      41              :                        const tnsr::I<DataVector, 3, Frame::Inertial>& shift,
      42              :                        const tnsr::ii<DataVector, 3>& spatial_metric);
      43              : };
      44              : 
      45              : }  // namespace Particles::MonteCarlo