-
Notifications
You must be signed in to change notification settings - Fork 555
Figure of Merit implementation #3363
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Conversation
Merge remote-tracking branch 'origin/stat_testing' into stat_testing
Merge remote-tracking branch 'upstream/develop' into stat_testing
I had a quick skim through the code, congrats on this work, I am keen to see this feature added. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@Grego01-biot I just did a bit of clean up here and I'd say it's ready to go. Thanks again for the PR!
@eepeterson let me know if you're good with the changes here now. If so, please approve so we can merge. Thanks! |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
yep all good now.
Description
This PR introduces a new statistical capability in OpenMC: the ability to compute the figure of merit for each tally bin. This could provide the user with more information on how quickly the population variance of a tally bin decreases. This tally reliability indicator is already implemented in MCNP and can be of interest for many users. Indeed, in order to optimize the tally efficiency, the user could run several simulations using different variance reduction techniques and choose the one with the largest FOM.
Fixes # (issue)
At the moment, this functionality is simply a flag defined in the settings of the tally that the user can activate in order to get it for every tally bin. Most of the modifications are made in the Python API since the mean, the variance and the wall time are already available in the statepoint file.
Checklist