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

Sample isotropic angular distributions as before so that test results match.

Browse source
This commit is contained in:
Paul Romano 2015-11-10 08:16:36 -06:00
parent 2abc9f1e29
commit 414610307d
1 changed files with 17 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

24
src/distribution_header.F90
View file

@ -1,6 +1,6 @@
module distribution_header
use constants, only: ZERO, HALF, ONE, HISTOGRAM, LINEAR_LINEAR
use constants, only: ZERO, HALF, ONE, TWO, PI, HISTOGRAM, LINEAR_LINEAR
use error, only: fatal_error
use math, only: rotate_angle, maxwell_spectrum, watt_spectrum
use random_lcg, only: prn
@ -203,14 +203,24 @@ contains
class(AngleDistribution), intent(in) :: this
real(8) :: uvw(3)
real(8) :: phi
real(8) :: mu
mu = this%mu%sample()
if (mu == ONE) then
uvw(:) = this%reference_uvw
else
uvw(:) = rotate_angle(this%reference_uvw, mu)
end if
select type (polar_cos => this%mu)
type is (Uniform)
phi = TWO*PI*prn()
mu = TWO*prn() - ONE
uvw(1) = mu
uvw(2) = sqrt(ONE - mu*mu) * cos(phi)
uvw(3) = sqrt(ONE - mu*mu) * sin(phi)
class default
mu = polar_cos%sample()
if (mu == ONE) then
uvw(:) = this%reference_uvw
else
uvw(:) = rotate_angle(this%reference_uvw, mu)
end if
end select
end function angle_sample
end module distribution_header