Skip to content
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

cover a previously missed quadric edge case #1812

Merged
merged 2 commits into from
Apr 8, 2021

Conversation

gridley
Copy link
Contributor

@gridley gridley commented Apr 7, 2021

Closes #1811

There is a measure zero set of quadric/particle orientations that lead to division by zero in the quadric intersection routine. One likely example is if a user inputs a general plane as a quadric, which will cause failure of this routine every time at the moment and segfault OpenMC, as shown in #1811. the code comment explains a little more. I think there are other edge cases this covers, but haven't put time into figuring what exactly they are.

Copy link
Contributor

@paulromano paulromano left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Makes sense to me. Thanks @gridley! I'll leave it open for a little bit in case others have thoughts.

Copy link
Contributor

@pshriwise pshriwise left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This approach and explanation look solid to me. Thanks for documenting it well @gridley. Just one comment about control flow and I'm happy.

src/surface.cpp Show resolved Hide resolved
@paulromano paulromano merged commit d11c46b into openmc-dev:develop Apr 8, 2021
pshriwise pushed a commit to pshriwise/openmc that referenced this pull request Apr 24, 2021
cover a previously missed quadric edge case
pshriwise pushed a commit to pshriwise/openmc that referenced this pull request Apr 24, 2021
cover a previously missed quadric edge case
pshriwise pushed a commit to pshriwise/openmc that referenced this pull request Apr 24, 2021
cover a previously missed quadric edge case
pshriwise pushed a commit to pshriwise/openmc that referenced this pull request Apr 26, 2021
cover a previously missed quadric edge case
paulromano added a commit that referenced this pull request Apr 26, 2021
cover a previously missed quadric edge case
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

segfault with particle first crossing quadric
3 participants