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

stop division by zero in keff trigger if no trigger given

Browse source
This commit is contained in:
Gavin Ridley 2019-06-20 14:03:51 -04:00
parent 8555f914fa
commit d13a84c36b
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;