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 "Domain/Tags.hpp"
       8              : #include "Evolution/DgSubcell/Tags/Coordinates.hpp"
       9              : #include "Evolution/DgSubcell/Tags/Mesh.hpp"
      10              : #include "Evolution/Particles/MonteCarlo/Tags.hpp"
      11              : #include "PointwiseFunctions/GeneralRelativity/Tags.hpp"
      12              : #include "Utilities/Gsl.hpp"
      13              : 
      14              : /// \cond
      15              : namespace gsl {
      16              : template <typename T>
      17              : class not_null;
      18              : }  // namespace gsl
      19              : 
      20              : class DataVector;
      21              : 
      22              : template <size_t Dim>
      23              : class Mesh;
      24              : /// \endcond
      25              : 
      26              : namespace Particles::MonteCarlo {
      27              : 
      28            0 : void cell_light_crossing_time(
      29              :     gsl::not_null<Scalar<DataVector>*> cell_light_crossing_time,
      30              :     const Mesh<3>& mesh,
      31              :     const tnsr::I<DataVector, 3, Frame::Inertial>& inertial_coordinates,
      32              :     const Scalar<DataVector>& lapse,
      33              :     const tnsr::I<DataVector, 3, Frame::Inertial>& shift,
      34              :     const tnsr::II<DataVector, 3, Frame::Inertial>& inv_spatial_metric);
      35              : 
      36            0 : struct CellLightCrossingTimeCompute : Tags::CellLightCrossingTime<DataVector>,
      37              :                                            db::ComputeTag {
      38            0 :   using base = Tags::CellLightCrossingTime<DataVector>;
      39            0 :   using return_type = typename base::type;
      40            0 :   using argument_tags = tmpl::list<
      41              :     evolution::dg::subcell::Tags::Mesh<3>,
      42              :     evolution::dg::subcell::Tags::Coordinates<3, Frame::Inertial>,
      43              :     gr::Tags::Lapse<DataVector>,
      44              :     gr::Tags::Shift<DataVector, 3, Frame::Inertial>,
      45              :     gr::Tags::InverseSpatialMetric<DataVector, 3, Frame::Inertial>>;
      46              : 
      47            0 :   static void function(
      48              :     gsl::not_null<return_type*> cell_light_crossing_time_,
      49              :     const Mesh<3>& mesh,
      50              :     const tnsr::I<DataVector, 3, Frame::Inertial>& inertial_coordinates,
      51              :     const Scalar<DataVector>& lapse,
      52              :     const tnsr::I<DataVector, 3, Frame::Inertial>& shift,
      53              :     const tnsr::II<DataVector, 3, Frame::Inertial>& inv_spatial_metric){
      54              :     cell_light_crossing_time(cell_light_crossing_time_, mesh,
      55              :         inertial_coordinates, lapse, shift, inv_spatial_metric);
      56              :   }
      57              : };
      58              : 
      59              : 
      60              : }  // namespace Particles::MonteCarlo