Use ufl_shape instead of shape

python/dune/codegen/pdelab/tensors.py

@@ -177,7 +177,7 @@ def pymbolic_matrix_inverse(o, visitor):
     expr = o.ufl_operands[0]
     name = name_assembled_tensor(expr, visitor)
 
-    if expr.shape[0] <= 3:
+    if expr.ufl_shape[0] <= 3:
         name = name_matrix_inverse(name, expr.ufl_shape, visitor)
     else:
         instruction(code="{}.invert();".format(name),