From e02429f5dc74b60b83d31936cbabc77604a761aa Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= <atgeirr@sintef.no>
Date: Mon, 23 Apr 2012 11:07:43 +0200
Subject: [PATCH 1/2] Using LinearSolverFactory now.

---
 examples/polymer_reorder.cpp | 13 +++----------
 1 file changed, 3 insertions(+), 10 deletions(-)

diff --git a/examples/polymer_reorder.cpp b/examples/polymer_reorder.cpp
index 9556f2f93..eb5358470 100644
--- a/examples/polymer_reorder.cpp
+++ b/examples/polymer_reorder.cpp
@@ -42,11 +42,7 @@
 #include <opm/core/fluid/IncompPropertiesFromDeck.hpp>
 #include <opm/core/fluid/RockCompressibility.hpp>
 
-#include <opm/core/linalg/LinearSolverUmfpack.hpp>
-// #define EXPERIMENT_ISTL
-#ifdef EXPERIMENT_ISTL
-#include <opm/core/linalg/LinearSolverIstl.hpp>
-#endif
+#include <opm/core/linalg/LinearSolverFactory.hpp>
 
 #include <opm/core/transport/transport_source.h>
 #include <opm/core/transport/CSRMatrixUmfpackSolver.hpp>
@@ -483,12 +479,9 @@ main(int argc, char** argv)
     }
 
     // Solvers init.
+    // Linear solver.
+    Opm::LinearSolverFactory linsolver(param);
     // Pressure solver.
-#ifdef EXPERIMENT_ISTL
-    Opm::LinearSolverIstl linsolver(param);
-#else
-    Opm::LinearSolverUmfpack linsolver;
-#endif // EXPERIMENT_ISTL
     const double *grav = use_gravity ? &gravity[0] : 0;
     Opm::IncompTpfa psolver(*grid->c_grid(), props->permeability(), grav, linsolver);
     // Reordering solver.

From 4ae13394603ba84d6412c41db8e15df9cf02e624 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= <atgeirr@sintef.no>
Date: Mon, 23 Apr 2012 11:09:09 +0200
Subject: [PATCH 2/2] Catch exceptions from boost::create_directories() to give
 useful error message.

---
 examples/polymer_reorder.cpp | 7 ++++++-
 1 file changed, 6 insertions(+), 1 deletion(-)

diff --git a/examples/polymer_reorder.cpp b/examples/polymer_reorder.cpp
index eb5358470..ac6593c82 100644
--- a/examples/polymer_reorder.cpp
+++ b/examples/polymer_reorder.cpp
@@ -307,7 +307,12 @@ main(int argc, char** argv)
         output_dir = param.getDefault("output_dir", std::string("output"));
         // Ensure that output dir exists
         boost::filesystem::path fpath(output_dir);
-        create_directories(fpath);
+        try {
+            create_directories(fpath);
+        }
+        catch (...) {
+            THROW("Creating directories failed: " << fpath);
+        }
         output_interval = param.getDefault("output_interval", output_interval);
     }
     const int num_transport_substeps = param.getDefault("num_transport_substeps", 1);