Skip to content
Snippets Groups Projects
Commit a2850e7f authored by Norbert Grunwald's avatar Norbert Grunwald
Browse files

fatal error to prevent div by zero in relPermBrooksCorey

parent 7c4c6316
No related branches found
No related tags found
No related merge requests found
......@@ -95,6 +95,10 @@ inline std::unique_ptr<RelPermBrooksCorey> createRelPermBrooksCorey(
auto const exponent =
//! \ogs_file_param{prj__media__medium__properties__property__RelPermBrooksCorey__lambda}
config.getConfigParameter<double>("lambda");
if (exponent == 0.)
{
OGS_FATAL("Exponent 'lambda' must be positive.");
}
return std::make_unique<RelPermBrooksCorey>(
residual_liquid_saturation,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment