OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-29 20:43:49 -06:00
Code Issues Packages Projects Releases Wiki Activity

Catch up to IncompTpfa needing a linear solver.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-02-22 09:28:57 +01:00
parent 5516ee7e04
commit 5b01436db9
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
examples/polymer_reorder.cpp
View File

@ -36,7 +36,7 @@
#include <opm/core/fluid/IncompPropertiesBasic.hpp>
#include <opm/core/fluid/IncompPropertiesFromDeck.hpp>
#include <opm/core/transport/CSRMatrixUmfpackSolver.hpp>
#include <opm/core/linalg/LinearSolverUmfpack.hpp>
#include <opm/polymer/polymertransport.hpp>
#include <opm/polymer/polymermodel.hpp>
@ -326,7 +326,8 @@ main(int argc, char** argv)
double tot_porevol = std::accumulate(porevol.begin(), porevol.end(), 0.0);
// Solvers init.
Opm::IncompTpfa psolver(*grid->c_grid(), props->permeability(), 0);
Opm::LinearSolverUmfpack linsolver;
Opm::IncompTpfa psolver(*grid->c_grid(), props->permeability(), 0, linsolver);
const int method = param.getDefault("method", 1);
const double nltol = param.getDefault("nl_tolerance", 1e-9);