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

Revert "check for absolute precision issues when neutron distances to lattice surface and global surface are very small and close"

Browse source 
This reverts commit d51d3f72de.
This commit is contained in:
Bryan Herman 2013-04-18 07:27:31 -07:00 • committed by Paul Romano
parent 7ba985812f
commit 10d3327bc4
1 changed files with 4 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

11
src/geometry.F90
View file

@ -1114,7 +1114,7 @@ contains
! Check is calculated distance is new minimum
if (d < dist) then
if (abs(d - dist)/dist >= FP_REL_PRECISION) then
if (abs(d - dist)/dist >= FP_PRECISION) then
dist = d
surface_crossed = -cl % surfaces(i)
lattice_crossed = NONE
@ -1156,8 +1156,7 @@ contains
! point precision.
if (d < dist) then
if (abs(d - dist)/dist >= FP_REL_PRECISION &
.and. abs(d - dist) >= FP_PRECISION) then
if (abs(d - dist)/dist >= FP_REL_PRECISION) then
dist = d
if (u > 0) then
lattice_crossed = LATTICE_RIGHT
@ -1178,8 +1177,7 @@ contains
end if
if (d < dist) then
if (abs(d - dist)/dist >= FP_REL_PRECISION &
.and. abs(d - dist) >= FP_PRECISION) then
if (abs(d - dist)/dist >= FP_REL_PRECISION) then
dist = d
if (v > 0) then
lattice_crossed = LATTICE_FRONT
@ -1203,8 +1201,7 @@ contains
end if
if (d < dist) then
if (abs(d - dist)/dist >= FP_REL_PRECISION &
.and. abs(d - dist) >= FP_PRECISION) then
if (abs(d - dist)/dist >= FP_REL_PRECISION) then
dist = d
if (w > 0) then
lattice_crossed = LATTICE_TOP