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

Store k_batch at beginning of finalize_batch.

Browse source
This commit is contained in:
Paul Romano 2012-12-04 11:02:53 -05:00
parent a5057a424d
commit 219a34b91f
2 changed files with 5 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

7
src/eigenvalue.F90
View file

@ -152,6 +152,11 @@ contains
integer :: i ! loop index for state point batches
! Before accumulating results for global_tallies, we need to get the
! current batch estimate of k_analog for displaying to output
if (run_mode == MODE_EIGENVALUE) k_batch(current_batch) = &
global_tallies(K_ANALOG) % value
! Collect tallies
if (tallies_on) then
call timer_start(time_tallies)
@ -545,8 +550,6 @@ contains
! =========================================================================
! SINGLE-BATCH ESTIMATE OF K-EFFECTIVE
if (.not. active_batches) k_batch(current_batch) = global_tallies(K_ANALOG) % value
#ifdef MPI
if ((.not. active_batches) .or. (.not. reduce_tallies)) then
! Reduce value of k_batch if running in parallel

6
src/tally.F90
View file

@ -1506,12 +1506,6 @@ contains
curr_ptr => curr_ptr % next
end do
if (run_mode == MODE_EIGENVALUE) then
! Before accumulating results for global_tallies, we need to get the
! current batch estimate of k_analog for displaying to output
if (active_batches) k_batch(current_batch) = global_tallies(K_ANALOG) % value
end if
! Accumulate results for global tallies
if (active_batches) call accumulate_result(global_tallies)
end if