crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 14:15:42 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Merge pull request #1267 from gridley/fix_keff_trigger

Browse source 
stop division by zero in keff trigger if no trigger given
This commit is contained in:
Paul Romano 2019-06-25 13:03:47 -05:00 • committed by GitHub
commit 7c82e7c9bc
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
1 changed files with 5 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

8
src/tallies/trigger.cpp
View file

@ -111,7 +111,7 @@ check_tally_triggers(double& ratio, int& tally_id, int& score)
double
check_keff_trigger()
{
if (settings::run_mode != RUN_MODE_EIGENVALUE) return 0.;
if (settings::run_mode != RUN_MODE_EIGENVALUE) return 0.0;
double k_combined[2];
openmc_get_keff(k_combined);
@ -127,8 +127,10 @@ check_keff_trigger()
case TriggerMetric::relative_error:
uncertainty = k_combined[1] / k_combined[0];
break;
case TriggerMetric::not_active:
break;
default:
// If it's an unrecognized TriggerMetric or no keff trigger is on,
// return 0 to stop division by zero where "ratio" is calculated.
return 0.0;
}
double ratio = uncertainty / settings::keff_trigger.threshold;