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

Moved line in ACE law 44 calculation to save a few flops.

Browse source
This commit is contained in:
Paul Romano 2012-07-25 18:09:13 -04:00
parent 36d4a8a998
commit 679e6732ae
1 changed files with 1 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

2
src/physics.F90
View file

@ -1821,8 +1821,8 @@ contains
! Sampled correlated angle from Kalbach-Mann parameters
r3 = prn()
r4 = prn()
T = (TWO*r4 - ONE) * sinh(KM_A)
if (r3 > KM_R) then
T = (TWO*r4 - ONE) * sinh(KM_A)
mu_out = log(T + sqrt(T*T + ONE))/KM_A
else
mu_out = log(r4*exp(KM_A) + (ONE - r4)*exp(-KM_A))/KM_A