From 0a3c4e489ae7ca01b20d3686cec969ecc9f27987 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Thu, 6 Aug 2020 09:03:13 -0500 Subject: [PATCH] If source rate is 0, don't run OpenMC transport solve --- openmc/deplete/operator.py | 11 ++++++++++- 1 file changed, 10 insertions(+), 1 deletion(-) diff --git a/openmc/deplete/operator.py b/openmc/deplete/operator.py index a11f6d95d..fd9f1bf31 100644 --- a/openmc/deplete/operator.py +++ b/openmc/deplete/operator.py @@ -276,6 +276,16 @@ class Operator(TransportOperator): Eigenvalue and reaction rates resulting from transport operator """ + # Reset results in OpenMC + openmc.lib.reset() + + # If the source rate is zero, return zero reaction rates without running + # a transport solve + if power == 0.0: + rates = self.reaction_rates.copy() + rates.fill(0.0) + return OperatorResult(ufloat(0.0, 0.0), rates) + # Prevent OpenMC from complaining about re-creating tallies openmc.reset_auto_ids() @@ -290,7 +300,6 @@ class Operator(TransportOperator): self._yield_helper.update_tally_nuclides(nuclides) # Run OpenMC - openmc.lib.reset() openmc.lib.run() openmc.lib.reset_timers()