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

Added check for zero cross section when calculating distance to collision.

Browse source
This commit is contained in:
Paul Romano 2012-01-18 18:26:16 -05:00
parent 97b8196c46
commit 7ba2b4d9ac
1 changed files with 5 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/physics.F90
View file

@ -79,7 +79,11 @@ contains
call distance_to_boundary(p, d_boundary, surface_crossed, lattice_crossed)
! Sample a distance to collision
d_collision = -log(prn()) / material_xs % total
if (material_xs % total == ZERO) then
d_collision = INFINITY
else
d_collision = -log(prn()) / material_xs % total
end if
! Select smaller of the two distances
distance = min(d_boundary, d_collision)