crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Improved distance to lattice calculation for cases with coincident surfaces.

Browse source
This commit is contained in:
Paul Romano 2012-01-18 20:00:48 -05:00
parent 059784c2db
commit 64a6c6d798
2 changed files with 15 additions and 14 deletions
Download patch file Download diff file Expand all files Collapse all files

3
src/constants.F90
View file

@ -22,7 +22,8 @@ module constants
real(8), parameter :: TINY_BIT = 1e-8_8
! User for precision in geometry
real(8), parameter :: FP_PRECISION = 1e-5_8
real(8), parameter :: FP_PRECISION = 1e-14_8
real(8), parameter :: FP_REL_PRECISION = 1e-5_8
! Maximum number of collisions/crossings
integer, parameter :: MAX_EVENTS = 10000

26
src/geometry.F90
View file

@ -847,7 +847,7 @@ contains
! Check is calculated distance is new minimum
if (d < dist) then
if (abs(d - dist)/dist >= FP_PRECISION) then
if (abs(d - dist)/dist >= FP_REL_PRECISION) then
dist = d
surface_crossed = -cl % surfaces(i)
lattice_crossed = NONE
@ -869,16 +869,16 @@ contains
z = coord % xyz(3)
! determine oncoming edge
x0 = lat % width_x * 0.5_8
y0 = lat % width_y * 0.5_8
x0 = sign(lat % width_x * 0.5_8, u)
y0 = sign(lat % width_y * 0.5_8, v)
! left and right sides
if (u == ZERO) then
if (abs(x - x0) < FP_PRECISION) then
d = INFINITY
elseif (u == ZERO) then
d = INFINITY
elseif (u > 0) then
d = (x0 - x)/u
else
d = -(x + x0)/u
d = (x0 - x)/u
end if
! If the lattice boundary is coincident with the parent cell boundary,
@ -889,7 +889,7 @@ contains
! point precision.
if (d < dist) then
if (abs(d - dist)/dist >= FP_PRECISION) then
if (abs(d - dist)/dist >= FP_REL_PRECISION) then
dist = d
if (u > 0) then
lattice_crossed = LATTICE_RIGHT
@ -901,16 +901,16 @@ contains
end if
! top and bottom sides
if (v == ZERO) then
if (abs(y - y0) < FP_PRECISION) then
d = INFINITY
elseif (v == ZERO) then
d = INFINITY
elseif (v > 0) then
d = (y0 - y)/v
else
d = -(y + y0)/v
d = (y0 - y)/v
end if
if (d < dist) then
if (abs(d - dist)/dist >= FP_PRECISION) then
if (abs(d - dist)/dist >= FP_REL_PRECISION) then
dist = d
if (v > 0) then
lattice_crossed = LATTICE_TOP