Made the solve() methods const.

opm/core/linalg/LinearSolverInterface.cpp

@@ -35,7 +35,7 @@ namespace Opm
     LinearSolverInterface::LinearSolverReport
     LinearSolverInterface::solve(const CSRMatrix* A,
                                  const double* rhs,
-                                 double* solution)
+                                 double* solution) const
     {
         return solve(A->m, A->nnz, A->ia, A->ja, A->sa, rhs, solution);
     }

opm/core/linalg/LinearSolverInterface.hpp

@@ -52,7 +52,7 @@ namespace Opm
         /// Note: this method is a convenience method that calls the virtual solve() method.
         LinearSolverReport solve(const CSRMatrix* A,
                                  const double* rhs,
-                                 double* solution);
+                                 double* solution) const;
 
         /// Solve a linear system, with a matrix given in compressed sparse row format.
         /// \param[in] size        # of rows in matrix
@@ -69,7 +69,7 @@ namespace Opm
                                          const int* ja,
                                          const double* sa,
                                          const double* rhs,
-                                         double* solution) = 0;
+                                         double* solution) const = 0;
 
     };
 

opm/core/linalg/LinearSolverIstl.cpp

@@ -109,7 +109,7 @@ namespace Opm
                             const int* ja,
                             const double* sa,
                             const double* rhs,
-                            double* solution)
+                            double* solution) const
     {
         // Build Istl structures from input.
         // System matrix

opm/core/linalg/LinearSolverIstl.hpp

@@ -70,7 +70,7 @@ namespace Opm
                                          const int* ja,
                                          const double* sa,
                                          const double* rhs,
-                                         double* solution);
+                                         double* solution) const;
     private:
         double linsolver_residual_tolerance_;
         int linsolver_verbosity_;

opm/core/linalg/LinearSolverUmfpack.cpp

@@ -45,7 +45,7 @@ namespace Opm
                                const int* ja,
                                const double* sa,
                                const double* rhs,
-                               double* solution)
+                               double* solution) const
     {
         CSRMatrix A  = {
             size,

opm/core/linalg/LinearSolverUmfpack.hpp

@@ -76,7 +76,7 @@ namespace Opm
                                          const int* ja,
                                          const double* sa,
                                          const double* rhs,
-                                         double* solution);
+                                         double* solution) const;
     };