update test to use fewer batches to try and avoid indeterminisim on civet

Author: aprilnovak

src/criticality/CriticalitySearchBase.C

@@ -74,6 +74,13 @@ CriticalitySearchBase::searchForCriticality()
            << std::to_string(_minimum) << " - " << std::to_string(_maximum) << " " << units() << " "
            << std::endl;
 
+  VariadicTable<int, Real, Real, Real> vt(
+      {"Iteration", quantity() + " " + units(), "  k (mean)  ", " k (std dev) "});
+  vt.setColumnFormat({VariadicTableColumnFormat::AUTO,
+                      VariadicTableColumnFormat::SCIENTIFIC,
+                      VariadicTableColumnFormat::SCIENTIFIC,
+                      VariadicTableColumnFormat::SCIENTIFIC});
+
   std::function<Real(Real)> func;
   func = [&](Real x)
   {
@@ -92,6 +99,9 @@ CriticalitySearchBase::searchForCriticality()
     _k_values.push_back(k);
     _k_std_dev_values.push_back(k_std_dev);
 
+    vt.addRow(_k_values.size() - 1, x, k, k_std_dev);
+    vt.print(_console);
+
     if (_tolerance < 3 * k_std_dev)
       mooseDoOnce(mooseWarning(
           "The 'tolerance' for the criticality search (" + std::to_string(_tolerance) +
@@ -109,16 +119,6 @@ CriticalitySearchBase::searchForCriticality()
     mooseError("Failed to converge criticality search! This may happen if your tolerance is too "
                "tight given the statistical error in the computation of k.");
 
-  VariadicTable<int, Real, Real, Real> vt(
-      {"Iteration", quantity() + " " + units(), "  k (mean)  ", " k (std dev) "});
-  vt.setColumnFormat({VariadicTableColumnFormat::AUTO,
-                      VariadicTableColumnFormat::SCIENTIFIC,
-                      VariadicTableColumnFormat::SCIENTIFIC,
-                      VariadicTableColumnFormat::SCIENTIFIC});
-  for (int i = 0; i < _inputs.size(); ++i)
-    vt.addRow(i, _inputs[i], _k_values[i], _k_std_dev_values[i]);
-  vt.print(_console);
-
   // fill the converged value into a postprocessor
   _openmc_problem->setPostprocessorValueByName(_pp_name, _inputs.back());
 }

test/tests/criticality/borated_water/make_openmc_model.py

@@ -50,8 +50,8 @@
 model.geometry = openmc.Geometry([box1, box2])
 
 # Finally, define some run settings
-model.settings.batches = 50
-model.settings.inactive = 10
+model.settings.batches = 10
+model.settings.inactive = 3
 model.settings.particles = 1000
 
 lower_left = (0, 0, 0)

test/tests/criticality/borated_water/openmc.i

@@ -34,7 +34,7 @@
     material_id = 2
     minimum = 0
     maximum = 1000
-    tolerance = 5e-2
+    tolerance = 5e-1
   []
 []
 

test/tests/criticality/borated_water/settings.xml

@@ -2,8 +2,8 @@
 <settings>
   <run_mode>eigenvalue</run_mode>
   <particles>1000</particles>
-  <batches>50</batches>
-  <inactive>10</inactive>
+  <batches>10</batches>
+  <inactive>3</inactive>
   <source particle="neutron" strength="1.0" type="independent">
     <space type="box">
       <parameters>0 0 0 10.0 10.0 10.0</parameters>