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

Merge branch 'master' into cmfd_parallel

Browse source
This commit is contained in:
Bryan Herman 2012-02-12 12:40:02 -08:00
commit 80777a7e36
1 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/physics.F90
View file

@ -1163,7 +1163,7 @@ contains
if (frac == ZERO) then
mu = mu0 + (xi - c_k)/p0
else
mu = mu0 + (sqrt(p0*p0 + 2*frac*(xi - c_k))-p0)/frac
mu = mu0 + (sqrt(max(ZERO, p0*p0 + 2*frac*(xi - c_k))) - p0)/frac
end if
else
message = "Unknown interpolation type: " // trim(to_str(interp))
@ -1213,8 +1213,8 @@ contains
! Precompute factors to save flops
sinphi = sin(phi)
cosphi = cos(phi)
a = sqrt(ONE - mu*mu)
b = sqrt(ONE - w0*w0)
a = sqrt(max(ZERO, ONE - mu*mu))
b = sqrt(max(ZERO, ONE - w0*w0))
! Need to treat special case where sqrt(1 - w**2) is close to zero by
! expanding about the v component rather than the w component