From 5b01436db9603e767c967e08d15f7ce66543b8ce Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= Date: Wed, 22 Feb 2012 09:28:57 +0100 Subject: [PATCH] Catch up to IncompTpfa needing a linear solver. --- examples/polymer_reorder.cpp | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/examples/polymer_reorder.cpp b/examples/polymer_reorder.cpp index caf215a92..3141d946c 100644 --- a/examples/polymer_reorder.cpp +++ b/examples/polymer_reorder.cpp @@ -36,7 +36,7 @@ #include #include -#include +#include #include #include @@ -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);