Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <vector>
       7             : 
       8             : #include "DataStructures/DataVector.hpp"
       9             : #include "DataStructures/Tensor/Tensor.hpp"
      10             : #include "NumericalAlgorithms/Spectral/Mesh.hpp"
      11             : #include "Utilities/Gsl.hpp"
      12             : 
      13             : namespace Frame {
      14             : struct Fluid;
      15             : }  // namespace Frame
      16             : 
      17             : /// Items related to the evolution of particles
      18             : /// Items related to Monte-Carlo radiation transport
      19             : namespace Particles::MonteCarlo {
      20             : 
      21             : struct Packet;
      22             : 
      23             : namespace detail {
      24             : 
      25             : // Time derivative of the spatial component of the momentum one-form on a null
      26             : // geodesic
      27             : void time_derivative_momentum_geodesic(
      28             :     gsl::not_null<std::array<double, 3>*> dt_momentum, const Packet& packet,
      29             :     const Scalar<DataVector>& lapse,
      30             :     const tnsr::i<DataVector, 3, Frame::Inertial>& d_lapse,
      31             :     const tnsr::iJ<DataVector, 3, Frame::Inertial>& d_shift,
      32             :     const tnsr::iJJ<DataVector, 3, Frame::Inertial>& d_inv_spatial_metric);
      33             : 
      34             : }  // namespace detail
      35             : 
      36             : // Advances a single packet by time time_step along a geodesic
      37           0 : void evolve_single_packet_on_geodesic(
      38             :     gsl::not_null<Packet*> packet, double time_step,
      39             :     const Scalar<DataVector>& lapse,
      40             :     const tnsr::I<DataVector, 3, Frame::Inertial>& shift,
      41             :     const tnsr::i<DataVector, 3, Frame::Inertial>& d_lapse,
      42             :     const tnsr::iJ<DataVector, 3, Frame::Inertial>& d_shift,
      43             :     const tnsr::iJJ<DataVector, 3, Frame::Inertial>& d_inv_spatial_metric,
      44             :     const tnsr::II<DataVector, 3, Frame::Inertial>& inv_spatial_metric,
      45             :     const std::optional<tnsr::I<DataVector, 3, Frame::Inertial>>& mesh_velocity,
      46             :     const InverseJacobian<DataVector, 3, Frame::ElementLogical,
      47             :                           Frame::Inertial>&
      48             :         inverse_jacobian_logical_to_inertial);
      49             : 
      50             : }  // namespace Particles::MonteCarlo