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

Added failure checking for solve method.

Browse Source
This commit is contained in:
Kjetil Olsen Lye 2012-05-09 11:28:56 +02:00
parent 5253ce7fbe
commit b5d474beca
1 changed files with 6 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
opm/core/pressure/IncompTpfa.cpp
View File

@ -210,13 +210,16 @@ namespace Opm
F.bc = bcs; F.bc = bcs;
F.totmob = &totmob[0]; F.totmob = &totmob[0];
F.wdp = &wdp[0]; F.wdp = &wdp[0];
int ok = true;
if (rock_comp.empty()) { if (rock_comp.empty()) {
ifs_tpfa_assemble(gg, &F, &trans_[0], &gpress_omegaweighted_[0], h_); ok = ifs_tpfa_assemble(gg, &F, &trans_[0], &gpress_omegaweighted_[0], h_);
} else { } else {
ifs_tpfa_assemble_comprock(gg, &F, &trans_[0], &gpress_omegaweighted_[0], ok = ifs_tpfa_assemble_comprock(gg, &F, &trans_[0], &gpress_omegaweighted_[0],
&porevol[0], &rock_comp[0], dt, &pressure[0], h_); &porevol[0], &rock_comp[0], dt, &pressure[0], h_);
} }
if (!ok) {
THROW("Failed assembling pressure system.");
}
linsolver_.solve(h_->A, h_->b, h_->x); linsolver_.solve(h_->A, h_->b, h_->x);