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

Use N-3 degrees of freedom for combined k estimator.

Browse source
This commit is contained in:
Paul Romano 2013-06-18 21:35:38 -04:00
parent 8b6f059c0b
commit 19cf65fe3b
1 changed files with 4 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

5
src/output.F90
View file

@ -1461,7 +1461,10 @@ contains
global_tallies(:) % sum_sq = t_value * global_tallies(:) % sum_sq
! Adjust combined estimator
k_combined(2) = t_value * k_combined(2)
if (n_realizations > 3) then
t_value = t_percentile(ONE - alpha/TWO, n_realizations - 3)
k_combined(2) = t_value * k_combined(2)
end if
end if
! write global tallies