Add optically thick sphere test

sxs-collaboration · Jun 5, 2024 · 72c504a · 72c504a
1 parent 363e582
commit 72c504a
Show file tree

Hide file tree

Showing 4 changed files with 484 additions and 3 deletions.
diff --git a/src/Evolution/Particles/MonteCarlo/EvolvePacketsInElement.tpp b/src/Evolution/Particles/MonteCarlo/EvolvePacketsInElement.tpp
@@ -411,7 +411,7 @@ void TemplatedLocalFunctions<EnergyBins, NeutrinoSpecies>::evolve_packets(
     // only the opacities are updated mid-step. Decide whether to do
     // the same here once we handle ghost zone (the main reason not to
     // update is to limit the amount of ghost zone information we
-    // store.
+    // store).
     packet.index_of_closest_grid_point =
         closest_point_index_3d[0] +
         extents[0] * (closest_point_index_3d[1] +

diff --git a/tests/Unit/Evolution/Particles/MonteCarlo/CMakeLists.txt b/tests/Unit/Evolution/Particles/MonteCarlo/CMakeLists.txt
@@ -10,6 +10,7 @@ set(LIBRARY_SOURCES
   Test_InterpolateOpacities.cpp
   Test_InverseJacobianInertialToFluid.cpp
   Test_NeutrinoInteractionTable.cpp
+  Test_OpticallyThickSphere.cpp
   Test_Packet.cpp
   Test_Scattering.cpp
   Test_TakeTimeStep.cpp

diff --git a/tests/Unit/Evolution/Particles/MonteCarlo/Test_NeutrinoInteractionTable.cpp b/tests/Unit/Evolution/Particles/MonteCarlo/Test_NeutrinoInteractionTable.cpp
@@ -1,10 +1,11 @@
 // Distributed under the MIT License.
 // See LICENSE.txt for details.
 
+#include "Framework/TestingFramework.hpp"
+
 #include "DataStructures/DataVector.hpp"
 #include "DataStructures/Tensor/Tensor.hpp"
 #include "Evolution/Particles/MonteCarlo/NeutrinoInteractionTable.hpp"
-#include "Framework/TestingFramework.hpp"
 #include "Informer/InfoFromBuild.hpp"
 
 namespace {
@@ -183,6 +184,5 @@ SPECTRE_TEST_CASE("Unit.Evolution.Particles.MonteCarloInteractionTable",
                      gsl::at(gsl::at(scattering_opacity, ns), ng)[2]));
     }
   }
-
   test_explicit_interaction_table();
 }