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

Make sure that the convergece report is stored before throw

Browse Source
This commit is contained in:
Tor Harald Sandve 2019-04-25 10:33:39 +02:00
parent cda7d1906e
commit 8a4117441c
2 changed files with 1 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
opm/autodiff/BlackoilModelEbos.hpp
View File

@ -795,14 +795,6 @@ namespace Opm {
auto report = getReservoirConvergence(timer.currentStepLength(), iteration, B_avg, residual_norms); auto report = getReservoirConvergence(timer.currentStepLength(), iteration, B_avg, residual_norms);
report += wellModel().getWellConvergence(B_avg); report += wellModel().getWellConvergence(B_avg);
// Throw if any NaN or too large residual found.
ConvergenceReport::Severity severity = report.severityOfWorstFailure();
if (severity == ConvergenceReport::Severity::NotANumber) {
OPM_THROW(Opm::NumericalIssue, "NaN residual found!");
} else if (severity == ConvergenceReport::Severity::TooLarge) {
OPM_THROW(Opm::NumericalIssue, "Too large residual found!");
}
return report; return report;
} }

1
opm/simulators/timestepping/AdaptiveTimeSteppingEbos.hpp
View File

@ -491,6 +491,7 @@ namespace Opm {
std::ostringstream msg; std::ostringstream msg;
msg << " Excessive chopping detected in report step " msg << " Excessive chopping detected in report step "
<< sr.back().report_step << ", substep " << sr.back().current_step << "\n"; << sr.back().report_step << ", substep " << sr.back().current_step << "\n";
assert(!sr.back().report.empty());
const auto& wfs = sr.back().report.back().wellFailures(); const auto& wfs = sr.back().report.back().wellFailures();
for (const auto& wf : wfs) { for (const auto& wf : wfs) {
msg << " Well that failed: " << wf.wellName() << "\n"; msg << " Well that failed: " << wf.wellName() << "\n";