diff --git a/python/dune/codegen/blockstructured/basis.py b/python/dune/codegen/blockstructured/basis.py
index 073c9de35b1be48a4ea1e91d915922538f74704a..ea90582e5f991f94d5044ed47e9711f060b730e8 100644
--- a/python/dune/codegen/blockstructured/basis.py
+++ b/python/dune/codegen/blockstructured/basis.py
@@ -177,52 +177,3 @@ def name_localbasis(leaf_element):
     globalarg(name)
     define_localbasis(leaf_element, name)
     return name
-
-
-# @kernel_cached
-# def evaluate_basis(leaf_element, name, restriction):
-#     temporary_variable(name, shape=((leaf_element.degree() + 1)**world_dimension(), 1),
-#                        decl_method=declare_cache_temporary(leaf_element, restriction, 'Function'))
-#     cache = name_localbasis_cache(leaf_element)
-#     qp = pymbolic_quadrature_position_in_cell(restriction)
-#     localbasis = name_localbasis(leaf_element)
-#     instruction(inames=get_backend("quad_inames")(),
-#                 code='{} = {}.evaluateFunction({}, {});'.format(name, cache, str(qp), localbasis),
-#                 assignees=frozenset({name}),
-#                 read_variables=frozenset({get_pymbolic_basename(qp)}),
-#                 )
-#
-#
-# def pymbolic_basis(leaf_element, restriction, number, context=''):
-#     assert not isinstance(leaf_element, MixedElement)
-#     name = "phi_{}".format(FEM_name_mangling(leaf_element))
-#     name = restricted_name(name, restriction)
-#     evaluate_basis(leaf_element, name, restriction)
-#     inames = lfs_inames(leaf_element, restriction, number, context=context)
-#
-#     return prim.Subscript(prim.Variable(name), (tensor_index_to_sequential_index(inames, leaf_element.degree() + 1), 0))
-
-#
-# @kernel_cached
-# def evaluate_reference_gradient(leaf_element, name, restriction):
-#     from pudb import set_trace; set_trace()
-#     temporary_variable(name, shape=((leaf_element.degree() + 1)**world_dimension(), 1, world_dimension()),
-#                        decl_method=declare_cache_temporary(leaf_element, restriction, 'Jacobian'))
-#     cache = name_localbasis_cache(leaf_element)
-#     qp = pymbolic_quadrature_position_in_cell(restriction)
-#     localbasis = name_localbasis(leaf_element)
-#     instruction(inames=get_backend("quad_inames")(),
-#                 code='{} = {}.evaluateJacobian({}, {});'.format(name, cache, str(qp), localbasis),
-#                 assignees=frozenset({name}),
-#                 read_variables=frozenset({get_pymbolic_basename(qp)}),
-#                 )
-#
-#
-# def pymbolic_reference_gradient(leaf_element, restriction, number, context=''):
-#     assert not isinstance(leaf_element, MixedElement)
-#     name = "js_{}".format(FEM_name_mangling(leaf_element))
-#     name = restricted_name(name, restriction)
-#     evaluate_reference_gradient(leaf_element, name, restriction)
-#     inames = lfs_inames(leaf_element, restriction, number, context=context)
-#
-#     return prim.Subscript(prim.Variable(name), (tensor_index_to_sequential_index(inames, leaf_element.degree() + 1), 0))