From c39a4e02a8cc492a250ce444e88db0313b34e46f Mon Sep 17 00:00:00 2001
From: Dominic Kempf <dominic.kempf@iwr.uni-heidelberg.de>
Date: Mon, 28 Aug 2017 10:19:46 +0200
Subject: [PATCH] Correctly multuiply blocks and degree
---
python/dune/perftool/pdelab/driver/gridfunctionspace.py | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/python/dune/perftool/pdelab/driver/gridfunctionspace.py b/python/dune/perftool/pdelab/driver/gridfunctionspace.py
index 67b7e64e..2d9b06ba 100644
--- a/python/dune/perftool/pdelab/driver/gridfunctionspace.py
+++ b/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())
--
GitLab