From 31ddb2f098a5889f93bd141c5322507d2fa35754 Mon Sep 17 00:00:00 2001
From: Liu Ming <miliu@statoil.com>
Date: Wed, 26 Mar 2014 16:02:39 +0800
Subject: [PATCH] let shear vef return derivative.

---
 opm/polymer/PolymerProperties.cpp | 7 +++++++
 opm/polymer/PolymerProperties.hpp | 2 ++
 2 files changed, 9 insertions(+)

diff --git a/opm/polymer/PolymerProperties.cpp b/opm/polymer/PolymerProperties.cpp
index 6ac5f125e..c2e1c4438 100644
--- a/opm/polymer/PolymerProperties.cpp
+++ b/opm/polymer/PolymerProperties.cpp
@@ -71,6 +71,13 @@ namespace Opm
         return Opm::linearInterpolation(water_vel_vals_, shear_vrf_vals_, velocity);
     }
 
+    double
+    PolymerProperties::shearVrfWithDer(const double velocity, double& der) const
+    {
+        der =  Opm::linearInterpolationDerivative(water_vel_vals_, shear_vrf_vals_, velocity);
+        return Opm::linearInterpolation(water_vel_vals_, shear_vrf_vals_, velocity);
+    }
+
     double PolymerProperties::viscMult(double c) const
     {
         return Opm::linearInterpolation(c_vals_visc_, visc_mult_vals_, c);
diff --git a/opm/polymer/PolymerProperties.hpp b/opm/polymer/PolymerProperties.hpp
index cfc01a835..0482fd46b 100644
--- a/opm/polymer/PolymerProperties.hpp
+++ b/opm/polymer/PolymerProperties.hpp
@@ -168,6 +168,8 @@ namespace Opm
 
         double shearVrf(const double velocity) const;
 
+        double shearVrfWithDer(const double velocity, double& der) const;
+
         double viscMult(double c) const;
 
         double viscMultWithDer(double c, double* der) const;