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 "Utilities/Gsl.hpp"
       8             : 
       9             : /// \cond
      10             : namespace gsl {
      11             :   template <typename T>
      12             :   class not_null;
      13             : }  // namespace gsl
      14             : 
      15             : class DataVector;
      16             : 
      17             : template<size_t Dim>
      18             : class Mesh;
      19             : /// \endcond
      20             : 
      21             : namespace Particles::MonteCarlo {
      22             : 
      23             : /// Proper 4-volume of a cell for a time step time_step.
      24             : /// We assume that determinant_spatial_metric is given
      25             : /// in inertial coordinates, hence the need for
      26             : /// det_jacobian_logical_to_inertial
      27           1 : void cell_proper_four_volume_finite_difference(
      28             :     gsl::not_null<Scalar<DataVector>*> cell_proper_four_volume,
      29             :     const Scalar<DataVector>& lapse,
      30             :     const Scalar<DataVector>& sqrt_determinant_spatial_metric, double time_step,
      31             :     const Mesh<3>& mesh,
      32             :     const Scalar<DataVector>& det_jacobian_logical_to_inertial);
      33             : 
      34             : /// 3-volume of a cell in inertial coordinate. Note that this is
      35             : /// the coordinate volume, not the proper volume. This quantity
      36             : /// is needed as a normalization factor for the terms coupling
      37             : /// Monte-Carlo transport to the fluid evolution (as we evolved
      38             : /// densitized fluid variables in inertial coordinates).
      39           1 : void cell_inertial_coordinate_three_volume_finite_difference(
      40             :   gsl::not_null<Scalar<DataVector>* > cell_inertial_three_volume,
      41             :   const Mesh<3>& mesh,
      42             :   const Scalar<DataVector>& det_jacobian_logical_to_inertial);
      43             : 
      44             : } // namespace Particles::MonteCarlo