From d557b91db059f93df5938dfa007b139b1ca6f54e Mon Sep 17 00:00:00 2001
From: Kjetil Olsen Lye <kjetil.olsen.lye@sintef.no>
Date: Mon, 16 Apr 2012 12:18:37 +0200
Subject: [PATCH] Made use of the new LinearSolverFactory in wells_example

---
 examples/wells_example.cpp | 14 ++++----------
 1 file changed, 4 insertions(+), 10 deletions(-)

diff --git a/examples/wells_example.cpp b/examples/wells_example.cpp
index 6c8d95509..8059cb8e6 100644
--- a/examples/wells_example.cpp
+++ b/examples/wells_example.cpp
@@ -8,16 +8,12 @@
 #include <opm/core/GridManager.hpp>
 #include <opm/core/pressure/IncompTpfa.hpp>
 #include <opm/core/fluid/IncompPropertiesFromDeck.hpp>
-#include <opm/core/linalg/LinearSolverUmfpack.hpp>
 #include <opm/core/newwells.h>
 #include <opm/core/grid.h>
 #include <opm/core/utility/miscUtilities.hpp>
 #include <opm/core/TwophaseState.hpp>
 #include <opm/core/pressure/FlowBCManager.hpp>
-
-#ifdef EXPERIMENT_ISTL
-#include <opm/core/linalg/LinearSolverIstl.hpp>
-#endif
+#include <opm/core/linalg/LinearSolverFactory.hpp>
 int main(int argc, char** argv) {
 
     using namespace Opm::parameter;
@@ -39,11 +35,9 @@ int main(int argc, char** argv) {
     double gravity[3] = {0.0, 0.0, parameters.getDefault<double>("gravity", 0.0)};
     IncompPropertiesFromDeck incomp_properties(parser, global_cells);
 
-#ifdef EXPERIMENT_ISTL
-    Opm::LinearSolverIstl linsolver(parameters);
-#else
-    Opm::LinearSolverUmfpack linsolver;
-#endif // EXPERIMENT_ISTL
+    Opm::LinearSolverFactory linsolver(parameters);
+
+    // EXPERIMENT_ISTL
     IncompTpfa pressure_solver(*grid.c_grid(), incomp_properties.permeability(), 
                                gravity, linsolver,  wells.c_wells());