diff --git a/src/compute_state.cc b/src/compute_state.cc
index 94e1a0d7156b1772de7015c173add3ca87fb3f86..24e420a86155e85de066c0ae9b371ae1b755fb99 100644
--- a/src/compute_state.cc
+++ b/src/compute_state.cc
@@ -15,8 +15,6 @@ class DecayingExponential {
   typedef Dune::FieldVector<double, 1> VectorType;
   typedef Dune::FieldMatrix<double, 1, 1> MatrixType;
 
-  double operator()(VectorType const &u) const { return std::exp(-u[0]); }
-
   void directionalSubDiff(VectorType const &u, VectorType const &v,
                           Interval<double> &D) const {
     D[0] = D[1] = v[0] * (-std::exp(-u[0]));