diff --git a/python/dune/perftool/pdelab/basis.py b/python/dune/perftool/pdelab/basis.py
index 3fd62f3c4118f11635bc97d2c562e942c665ff7b..2bda1822594203be1e406e6452154fb13fd8f8f5 100644
--- a/python/dune/perftool/pdelab/basis.py
+++ b/python/dune/perftool/pdelab/basis.py
@@ -170,7 +170,7 @@ def evaluate_coefficient(visitor, element, name, container, restriction, index):
 
     if get_option("blockstructured"):
         from dune.perftool.blockstructured.argument import pymbolic_coefficient
-        coeff = pymbolic_coefficient(container, lfs, element, basisindex)
+        coeff = pymbolic_coefficient(container, lfs, sub_element, basisindex)
     else:
         from dune.perftool.pdelab.argument import pymbolic_coefficient
         coeff = pymbolic_coefficient(container, lfs, basisindex)
@@ -209,7 +209,7 @@ def evaluate_coefficient_gradient(visitor, element, name, container, restriction
 
     if get_option("blockstructured"):
         from dune.perftool.blockstructured.argument import pymbolic_coefficient
-        coeff = pymbolic_coefficient(container, lfs, element, basisindex)
+        coeff = pymbolic_coefficient(container, lfs, sub_element, basisindex)
     else:
         from dune.perftool.pdelab.argument import pymbolic_coefficient
         coeff = pymbolic_coefficient(container, lfs, basisindex)
diff --git a/test/blockstructured/stokes/stokes.mini b/test/blockstructured/stokes/stokes.mini
index 9a4903b89041ce50ecbbfabc0fbf2264bef53d1d..c44cc189861c3a1e23963ff21fd56568fdaf4173 100644
--- a/test/blockstructured/stokes/stokes.mini
+++ b/test/blockstructured/stokes/stokes.mini
@@ -13,7 +13,4 @@ extension = vtu
 numerical_jacobian = 0, 1 | expand num
 compare_l2errorsquared = 1e-10
 blockstructured = 1
-number_of_blocks = 3
-
-# Ignore the numdiff case until the associated issue is fixed!
-{__exec_suffix} == numdiff | exclude
\ No newline at end of file
+number_of_blocks = 3
\ No newline at end of file