diff --git a/opm/autodiff/FullyImplicitBlackoilSolver.hpp b/opm/autodiff/FullyImplicitBlackoilSolver.hpp
index ea45b2fa6..5fd4dc046 100644
--- a/opm/autodiff/FullyImplicitBlackoilSolver.hpp
+++ b/opm/autodiff/FullyImplicitBlackoilSolver.hpp
@@ -289,9 +289,6 @@ namespace Opm {
 
         void applyThresholdPressures(ADB& dp);
 
-        double
-        residualNorm() const;
-
         /// \brief Compute the residual norms of the mass balance for each phase,
         /// the well flux, and the well equation.
         /// \return a vector that contains for each phase the norm of the mass balance
diff --git a/opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp b/opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
index 66a94f710..d92b4aed0 100644
--- a/opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
+++ b/opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
@@ -1722,34 +1722,6 @@ namespace detail {
         dp = keep_high_potential * (dp - threshold_modification);
     }
 
-
-
-
-
-    template<class T>
-    double
-    FullyImplicitBlackoilSolver<T>::residualNorm() const
-    {
-        double globalNorm = 0;
-        std::vector<ADB>::const_iterator quantityIt = residual_.material_balance_eq.begin();
-        const std::vector<ADB>::const_iterator endQuantityIt = residual_.material_balance_eq.end();
-        for (; quantityIt != endQuantityIt; ++quantityIt) {
-            const double quantityResid = (*quantityIt).value().matrix().norm();
-            if (!std::isfinite(quantityResid)) {
-                const int trouble_phase = quantityIt - residual_.material_balance_eq.begin();
-                OPM_THROW(Opm::NumericalProblem,
-                          "Encountered a non-finite residual in material balance equation "
-                          << trouble_phase);
-            }
-            globalNorm = std::max(globalNorm, quantityResid);
-        }
-        globalNorm = std::max(globalNorm, residual_.well_flux_eq.value().matrix().norm());
-        globalNorm = std::max(globalNorm, residual_.well_eq.value().matrix().norm());
-
-        return globalNorm;
-    }
-
-
     template<class T>
     std::vector<double>
     FullyImplicitBlackoilSolver<T>::computeResidualNorms() const