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

Laplace DG is working!

parent 28381e6e
No related branches found
No related tags found
No related merge requests found
......@@ -260,7 +260,7 @@ def typedef_fem(expr, name):
if isQuadrilateral(expr):
include_file("dune/pdelab/finiteelementmap/qkdg.hh", filetag="driver")
# TODO allow switching the basis here!
return "typedef Dune::PDELab::QkDGLocalFiniteElementMap<{}, {}, {}, {}> {}".format(df, r, expr._degree, dim, name)
return "typedef Dune::PDELab::QkDGLocalFiniteElementMap<{}, {}, {}, {}> {};".format(df, r, expr._degree, dim, name)
if isSimplical(expr):
include_file("dune/pdelab/finiteelementmap/opbfem.hh", filetag="driver")
return "typedef Dune::PDELab::OPBLocalFiniteElementMap<{}, {}, {}, {}, Dune::GeometryType::simplex> {};".format(df, r, expr._degree, dim, name)
......
......@@ -216,7 +216,7 @@ def apply_in_cell_transformation(name, local, restriction):
def name_in_cell_coordinates(local, basename, restriction):
name = "{}_in_{}side".format(basename, "in" if restriction is Restriction.NEGATIVE else "out")
temporary_variable(name, shape=(name_dimension(),), shape_impl=("fv",))
apply_in_cell_transformation(name, local, restriction=Restriction.NEGATIVE)
apply_in_cell_transformation(name, local, restriction)
return name
......
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