Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <array>
       7              : 
       8              : #include "DataStructures/Tensor/IndexType.hpp"
       9              : #include "DataStructures/Tensor/TypeAliases.hpp"
      10              : 
      11              : /// \cond
      12              : namespace gsl {
      13              :   template <typename T>
      14              :   class not_null;
      15              : }  // namespace gsl
      16              : 
      17              : class DataVector;
      18              : namespace Particles::MonteCarlo {
      19              :   class Packet;
      20              : } // namespace Particles::MonteCarlo
      21              : /// \endcond
      22              : 
      23              : namespace Particles::MonteCarlo {
      24              : 
      25              : /// \brief Contribution to the neutrino-matter coupling
      26              : /// terms from evolving a packet by time dt, given
      27              : /// absorption and scattering opacities.
      28              : ///
      29              : /// \details We calculate total momentum exchanges, not
      30              : /// densities; thus when coupling to the evolution of the
      31              : /// energy/momentum density, division by the spatial volume is
      32              : /// necessary. We calculate the source terms of Eqs (62-64) of
      33              : /// \cite Foucart2021mcb, with integrals calculated as in
      34              : /// Eqs (6-10) of that manuscript (except that we do not
      35              : /// divide by the coordinate volume V).
      36            1 : void AddCouplingTermsForPropagation(
      37              :     gsl::not_null<Scalar<DataVector>*> coupling_tilde_tau,
      38              :     gsl::not_null<tnsr::i<DataVector, 3, Frame::Inertial>*> coupling_tilde_s,
      39              :     gsl::not_null<Scalar<DataVector>*> coupling_rho_ye, const Packet& packet,
      40              :     size_t extended_idx, double dt, double absorption_opacity,
      41              :     double scattering_opacity, double fluid_frame_energy, double lapse,
      42              :     double lorentz_factor,
      43              :     const std::array<double, 3>& lower_spatial_four_velocity_packet);
      44              : 
      45              : } // namespace Particles::MonteCarlo