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

Fix bug with torus distance when particle is coincident

Browse source
This commit is contained in:
Paul Romano 2022-05-02 22:09:06 -05:00
parent f9027893dd
commit 01f26d6071
1 changed files with 4 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/surface.cpp
View file

@ -1026,9 +1026,11 @@ double torus_distance(double x1, double x2, double x3, double u1, double u2,
double c1p = 2 * four_A2 * (u1 * x1 + u2 * x2);
double c0p = four_A2 * (x1 * x1 + x2 * x2);
// Coefficient for equation: a t^4 + b t^3 + c t^2 + d t + e = 0
// Coefficient for equation: a t^4 + b t^3 + c t^2 + d t + e = 0. If the point
// is coincident, the 'e' coefficient should be zero. Explicitly setting it to
// zero helps avoid numerical issues below with root finding.
double coeff[5];
coeff[0] = c0 * c0 - c0p;
coeff[0] = coincident ? 0.0 : c0 * c0 - c0p;
coeff[1] = 2 * c0 * c1 - c1p;
coeff[2] = c1 * c1 + 2 * c0 * c2 - c2p;
coeff[3] = 2 * c1 * c2;