crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 05:35:49 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Speed up reaction rate lookup for FluxCollapseHelper (#3724)

Browse source
This commit is contained in:
Jonathan Shimwell 2026-01-13 21:34:07 +01:00 • committed by GitHub
parent 0486e433d2
commit 84a413b130
No known key found for this signature in database
GPG key ID: B5690EEEBB952194
1 changed files with 4 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

8
openmc/deplete/helpers.py
View file

@ -379,17 +379,17 @@ class FluxCollapseHelper(ReactionRateHelper):
nuclides_direct = self._rate_tally.nuclides
shape = (len(nuclides_direct), len(self._reactions_direct))
rx_rates = self.rate_tally_means[mat_index].reshape(shape)
direct_rx_index = {score: i for i, score in enumerate(self._reactions_direct)}
direct_nuc_index = {nuc: i for i, nuc in enumerate(nuclides_direct)}
mat = self._materials[mat_index]
for name, i_nuc in zip(self.nuclides, nuc_index):
for mt, score, i_rx in zip(self._mts, self._scores, react_index):
if score in self._reactions_direct and name in nuclides_direct:
# Determine index in rx_rates
i_rx_direct = self._reactions_direct.index(score)
i_nuc_direct = nuclides_direct.index(name)
# Get reaction rate from tally
i_rx_direct = direct_rx_index[score]
i_nuc_direct = direct_nuc_index[name]
self._results_cache[i_nuc, i_rx] = rx_rates[i_nuc_direct, i_rx_direct]
else:
# Use flux to collapse reaction rate (per N)