Random Ray Missed Cell Policy Change for Adjoint Mode (#3434)

src/random_ray/flat_source_domain.cpp

@@ -298,7 +298,7 @@ int64_t FlatSourceDomain::add_source_to_scalar_flux()
         // in the cell we will use the previous iteration's flux estimate. This
         // injects a small degree of correlation into the simulation, but this
         // is going to be trivial when the miss rate is a few percent or less.
-        if (source_regions_.external_source_present(sr) && !adjoint_) {
+        if (source_regions_.external_source_present(sr)) {
           set_flux_to_old_flux(sr, g);
         } else {
           set_flux_to_source(sr, g);

tests/regression_tests/weightwindows_fw_cadis_mesh/linear/results_true.dat

@@ -60,7 +60,7 @@ Lower Bounds
 1.68e-01
 1.70e-01
 1.75e-01
-1.73e-01
+1.72e-01
 1.73e-01
 1.65e-01
 1.84e-01
@@ -281,7 +281,7 @@ Lower Bounds
 1.86e-01
 1.81e-01
 1.83e-01
-1.73e-01
+1.74e-01
 1.67e-01
 1.72e-01
 1.75e-01
@@ -722,7 +722,7 @@ Lower Bounds
 1.63e-01
 1.58e-01
 1.70e-01
-1.70e-01
+1.69e-01
 1.84e-01
 2.44e-01
 2.71e-01
@@ -775,7 +775,7 @@ Lower Bounds
 1.69e-01
 1.67e-01
 1.76e-01
-1.59e-01
+1.64e-01
 1.59e-01
 1.58e-01
 1.64e-01
@@ -835,7 +835,7 @@ Lower Bounds
 1.76e-01
 1.78e-01
 1.68e-01
-1.78e-01
+1.79e-01
 1.73e-01
 1.69e-01
 1.62e-01
@@ -850,11 +850,11 @@ Lower Bounds
 1.88e-01
 1.77e-01
 1.78e-01
-8.59e-02
+1.76e-01
 1.82e-01
 1.90e-01
 1.76e-01
-1.63e-01
+1.62e-01
 1.69e-01
 1.69e-01
 1.69e-01
@@ -1201,7 +1201,7 @@ Lower Bounds
 1.57e-01
 1.58e-01
 1.58e-01
-1.57e-01
+1.56e-01
 1.58e-01
 1.66e-01
 2.50e-01
@@ -1223,7 +1223,7 @@ Lower Bounds
 2.18e-01
 9.03e-02
 1.68e-01
-1.70e-01
+1.69e-01
 1.64e-01
 1.56e-01
 1.58e-01
@@ -1382,7 +1382,7 @@ Lower Bounds
 1.52e-01
 1.59e-01
 1.52e-01
-1.49e-01
+1.48e-01
 1.55e-01
 2.25e-01
 2.01e-01
@@ -1478,12 +1478,12 @@ Lower Bounds
 1.63e-01
 4.89e-02
 1.64e-01
-6.61e-02
+1.55e-01
 1.57e-01
 1.53e-01
 1.48e-01
 1.53e-01
-7.58e-02
+1.52e-01
 1.40e-01
 1.38e-01
 1.40e-01
@@ -1539,7 +1539,7 @@ Lower Bounds
 2.68e-02
 2.40e-01
 2.25e-01
-2.21e-01
+2.22e-01
 2.27e-01
 2.31e-01
 2.20e-01
@@ -1572,7 +1572,7 @@ Lower Bounds
 8.97e-02
 4.80e-02
 4.87e-02
-6.44e-02
+6.43e-02
 6.38e-02
 4.70e-02
 4.40e-02
@@ -1581,7 +1581,7 @@ Lower Bounds
 4.03e-02
 1.52e-02
 8.66e-03
-6.15e-03
+6.14e-03
 1.56e-01
 1.57e-01
 1.54e-01
@@ -1637,7 +1637,7 @@ Lower Bounds
 1.53e-01
 1.47e-01
 1.42e-01
-1.57e-01
+1.56e-01
 1.66e-01
 2.10e-01
 1.46e-01
@@ -1792,7 +1792,7 @@ Lower Bounds
 4.29e-02
 1.08e-02
 4.27e-03
-7.46e-02
+7.45e-02
 7.65e-02
 5.07e-02
 4.76e-02
@@ -1898,7 +1898,7 @@ Lower Bounds
 1.50e-01
 5.21e-02
 1.53e-01
-1.49e-01
+1.53e-01
 1.52e-01
 1.45e-01
 1.45e-01
@@ -1928,13 +1928,13 @@ Lower Bounds
 7.19e-02
 2.90e-02
 1.48e-01
-1.19e-01
+1.44e-01
 1.48e-01
 1.48e-01
 1.49e-01
 1.54e-01
 1.45e-01
-1.21e-01
+1.38e-01
 1.29e-01
 1.30e-01
 1.31e-01
@@ -1970,7 +1970,7 @@ Lower Bounds
 1.18e-01
 1.31e-01
 1.83e-01
-9.56e-02
+9.55e-02
 2.75e-02
 1.63e-01
 1.47e-01
@@ -2058,7 +2058,7 @@ Lower Bounds
 1.58e-01
 1.43e-01
 1.49e-01
-1.60e-01
+1.59e-01
 2.21e-01
 1.56e-01
 4.57e-02
@@ -2208,7 +2208,7 @@ Lower Bounds
 1.36e-01
 1.33e-01
 1.26e-01
-1.10e-01
+1.11e-01
 1.31e-01
 5.72e-02
 1.80e-02
@@ -2245,7 +2245,7 @@ Lower Bounds
 7.21e-02
 6.08e-02
 6.40e-02
-5.54e-02
+5.53e-02
 6.28e-02
 5.59e-02
 4.39e-02
@@ -2405,7 +2405,7 @@ Lower Bounds
 1.23e-01
 1.24e-01
 1.55e-01
-8.14e-02
+8.15e-02
 2.87e-02
 1.44e-01
 1.47e-01
@@ -2416,7 +2416,7 @@ Lower Bounds
 1.42e-01
 1.37e-01
 1.27e-01
-1.28e-01
+1.29e-01
 1.17e-01
 1.21e-01
 1.63e-01
@@ -2479,12 +2479,12 @@ Lower Bounds
 2.61e-02
 2.32e-02
 1.25e-02
-9.63e-03
+9.62e-03
 6.12e-03
 3.33e-03
 1.57e-01
 1.68e-01
-1.77e-01
+1.76e-01
 1.71e-01
 1.70e-01
 1.76e-01
@@ -2533,7 +2533,7 @@ Lower Bounds
 1.65e-01
 1.69e-01
 1.74e-01
-6.32e-02
+1.78e-01
 1.61e-01
 1.64e-01
 1.57e-01
@@ -2547,8 +2547,8 @@ Lower Bounds
 1.61e-01
 1.66e-01
 1.68e-01
-1.70e-01
-1.63e-01
+1.71e-01
+1.64e-01
 1.61e-01
 1.62e-01
 1.60e-01
@@ -2647,8 +2647,8 @@ Lower Bounds
 1.45e-01
 4.24e-02
 1.27e-02
-1.42e-01
-6.20e-02
+1.43e-01
+1.40e-01
 1.44e-01
 1.45e-01
 1.42e-01
@@ -2660,7 +2660,7 @@ Lower Bounds
 1.15e-01
 1.07e-01
 1.24e-01
-3.83e-02
+3.91e-02
 1.19e-02
 1.85e-01
 1.70e-01
@@ -2678,7 +2678,7 @@ Lower Bounds
 2.05e-02
 7.17e-03
 1.27e-01
-6.47e-02
+6.46e-02
 8.27e-02
 1.04e-01
 7.86e-02
@@ -3043,7 +3043,7 @@ Lower Bounds
 1.57e-01
 1.80e-01
 1.80e-01
-1.24e-01
+1.23e-01
 1.22e-01
 8.18e-02
 7.30e-02
@@ -3055,7 +3055,7 @@ Lower Bounds
 1.13e-01
 1.37e-01
 1.61e-01
-2.21e-01
+2.20e-01
 1.71e-01
 1.51e-01
 1.16e-01
@@ -3187,14 +3187,14 @@ Lower Bounds
 2.64e-02
 1.76e-02
 8.34e-03
-8.72e-02
+8.71e-02
 6.82e-02
 1.03e-01
 9.93e-02
 4.68e-02
 8.00e-02
 6.54e-02
-5.74e-02
+5.73e-02
 5.88e-02
 4.17e-02
 6.05e-02
@@ -3252,7 +3252,7 @@ Lower Bounds
 7.34e-02
 7.91e-02
 4.60e-02
-5.73e-02
+5.72e-02
 6.01e-02
 5.26e-02
 5.01e-02
@@ -3399,7 +3399,7 @@ Upper Bounds
 2.14e+00
 7.89e-01
 8.70e-01
-9.12e-01
+9.11e-01
 9.06e-01
 8.15e-01
 8.92e-01
@@ -3436,7 +3436,7 @@ Upper Bounds
 8.38e-01
 8.52e-01
 8.74e-01
-8.63e-01
+8.62e-01
 8.63e-01
 8.27e-01
 9.20e-01
@@ -3509,8 +3509,8 @@ Upper Bounds
 9.09e-01
 8.56e-01
 8.36e-01
-8.27e-01
-8.09e-01
+8.26e-01
+8.10e-01
 7.69e-01
 7.73e-01
 7.71e-01
@@ -3642,14 +3642,14 @@ Upper Bounds
 8.80e-01
 8.68e-01
 8.83e-01
-8.81e-01
+8.82e-01
 8.83e-01
 8.68e-01
 8.89e-01
 8.61e-01
 9.00e-01
 9.04e-01
-9.07e-01
+9.06e-01
 1.33e+00
 1.66e+00
 5.75e-01
@@ -3657,7 +3657,7 @@ Upper Bounds
 9.28e-01
 9.04e-01
 9.15e-01
-8.67e-01
+8.68e-01
 8.37e-01
 8.58e-01
 8.74e-01
@@ -3673,7 +3673,7 @@ Upper Bounds
 8.97e-01
 8.66e-01
 8.29e-01
-8.49e-01
+8.48e-01
 8.85e-01
 8.47e-01
 8.39e-01
@@ -3734,7 +3734,7 @@ Upper Bounds
 8.72e-01
 8.86e-01
 8.25e-01
-8.06e-01
+8.05e-01
 7.70e-01
 7.66e-01
 7.77e-01
@@ -4084,7 +4084,7 @@ Upper Bounds
 8.23e-01
 8.70e-01
 8.45e-01
-8.81e-01
+8.80e-01
 1.24e+00
 1.66e+00
 6.51e-01
@@ -4098,7 +4098,7 @@ Upper Bounds
 8.14e-01
 7.91e-01
 8.49e-01
-8.48e-01
+8.47e-01
 9.21e-01
 1.22e+00
 1.35e+00
@@ -4124,7 +4124,7 @@ Upper Bounds
 8.56e-01
 8.12e-01
 7.90e-01
-8.10e-01
+8.09e-01
 7.88e-01
 8.08e-01
 8.07e-01
@@ -4141,7 +4141,7 @@ Upper Bounds
 7.98e-01
 8.23e-01
 8.40e-01
-8.12e-01
+8.11e-01
 8.01e-01
 8.15e-01
 8.78e-01
@@ -4149,14 +4149,14 @@ Upper Bounds
 1.49e+00
 5.48e-01
 8.43e-01
-8.37e-01
+8.36e-01
 8.81e-01
-7.97e-01
+8.18e-01
 7.95e-01
 7.90e-01
 8.20e-01
 7.94e-01
-7.75e-01
+7.74e-01
 7.96e-01
 7.56e-01
 8.00e-01
@@ -4176,10 +4176,10 @@ Upper Bounds
 7.46e-01
 7.89e-01
 1.17e+00
-9.25e-01
+9.24e-01
 3.31e-01
 8.47e-01
-8.26e-01
+8.25e-01
 8.14e-01
 8.34e-01
 8.09e-01
@@ -4211,7 +4211,7 @@ Upper Bounds
 8.80e-01
 8.88e-01
 8.39e-01
-8.92e-01
+8.93e-01
 8.67e-01
 8.43e-01
 8.12e-01
@@ -4226,11 +4226,11 @@ Upper Bounds
 9.42e-01
 8.84e-01
 8.89e-01
-4.30e-01
+8.81e-01
 9.12e-01
 9.52e-01
 8.80e-01
-8.13e-01
+8.12e-01
 8.44e-01
 8.47e-01
 8.47e-01
@@ -4284,7 +4284,7 @@ Upper Bounds
 1.28e-01
 6.06e-02
 8.53e-01
-8.42e-01
+8.41e-01
 8.16e-01
 8.39e-01
 8.62e-01
@@ -4293,7 +4293,7 @@ Upper Bounds
 7.57e-01
 7.84e-01
 8.17e-01
-7.86e-01
+7.85e-01
 8.91e-01
 1.29e+00
 1.58e+00
@@ -4384,7 +4384,7 @@ Upper Bounds
 7.51e-01
 7.69e-01
 7.83e-01
-7.95e-01
+7.94e-01
 1.09e+00
 1.03e+00
 3.79e-01
@@ -4439,7 +4439,7 @@ Upper Bounds
 8.29e-01
 8.34e-01
 8.10e-01
-8.16e-01
+8.15e-01
 8.13e-01
 7.91e-01
 7.58e-01
@@ -4500,7 +4500,7 @@ Upper Bounds
 5.38e-01
 4.35e-01
 5.87e-01
-2.74e-01
+2.73e-01
 3.52e-01
 2.53e-01
 1.54e-01
@@ -4577,7 +4577,7 @@ Upper Bounds
 7.84e-01
 7.91e-01
 7.91e-01
-7.83e-01
+7.82e-01
 7.90e-01
 8.30e-01
 1.25e+00
@@ -4590,7 +4590,7 @@ Upper Bounds
 7.98e-01
 7.82e-01
 7.55e-01
-7.52e-01
+7.51e-01
 8.14e-01
 8.19e-01
 7.83e-01
@@ -4599,7 +4599,7 @@ Upper Bounds
 1.09e+00
 4.51e-01
 8.40e-01
-8.48e-01
+8.47e-01
 8.18e-01
 7.81e-01
 7.89e-01
@@ -4623,7 +4623,7 @@ Upper Bounds
 7.25e-01
 7.18e-01
 7.55e-01
-7.87e-01
+7.86e-01
 8.22e-01
 1.16e+00
 5.88e-01
@@ -4658,7 +4658,7 @@ Upper Bounds
 1.05e+00
 6.00e-01
 1.91e-01
-8.27e-01
+8.26e-01
 7.94e-01
 8.07e-01
 8.05e-01
@@ -4758,7 +4758,7 @@ Upper Bounds
 7.61e-01
 7.94e-01
 7.58e-01
-7.43e-01
+7.42e-01
 7.76e-01
 1.12e+00
 1.01e+00
@@ -4784,12 +4784,12 @@ Upper Bounds
 7.60e-01
 7.73e-01
 8.04e-01
-7.87e-01
+7.86e-01
 7.68e-01
 7.69e-01
 7.72e-01
 7.79e-01
-8.20e-01
+8.19e-01
 1.12e+00
 7.80e-01
 3.07e-01
@@ -4817,7 +4817,7 @@ Upper Bounds
 7.52e-01
 7.32e-01
 7.17e-01
-7.61e-01
+7.60e-01
 7.74e-01
 8.04e-01
 1.13e+00
@@ -4854,12 +4854,12 @@ Upper Bounds
 8.15e-01
 2.44e-01
 8.21e-01
-3.30e-01
+7.75e-01
 7.87e-01
 7.66e-01
 7.39e-01
 7.65e-01
-3.79e-01
+7.62e-01
 6.99e-01
 6.88e-01
 6.98e-01
@@ -4962,10 +4962,10 @@ Upper Bounds
 7.87e-01
 7.72e-01
 8.00e-01
-8.25e-01
-8.31e-01
+8.24e-01
+8.30e-01
 8.27e-01
-7.89e-01
+7.88e-01
 7.35e-01
 7.02e-01
 6.79e-01
@@ -4991,13 +4991,13 @@ Upper Bounds
 7.91e-01
 7.84e-01
 8.12e-01
-8.14e-01
+8.13e-01
 7.93e-01
 7.86e-01
 7.80e-01
 7.88e-01
 7.62e-01
-7.26e-01
+7.25e-01
 7.48e-01
 8.00e-01
 1.08e+00
@@ -5013,7 +5013,7 @@ Upper Bounds
 7.67e-01
 7.34e-01
 7.08e-01
-7.83e-01
+7.82e-01
 8.28e-01
 1.05e+00
 7.28e-01
@@ -5043,7 +5043,7 @@ Upper Bounds
 7.26e-01
 7.36e-01
 7.31e-01
-7.41e-01
+7.40e-01
 7.92e-01
 1.07e+00
 6.18e-01
@@ -5089,7 +5089,7 @@ Upper Bounds
 6.75e-01
 6.66e-01
 7.38e-01
-7.32e-01
+7.31e-01
 9.88e-01
 5.41e-01
 1.68e-01
@@ -5219,7 +5219,7 @@ Upper Bounds
 7.96e-01
 8.02e-01
 7.57e-01
-7.80e-01
+7.79e-01
 7.99e-01
 7.43e-01
 7.26e-01
@@ -5251,7 +5251,7 @@ Upper Bounds
 7.54e-01
 7.53e-01
 7.66e-01
-7.41e-01
+7.40e-01
 7.01e-01
 7.50e-01
 8.07e-01
@@ -5274,9 +5274,9 @@ Upper Bounds
 7.50e-01
 2.61e-01
 7.64e-01
-7.45e-01
+7.65e-01
 7.60e-01
-7.26e-01
+7.25e-01
 7.24e-01
 7.17e-01
 6.99e-01
@@ -5294,7 +5294,7 @@ Upper Bounds
 7.14e-01
 7.38e-01
 7.21e-01
-6.81e-01
+6.80e-01
 6.69e-01
 6.84e-01
 6.84e-01
@@ -5304,13 +5304,13 @@ Upper Bounds
 3.59e-01
 1.45e-01
 7.38e-01
-5.93e-01
+7.18e-01
 7.38e-01
 7.39e-01
 7.46e-01
 7.70e-01
 7.25e-01
-6.07e-01
+6.90e-01
 6.46e-01
 6.49e-01
 6.56e-01
@@ -5318,7 +5318,7 @@ Upper Bounds
 9.65e-01
 3.98e-01
 1.02e-01
-7.63e-01
+7.64e-01
 7.31e-01
 7.32e-01
 7.05e-01
@@ -5369,7 +5369,7 @@ Upper Bounds
 1.11e+00
 1.11e+00
 1.08e+00
-9.96e-01
+9.95e-01
 1.02e+00
 9.59e-01
 1.01e+00
@@ -5434,7 +5434,7 @@ Upper Bounds
 7.89e-01
 7.16e-01
 7.47e-01
-7.98e-01
+7.97e-01
 1.11e+00
 7.82e-01
 2.29e-01
@@ -5447,7 +5447,7 @@ Upper Bounds
 7.75e-01
 7.99e-01
 8.00e-01
-7.27e-01
+7.26e-01
 7.00e-01
 7.95e-01
 1.12e+00
@@ -5584,7 +5584,7 @@ Upper Bounds
 6.82e-01
 6.66e-01
 6.30e-01
-5.52e-01
+5.53e-01
 6.56e-01
 2.86e-01
 9.02e-02
@@ -5603,7 +5603,7 @@ Upper Bounds
 5.90e-01
 1.70e-01
 6.04e-02
-9.35e-01
+9.34e-01
 7.67e-01
 7.48e-01
 8.36e-01
@@ -5688,7 +5688,7 @@ Upper Bounds
 7.74e-01
 7.69e-01
 7.48e-01
-7.27e-01
+7.26e-01
 7.70e-01
 1.04e+00
 5.38e-01
@@ -5698,7 +5698,7 @@ Upper Bounds
 7.98e-01
 7.72e-01
 8.10e-01
-7.95e-01
+7.96e-01
 7.80e-01
 7.69e-01
 7.86e-01
@@ -5776,7 +5776,7 @@ Upper Bounds
 7.07e-01
 6.87e-01
 6.74e-01
-6.34e-01
+6.35e-01
 6.34e-01
 6.14e-01
 6.22e-01
@@ -5792,10 +5792,10 @@ Upper Bounds
 7.08e-01
 6.86e-01
 6.37e-01
-6.42e-01
+6.43e-01
 5.84e-01
 6.05e-01
-8.13e-01
+8.14e-01
 2.50e-01
 7.65e-02
 7.35e-01
@@ -5805,7 +5805,7 @@ Upper Bounds
 7.24e-01
 7.52e-01
 7.60e-01
-7.36e-01
+7.35e-01
 6.72e-01
 6.44e-01
 6.12e-01
@@ -5860,7 +5860,7 @@ Upper Bounds
 1.67e-02
 7.86e-01
 8.39e-01
-8.83e-01
+8.82e-01
 8.55e-01
 8.50e-01
 8.79e-01
@@ -5901,7 +5901,7 @@ Upper Bounds
 7.49e-01
 8.31e-01
 1.07e+00
-4.89e-01
+4.88e-01
 1.30e-01
 8.26e-01
 8.49e-01
@@ -5909,7 +5909,7 @@ Upper Bounds
 8.27e-01
 8.44e-01
 8.70e-01
-3.16e-01
+8.88e-01
 8.04e-01
 8.19e-01
 7.85e-01
@@ -5923,8 +5923,8 @@ Upper Bounds
 8.03e-01
 8.28e-01
 8.39e-01
-8.52e-01
-8.17e-01
+8.53e-01
+8.18e-01
 8.07e-01
 8.12e-01
 8.02e-01
@@ -6023,8 +6023,8 @@ Upper Bounds
 7.24e-01
 2.12e-01
 6.37e-02
-7.12e-01
-3.10e-01
+7.13e-01
+6.99e-01
 7.21e-01
 7.24e-01
 7.11e-01
@@ -6036,7 +6036,7 @@ Upper Bounds
 5.74e-01
 5.37e-01
 6.18e-01
-1.92e-01
+1.95e-01
 5.93e-02
 9.23e-01
 8.49e-01
@@ -6053,8 +6053,8 @@ Upper Bounds
 3.58e-01
 1.02e-01
 3.58e-02
-6.34e-01
-3.24e-01
+6.33e-01
+3.23e-01
 4.14e-01
 5.18e-01
 3.93e-01
@@ -6067,7 +6067,7 @@ Upper Bounds
 2.08e-01
 1.27e-01
 3.00e-02
-9.84e-03
+9.83e-03
 2.53e-01
 1.24e-01
 1.72e-01
@@ -6211,7 +6211,7 @@ Upper Bounds
 1.03e+00
 1.15e+00
 1.11e+00
-9.64e-01
+9.63e-01
 1.01e+00
 9.05e-01
 8.67e-01
@@ -6229,8 +6229,8 @@ Upper Bounds
 8.93e-01
 8.70e-01
 8.58e-01
-8.56e-01
-8.00e-01
+8.57e-01
+8.01e-01
 2.31e-01
 7.13e-02
 9.12e-01
@@ -6249,7 +6249,7 @@ Upper Bounds
 2.01e-01
 7.22e-02
 7.80e-01
-8.00e-01
+7.99e-01
 9.64e-01
 9.33e-01
 9.70e-01
@@ -6277,7 +6277,7 @@ Upper Bounds
 5.08e-01
 3.32e-01
 9.24e-02
-4.15e-02
+4.14e-02
 3.42e-01
 2.04e-01
 2.56e-01
@@ -6414,12 +6414,12 @@ Upper Bounds
 1.17e-01
 7.24e-02
 5.31e-01
-7.17e-01
+7.16e-01
 8.03e-01
 7.87e-01
 9.01e-01
 9.00e-01
-6.18e-01
+6.17e-01
 6.08e-01
 4.09e-01
 3.65e-01
@@ -6545,7 +6545,7 @@ Upper Bounds
 2.87e-01
 1.88e-01
 1.85e-01
-1.80e-01
+1.79e-01
 1.07e-01
 3.77e-02
 5.65e-01
@@ -6578,13 +6578,13 @@ Upper Bounds
 1.41e-01
 5.88e-02
 6.26e-02
-4.47e-01
+4.46e-01
 3.47e-01
 5.27e-01
 6.11e-01
 3.91e-01
 4.22e-01
-4.77e-01
+4.76e-01
 4.54e-01
 4.31e-01
 3.19e-01
@@ -6727,14 +6727,14 @@ Upper Bounds
 5.73e-02
 5.81e-02
 1.92e-02
-6.10e-03
+6.09e-03
 2.72e-02
 4.77e-02
 4.74e-02
 3.34e-02
 5.41e-02
 7.93e-02
-6.70e-02
+6.69e-02
 7.67e-02
 9.31e-02
 5.17e-02