diff --git a/opm/core/simulator/SimulatorTwophase.cpp b/opm/core/simulator/SimulatorTwophase.cpp
index 1a617a4b..71586468 100644
--- a/opm/core/simulator/SimulatorTwophase.cpp
+++ b/opm/core/simulator/SimulatorTwophase.cpp
@@ -65,7 +65,7 @@ namespace Opm
              const Wells* wells,
              const std::vector<double>& src,
              const FlowBoundaryConditions* bcs,
-             const LinearSolverInterface& linsolver,
+             LinearSolverInterface& linsolver,
              const double* gravity);
 
         void run(SimulatorTimer& timer,
@@ -115,7 +115,7 @@ namespace Opm
                                          const Wells* wells,
                                          const std::vector<double>& src,
                                          const FlowBoundaryConditions* bcs,
-                                         const LinearSolverInterface& linsolver,
+                                         LinearSolverInterface& linsolver,
                                          const double* gravity)
     {
         pimpl_.reset(new Impl(param, grid, props, rock_comp, wells, src, bcs, linsolver, gravity));
@@ -204,7 +204,7 @@ namespace Opm
                                   const Wells* wells,
                                   const std::vector<double>& src,
                                   const FlowBoundaryConditions* bcs,
-                                  const LinearSolverInterface& linsolver,
+                                  LinearSolverInterface& linsolver,
                                   const double* gravity)
         : grid_(grid),
           props_(props),
diff --git a/opm/core/simulator/SimulatorTwophase.hpp b/opm/core/simulator/SimulatorTwophase.hpp
index 053a3196..7b885a1b 100644
--- a/opm/core/simulator/SimulatorTwophase.hpp
+++ b/opm/core/simulator/SimulatorTwophase.hpp
@@ -71,7 +71,7 @@ namespace Opm
                          const Wells* wells,
                          const std::vector<double>& src,
                          const FlowBoundaryConditions* bcs,
-                         const LinearSolverInterface& linsolver,
+                         LinearSolverInterface& linsolver,
                          const double* gravity);
 
         /// Run the simulation.