mcdc-project / mcdc-project/mcdc
mesh tally bug with time census
Nobody has claimed this yet.
- Dominant language
- Python
- Stars
- 61
- Forks
- 38
- Avg merge
- 1d 13h
- Merged PRs (30d)
- 10
Description
I was adding a time census version of the inf_shem361 problem in our regression suite for testing some of the cell tally changes I made, but the mesh tally in the problem errors when census is enabled.
The stdev in the final energy group of the mesh tally is a NAN. The mean looks reasonable.
error printed at the end of the final census (it looks like this code happens as tally post-processing):
MCDC\mcdc\kernel.py:2340: RuntimeWarning: invalid value encountered in sqrt
tally[TALLY_SUM_SQ] = np.sqrt(
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start with the time-census version of the inf_shem361 regression problem and inspect tally post-processing around mcdc/kernel.py:2340, where the invalid sqrt warning occurs. Compare the mesh tally's final energy-group mean and standard deviation with census disabled and determine why the variance becomes NaN. Done means the time-census regression runs without the warning and reports a valid mesh-tally standard deviation.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- hpc
- Issue type
- Bug
- Difficulty
- 3/5
- Estimated time
- 1-2 days
- Activity status
- Stale
- Clarity
- Mostly clear
- Newbie friendliness
- 45/100