diff --git a/python/dune/perftool/pdelab/driver/vtk.py b/python/dune/perftool/pdelab/driver/vtk.py
index 2f68140e15c68e572616cd30d2e8a288ef61ebda..15c39c3d22f2a7542c0c08e59cc8ac1c9838ab31 100644
--- a/python/dune/perftool/pdelab/driver/vtk.py
+++ b/python/dune/perftool/pdelab/driver/vtk.py
@@ -4,8 +4,10 @@ from dune.perftool.generation import (include_file,
 from dune.perftool.pdelab.driver import (get_formdata,
                                          get_trial_element,
                                          name_initree,
+                                         preprocess_leaf_data,
                                          )
 from dune.perftool.pdelab.driver.gridfunctionspace import (name_leafview,
+                                                           name_gfs,
                                                            name_trial_gfs,
                                                            type_leafview,
                                                            )
@@ -89,7 +91,6 @@ def name_predicate():
     return "predicate"
 
 
-
 @preamble
 def typedef_vtk_sequence_writer(name):
     include_file("dune/grid/io/file/vtk/vtksequencewriter.hh", filetag="driver")
@@ -122,7 +123,7 @@ def visualize_initial_condition():
     vtkwriter = name_vtk_sequence_writer()
     element = get_trial_element()
     define_gfs_names(element)
-    gfs = name_gfs(element)
+    gfs = name_trial_gfs()
     vector = name_vector(get_formdata())
     predicate = name_predicate()
     from dune.perftool.pdelab.driver.instationary import name_time
diff --git a/test/hyperbolic/linearacoustics.ufl b/test/hyperbolic/linearacoustics.ufl
index 97627f232484863617400eec74cdedfaa99bcae9..3fd12b2caa109755f79074d439362f42b9b98da0 100644
--- a/test/hyperbolic/linearacoustics.ufl
+++ b/test/hyperbolic/linearacoustics.ufl
@@ -1,7 +1,8 @@
 cell = quadrilateral
 x = SpatialCoordinate(cell)
 
-MV = VectorElement("DG", cell=cell, degree=1, dim=3)
+V = FiniteElement("DG", cell, 1)
+MV = MixedElement(V, V, V)
 
 n = FacetNormal(cell)('+')
 
@@ -20,3 +21,4 @@ r = -1.*inner(flux, grad(v))*dx \
   + inner(u, v)*ds
 
 forms = [mass, r]
+dirichlet_expression =  sin(pi*(x[0]-0.45)/0.1)*sin(pi*(x[0]-0.45)/0.1)*sin(pi*(x[1]-0.3)/0.1)*sin(pi*(x[1]-0.3)/0.1), 0.0, 0.0