diff --git a/python/dune/codegen/sumfact/quadrature.py b/python/dune/codegen/sumfact/quadrature.py
index b68da208360fe12aa7b85894759b9dfa53f4136f..92a598d958096601afe9f74b59533a84fe0f2d80 100644
--- a/python/dune/codegen/sumfact/quadrature.py
+++ b/python/dune/codegen/sumfact/quadrature.py
@@ -179,63 +179,6 @@ def constructor_quadrature_inames(dim):
     return tuple(constructor_quad_iname(name, d, local_qps_per_dir[d]) for d in range(local_dimension()))
 
 
-def define_recursive_quadrature_weight(visitor, name, direction):
-    info = visitor.current_info[1]
-    if info is None:
-        element = None
-    else:
-        element = info.element
-    iname = quadrature_inames(element)[direction]
-    temporary_variable(name, shape=(), shape_impl=())
-
-    qps_per_dir = quadrature_points_per_direction()
-    qp_bound = qps_per_dir[direction]
-    instruction(expression=Product((recursive_quadrature_weight(direction + 1),
-                                    Subscript(Variable(name_oned_quadrature_weights(qp_bound)),
-                                              (Variable(iname),)
-                                              ))
-                                   ),
-                assignee=Variable(name),
-                forced_iname_deps=frozenset(quadrature_inames()[direction:]),
-                forced_iname_deps_is_final=True,
-                tags=frozenset({"quad"}),
-                )
-
-
-def recursive_quadrature_weight(visitor, direction=0):
-    if direction == local_dimension():
-        return pymbolic_base_weight()
-    else:
-        name = 'weight_{}'.format(direction)
-        define_recursive_quadrature_weight(visitor, name, direction)
-        return Variable(name)
-
-
-def define_quadrature_position(name, local_index):
-    # TODO: This whole function looks quite suspicious:
-    # - qps_per_dir not used
-    # - probably rename index to local_index
-    #
-    # This whole function is probably not used anymore since we use sum
-    # factorization for evaluation of global coordinate, jacobian inverse
-    # transposed and determinant of the jacobian of the geometry mapping.
-    #
-    # Needs to be fixed. Or removed. Added to my to do list ;)
-    assert False
-
-    qps_per_dir = quadrature_points_per_direction()
-    local_qps_per_dir = local_quadrature_points_per_direction()
-    qp_bound = local_qps_per_dir[index]
-    expression = Subscript(Variable(name_oned_quadrature_points(qp_bound)),
-                           (Variable(quadrature_inames()[index]),))
-    instruction(expression=expression,
-                assignee=Subscript(Variable(name), (index,)),
-                forced_iname_deps=frozenset(quadrature_inames(element)),
-                forced_iname_deps_is_final=True,
-                tags=frozenset({"quad"}),
-                )
-
-
 def additional_inames(visitor):
     """Return inames for iterating over ansatz space in the case of jacobians
     """