OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

MultisegmentWellEval: use {fmt} to format error messages

Browse Source
This commit is contained in:
Arne Morten Kvarving 2022-12-20 12:40:37 +01:00
parent 9e789fc144
commit 3083525e8e
1 changed files with 16 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
opm/simulators/wells/MultisegmentWellEval.cpp
View File

@ -36,6 +36,8 @@
#include <opm/simulators/wells/WellInterfaceIndices.hpp> #include <opm/simulators/wells/WellInterfaceIndices.hpp>
#include <opm/simulators/wells/WellState.hpp> #include <opm/simulators/wells/WellState.hpp>
#include <fmt/format.h>
#include <algorithm> #include <algorithm>
#include <array> #include <array>
#include <cassert> #include <cassert>
@ -275,8 +277,11 @@ assembleICDPressureEq(const int seg,
icd_pressure_drop = segments_.pressureDropValve(seg); icd_pressure_drop = segments_.pressureDropValve(seg);
break; break;
default: { default: {
OPM_DEFLOG_THROW(std::runtime_error, "Segment " + std::to_string(this->segmentSet()[seg].segmentNumber()) OPM_DEFLOG_THROW(std::runtime_error,
+ " for well " + baseif_.name() + " is not of ICD type", deferred_logger); fmt::format("Segment {} for well {} is not of ICD type",
this->segmentSet()[seg].segmentNumber(),
baseif_.name()),
deferred_logger);
} }
} }
pressure_equation = pressure_equation - icd_pressure_drop; pressure_equation = pressure_equation - icd_pressure_drop;
@ -336,8 +341,8 @@ getFiniteWellResiduals(const std::vector<Scalar>& B_avg,
} }
} }
if (std::isnan(residual) || std::isinf(residual)) { if (std::isnan(residual) || std::isinf(residual)) {
deferred_logger.debug("nan or inf value for residal get for well " + baseif_.name() deferred_logger.debug(fmt::format("nan or inf value for residual for well {} segment {} eq_idx {}",
+ " segment " + std::to_string(seg) + " eq_idx " + std::to_string(eq_idx)); baseif_.name(), seg, eq_idx));
return {false, residuals}; return {false, residuals};
} }
@ -351,7 +356,7 @@ getFiniteWellResiduals(const std::vector<Scalar>& B_avg,
{ {
const double control_residual = std::abs(linSys_.residual()[0][numWellEq - 1]); const double control_residual = std::abs(linSys_.residual()[0][numWellEq - 1]);
if (std::isnan(control_residual) || std::isinf(control_residual)) { if (std::isnan(control_residual) || std::isinf(control_residual)) {
deferred_logger.debug("nan or inf value for control residal get for well " + baseif_.name()); deferred_logger.debug(fmt::format("nan or inf value for control residual for well {}", baseif_.name()));
return {false, residuals}; return {false, residuals};
} }
residuals[numWellEq] = control_residual; residuals[numWellEq] = control_residual;
@ -390,7 +395,9 @@ getControlTolerance(const WellState& well_state,
control_tolerance = tolerance_wells; control_tolerance = tolerance_wells;
break; break;
default: default:
OPM_DEFLOG_THROW(std::runtime_error, "Unknown well control control types for well " << baseif_.name(), deferred_logger); OPM_DEFLOG_THROW(std::runtime_error,
fmt::format("Unknown well control control types for well {}", baseif_.name()),
deferred_logger);
} }
} }
@ -416,7 +423,9 @@ getControlTolerance(const WellState& well_state,
control_tolerance = tolerance_wells; // smaller tolerance for rate control control_tolerance = tolerance_wells; // smaller tolerance for rate control
break; break;
default: default:
OPM_DEFLOG_THROW(std::runtime_error, "Unknown well control control types for well " << baseif_.name(), deferred_logger); OPM_DEFLOG_THROW(std::runtime_error,
fmt::format("Unknown well control control types for well {}", baseif_.name()),
deferred_logger);
} }
} }