diff --git a/python/dune/perftool/pdelab/driver.py b/python/dune/perftool/pdelab/driver.py
index e11e659fc0c8a4f2e3045cb13485c492c626bb2b..26de3b67de909ee8156d33aefafcaff3727c718d 100644
--- a/python/dune/perftool/pdelab/driver.py
+++ b/python/dune/perftool/pdelab/driver.py
@@ -1783,7 +1783,7 @@ def generate_driver(formdatas, data):
         # In case of operator conunting we only assemble the matrix and evaluate the residual
         #assemble_matrix_timer()
         evaluate_residual_timer()
-#        apply_jacobian_timer()
+        apply_jacobian_timer()
     elif is_stationary():
         # We could also use solve if we are not interested in visualization
         vtkoutput()
diff --git a/python/dune/perftool/sumfact/__init__.py b/python/dune/perftool/sumfact/__init__.py
index 520fa7d1c00ae4fc80ce0f6dcb8bc927439d6272..60b120ce95e9650b9a358632f6d3bad7249a4d25 100644
--- a/python/dune/perftool/sumfact/__init__.py
+++ b/python/dune/perftool/sumfact/__init__.py
@@ -37,6 +37,16 @@ class SumFactInterface(PDELabInterface):
         visitor.indices = indices
         return ret
 
+    def pymbolic_apply_function_gradient(self, element, restriction, component, visitor=None):
+        ret, indices = pymbolic_coefficient_gradient(element, restriction, component, name_applycontainer, visitor)
+        visitor.indices = indices
+        return ret
+
+    def pymbolic_apply_function(self, element, restriction, component, visitor=None):
+        ret, indices = pymbolic_coefficient(element, restriction, component, name_applycontainer, visitor)
+        visitor.indices = indices
+        return ret
+
     def quadrature_inames(self):
         return quadrature_inames()