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 omp reduction instead of atomic on global tallies

Browse source
This commit is contained in:
Paul Romano 2015-06-03 16:42:55 +07:00
parent 825f037dba
commit 52b31d8ac9
5 changed files with 34 additions and 16 deletions
Download patch file Download diff file Expand all files Collapse all files

26
src/eigenvalue.F90
View file

@ -75,7 +75,11 @@ contains
! ====================================================================
! LOOP OVER PARTICLES
!$omp parallel do schedule(static) firstprivate(p)
!$omp parallel do schedule(static) firstprivate(p) &
!$omp reduction(+:tally_tracklength) &
!$omp reduction(+:tally_collision) &
!$omp reduction(+:tally_leakage) &
!$omp reduction(+:tally_absorption)
PARTICLE_LOOP: do i_work = 1, work
current_work = i_work
@ -171,6 +175,22 @@ contains
subroutine finalize_generation()
! Update global tallies with the omp private accumulation variables
global_tallies(K_TRACKLENGTH) % value = &
global_tallies(K_TRACKLENGTH) % value + tally_tracklength
global_tallies(K_COLLISION) % value = &
global_tallies(K_COLLISION) % value + tally_collision
global_tallies(LEAKAGE) % value = &
global_tallies(LEAKAGE) % value + tally_leakage
global_tallies(K_ABSORPTION) % value = &
global_tallies(K_ABSORPTION) % value + tally_absorption
! reset private tallies
tally_tracklength = 0
tally_collision = 0
tally_leakage = 0
tally_absorption = 0
#ifdef _OPENMP
! Join the fission bank from each thread into one global fission bank
call join_bank_from_threads()
@ -210,6 +230,10 @@ contains
if (.not. active_batches) then
call reset_result(global_tallies)
n_realizations = 0
tally_tracklength = 0
tally_collision = 0
tally_leakage = 0
tally_absorption = 0
end if
! Perform CMFD calculation if on

4
src/geometry.F90
View file

@ -320,9 +320,7 @@ contains
! Score to global leakage tally
if (tallies_on) then
!$omp atomic
global_tallies(LEAKAGE) % value = &
global_tallies(LEAKAGE) % value + p % wgt
tally_leakage = tally_leakage + p % wgt
end if
! Display message

4
src/global.F90
View file

@ -116,6 +116,10 @@ module global
! 3) leakage fraction
type(TallyResult), allocatable, target :: global_tallies(:)
real(8) :: tally_tracklength = 0
real(8) :: tally_collision = 0
real(8) :: tally_leakage = 0
real(8) :: tally_absorption = 0
! Tally map structure
type(TallyMap), allocatable :: tally_maps(:)

8
src/physics.F90
View file

@ -253,18 +253,14 @@ contains
p % last_wgt = p % wgt
! Score implicit absorption estimate of keff
!$omp atomic
global_tallies(K_ABSORPTION) % value = &
global_tallies(K_ABSORPTION) % value + p % absorb_wgt * &
tally_absorption = tally_absorption + p % absorb_wgt * &
micro_xs(i_nuclide) % nu_fission / micro_xs(i_nuclide) % absorption
else
! See if disappearance reaction happens
if (micro_xs(i_nuclide) % absorption > &
prn() * micro_xs(i_nuclide) % total) then
! Score absorption estimate of keff
!$omp atomic
global_tallies(K_ABSORPTION) % value = &
global_tallies(K_ABSORPTION) % value + p % wgt * &
tally_absorption = tally_absorption + p % wgt * &
micro_xs(i_nuclide) % nu_fission / micro_xs(i_nuclide) % absorption
p % alive = .false.

8
src/tracking.F90
View file

@ -110,9 +110,7 @@ contains
call score_tracklength_tally(p, distance)
! Score track-length estimate of k-eff
!$omp atomic
global_tallies(K_TRACKLENGTH) % value = &
global_tallies(K_TRACKLENGTH) % value + p % wgt * distance * &
tally_tracklength = tally_tracklength + p % wgt * distance * &
material_xs % nu_fission
if (d_collision > d_boundary) then
@ -138,9 +136,7 @@ contains
! PARTICLE HAS COLLISION
! Score collision estimate of keff
!$omp atomic
global_tallies(K_COLLISION) % value = &
global_tallies(K_COLLISION) % value + p % wgt * &
tally_collision = tally_collision + p % wgt * &
material_xs % nu_fission / material_xs % total
! score surface current tallies -- this has to be done before the collision