Add a block assembler and solver for OPMPressure's CSRMatrix.

Solver based on UMFPACK.

src/CSRMatrixUmfpackSolver.hpp

@@ -0,0 +1,60 @@
+/*===========================================================================
+//
+// File: CSRMatrixUmfpackSolver.hpp
+//
+// Created: 2011-10-03 17:27:26+0200
+//
+// Authors: Ingeborg S. Ligaarden <Ingeborg.Ligaarden@sintef.no>
+//          Jostein R. Natvig     <Jostein.R.Natvig@sintef.no>
+//          Halvor M. Nilsen      <HalvorMoll.Nilsen@sintef.no>
+//          Atgeirr F. Rasmussen  <atgeirr@sintef.no>
+//          Bård Skaflestad       <Bard.Skaflestad@sintef.no>
+//
+//==========================================================================*/
+
+
+/*
+  Copyright 2011 SINTEF ICT, Applied Mathematics.
+  Copyright 2011 Statoil ASA.
+
+  This file is part of the Open Porous Media Project (OPM).
+
+  OPM is free software: you can redistribute it and/or modify
+  it under the terms of the GNU General Public License as published by
+  the Free Software Foundation, either version 3 of the License, or
+  (at your option) any later version.
+
+  OPM is distributed in the hope that it will be useful,
+  but WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+  GNU General Public License for more details.
+
+  You should have received a copy of the GNU General Public License
+  along with OPM.  If not, see <http://www.gnu.org/licenses/>.
+*/
+
+#ifndef OPM_CSRMATRIXUMFPACKSOLVER_HPP_HEADER
+#define OPM_CSRMATRIXUMFPACKSOLVER_HPP_HEADER
+
+#include <suitesparse/umfpack.h>
+
+#include <sparse_sys.h>
+#include <call_umfpack.h>
+
+namespace Opm {
+    namespace ImplicitTransportLinAlgSupport {
+        class CSRMatrixUmfpackSolver {
+        public:
+            template <class Vector>
+            void
+            solve(const struct CSRMatrix& A,
+                  const Vector&           b,
+                  Vector&                 x) {
+
+                call_UMFPACK(&A, &b[0], &x[0]);
+            }
+        };
+    }
+}
+
+#endif  /* OPM_CSRMATRIXUMFPACKSOLVER_HPP_HEADER */