diff --git a/python/dune/codegen/generation/mixins.py b/python/dune/codegen/generation/mixins.py
index fbde379c9d1c58a7b262d31f7aaa2681b378f320..06141dbe873f53d4c4644ee7863bb3fca9fad90b 100644
--- a/python/dune/codegen/generation/mixins.py
+++ b/python/dune/codegen/generation/mixins.py
@@ -12,7 +12,7 @@ def mixin_base(mixintype, name):
     def _dec(cls):
         _mixin_registry[mixintype][name] = cls
         return cls
-    
+
     return _dec
 
 
diff --git a/python/dune/codegen/pdelab/__init__.py b/python/dune/codegen/pdelab/__init__.py
index 4b999bd0bdcb65c2a7d9aa80059d2c65d61ecbcc..10175f38bb939a11cd48b3fc40206a8993645c38 100644
--- a/python/dune/codegen/pdelab/__init__.py
+++ b/python/dune/codegen/pdelab/__init__.py
@@ -3,6 +3,7 @@
 # Trigger some imports that are needed to have all backend implementations visible
 # to the selection mechanisms
 
+
 class PDELabInterface(object):
     def __init__(self):
         # The visitor instance will be registered by its init method
diff --git a/python/dune/codegen/pdelab/basis.py b/python/dune/codegen/pdelab/basis.py
index ed5469ce1b9b690479068400f38fc81806637fc8..583c3738f3e450f4839e852291a2ee868eb7b000 100644
--- a/python/dune/codegen/pdelab/basis.py
+++ b/python/dune/codegen/pdelab/basis.py
@@ -188,7 +188,7 @@ class GenericBasisMixin(BasisMixinBase):
         basis = self.implement_basis(sub_element, restriction, 0, context='trial')
         basisindex = get_pymbolic_indices(basis)[:-1]
 
-        #TODO get rid ot this!
+        # TODO get rid ot this!
         if get_form_option("blockstructured"):
             from dune.codegen.blockstructured.argument import pymbolic_coefficient
             coeff = pymbolic_coefficient(container, lfs, sub_element, basisindex)
diff --git a/python/dune/codegen/pdelab/geometry.py b/python/dune/codegen/pdelab/geometry.py
index 5bdc2fac640e5c6eb3d2363ad2d972c6560cafd2..4ade8a2284cdfff27120c396584982c675ca44d5 100644
--- a/python/dune/codegen/pdelab/geometry.py
+++ b/python/dune/codegen/pdelab/geometry.py
@@ -146,7 +146,7 @@ class GenericPDELabGeometryMixin(GeometryMixinBase):
         name = restricted_name("jit", restriction)
         self.define_jacobian_inverse_transposed(name, restriction)
 
-        return prim.Subscript(prim.Variable(name), (j,i))
+        return prim.Subscript(prim.Variable(name), (j, i))
 
     def define_jacobian_inverse_transposed(self, name, restriction):
         dim = world_dimension()
@@ -231,7 +231,6 @@ class AxiparallelGeometryMixin(GenericPDELabGeometryMixin):
 
         return jac
 
-
     def facet_area(self, o):
         # This is not 100% correct, but in practice structured grid implementations are not
         # embedded into higher dimensional world space.
diff --git a/python/dune/codegen/sumfact/__init__.py b/python/dune/codegen/sumfact/__init__.py
index fa05f8a0bd82b53ae35089e37a5e1850cf128e0e..ea07249ef76da9b53a16423a2a473b69e62aadf8 100644
--- a/python/dune/codegen/sumfact/__init__.py
+++ b/python/dune/codegen/sumfact/__init__.py
@@ -6,4 +6,4 @@ from dune.codegen.pdelab import PDELabInterface
 
 
 class SumFactInterface(PDELabInterface):
-    pass
\ No newline at end of file
+    pass
diff --git a/python/dune/codegen/sumfact/geometry.py b/python/dune/codegen/sumfact/geometry.py
index ade47c57918c025463109327531d79975d4c4f21..aa3963405c2b2915e2b7a1a6bc1792e5b71094cc 100644
--- a/python/dune/codegen/sumfact/geometry.py
+++ b/python/dune/codegen/sumfact/geometry.py
@@ -66,7 +66,7 @@ class SumfactMultiLinearGeometryMixin(GenericPDELabGeometryMixin):
         name = restricted_name(self._jacobian_inverse(), restriction)
         self.define_jacobian_inverse(name, restriction)
 
-        return prim.Subscript(prim.Variable(name), (i,j))
+        return prim.Subscript(prim.Variable(name), (i, j))
 
     def _jacobian_determinant(self):
         return "detjac"
@@ -312,6 +312,7 @@ class SumfactEqudistantGeometryMixin(EquidistantGeometryMixin, SumfactAxiParalle
         self.indices = None
         return prim.Subscript(prim.Variable(lowcorner), (index,)) + x * prim.Subscript(prim.Variable(meshwidth), (index,))
 
+
 @iname
 def global_corner_iname(restriction):
     name = get_counted_variable(restricted_name("global_corneriname", restriction))
@@ -471,4 +472,4 @@ def _name_jacobian(i, j, restriction, visitor):
 
 
 def normalize(expr, dim):
-    return prim.Call(prim.Variable("sqrt"), (prim.Sum(tuple(expr[i]*expr[i] for i in range(dim))),))
+    return prim.Call(prim.Variable("sqrt"), (prim.Sum(tuple(expr[i] * expr[i] for i in range(dim))),))