Correctly multuiply blocks and degree

python/dune/perftool/pdelab/driver/gridfunctionspace.py

@@ -123,6 +123,10 @@ def typedef_fem(element, name):
     df = type_domainfield()
     r = type_range()
     dim = get_dimension()
+    if get_option("blockstructured"):
+        include_file("dune/perftool/blockstructured/blockstructuredqkfem.hh", filetag="driver")
+        degree = element.degree() * get_option("number_of_blocks")
+        return "using {} = Dune::PDELab::BlockstructuredQkLocalFiniteElementMap<{}, {}, {}, {}>;".format(name, gv, df, r, degree)
     if isQk(element):
         include_file("dune/pdelab/finiteelementmap/qkfem.hh", filetag="driver")
         return "using {} = Dune::PDELab::QkLocalFiniteElementMap<{}, {}, {}, {}>;".format(name, gv, df, r, element.degree())