OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Use control parameters that don't lead to singular matrices.

Browse Source
This commit is contained in:
Bård Skaflestad 2011-10-06 15:07:00 +02:00
parent b0e35ca07a
commit c56215bdcd
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
examples/spu_2p.cpp
View File

@ -265,7 +265,6 @@ main()
Opm::ImplicitTransportDetails::NRReport rpt;
Opm::ImplicitTransportDetails::NRControl ctrl;
ctrl.max_it = 1;
using Opm::ImplicitTransportLinAlgSupport::CSRMatrixUmfpackSolver;
CSRMatrixUmfpackSolver linsolve;
@ -280,7 +279,8 @@ main()
TransportModel model (fluid, *grid, porevol);
TransportSolver tsolver(model);
double dt = 1;
double dt = 0.5;
ctrl.max_it = 10 ;
tsolver.solve(*grid, tsrc, dt, ctrl, state, linsolve, rpt);
std::cerr << state.saturation() << '\n';