avoid zero division if source rate of previous result is zero (#3169)

openmc-dev · Oct 14, 2024 · dcb2557 · dcb2557
1 parent 8263f05
commit dcb2557
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/openmc/deplete/abc.py b/openmc/deplete/abc.py
@@ -754,8 +754,9 @@ def _get_bos_data_from_restart(self, source_rate, bos_conc):
         rates = res.rates[0]
         k = ufloat(res.k[0, 0], res.k[0, 1])
 
-        # Scale reaction rates by ratio of source rates
-        rates *= source_rate / res.source_rate
+        if res.source_rate != 0.0:
+            # Scale reaction rates by ratio of source rates
+            rates *= source_rate / res.source_rate
         return bos_conc, OperatorResult(k, rates)
 
     def _get_start_data(self):