Skip to content
Snippets Groups Projects
Commit 3560297a authored by Dominic Kempf's avatar Dominic Kempf
Browse files

Fixup 1-dimensional case

IT is both simplicial and quadrilateral...
parent 9e195534
No related branches found
No related tags found
No related merge requests found
...@@ -154,20 +154,20 @@ def typedef_fem(element, name): ...@@ -154,20 +154,20 @@ def typedef_fem(element, name):
# Choose the correct finite element implementation # Choose the correct finite element implementation
if short == "CG": if short == "CG":
if isQuadrilateral(cell): if isSimplical(cell):
if dim in (2, 3) and degree < 3:
include_file("dune/pdelab/finiteelementmap/qkfem.hh", filetag="driver")
return "using {} = Dune::PDELab::QkLocalFiniteElementMap<{}, {}, {}, {}>;" \
.format(name, gv, df, r, degree)
else:
raise CodegenUnsupportedFiniteElementError(element)
elif isSimplical(cell):
if dim in (1, 2, 3): if dim in (1, 2, 3):
include_file("dune/pdelab/finiteelementmap/pkfem.hh", filetag="driver") include_file("dune/pdelab/finiteelementmap/pkfem.hh", filetag="driver")
return "using {} = Dune::PDELab::PkLocalFiniteElementMap<{}, {}, {}, {}>;" \ return "using {} = Dune::PDELab::PkLocalFiniteElementMap<{}, {}, {}, {}>;" \
.format(name, gv, df, r, degree) .format(name, gv, df, r, degree)
else: else:
raise CodegenUnsupportedFiniteElementError(element) raise CodegenUnsupportedFiniteElementError(element)
elif isQuadrilateral(cell):
if dim in (2, 3) and degree < 3:
include_file("dune/pdelab/finiteelementmap/qkfem.hh", filetag="driver")
return "using {} = Dune::PDELab::QkLocalFiniteElementMap<{}, {}, {}, {}>;" \
.format(name, gv, df, r, degree)
else:
raise CodegenUnsupportedFiniteElementError(element)
else: else:
raise CodegenUnsupportedFiniteElementError(element) raise CodegenUnsupportedFiniteElementError(element)
elif short == "DG": elif short == "DG":
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment