crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 14:15:42 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Fix distance_to_boundary for cones. Closes #143.

Browse source
This commit is contained in:
Paul Romano 2013-01-31 20:07:33 -05:00
parent 1a8abdf297
commit 8b2636ff32
1 changed files with 6 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

12
src/geometry.F90
View file

@ -960,9 +960,9 @@ contains
y = y - y0
z = z - z0
a = v*v + w*w - r*u*u
k = y*v + z*w + r*x*u
k = y*v + z*w - r*x*u
c = y*y + z*z - r*x*x
quad = k*k - c
quad = k*k - a*c
if (quad < ZERO .or. a == ZERO) then
! no intersection with cone
@ -1007,9 +1007,9 @@ contains
y = y - y0
z = z - z0
a = u*u + w*w - r*v*v
k = x*u + z*w + r*y*v
k = x*u + z*w - r*y*v
c = x*x + z*z - r*y*y
quad = k*k - c
quad = k*k - a*c
if (quad < ZERO .or. a == ZERO) then
! no intersection with cone
@ -1054,9 +1054,9 @@ contains
y = y - y0
z = z - z0
a = u*u + v*v - r*w*w
k = x*u + y*v + r*z*w
k = x*u + y*v - r*z*w
c = x*x + y*y - r*z*z
quad = k*k - c
quad = k*k - a*c
if (quad < ZERO .or. a == ZERO) then
! no intersection with cone