Port all the hyperbolic tests

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

@@ -118,7 +118,9 @@ def define_timesteppingmethod(name):
     if explicit:
         return "{} {};".format(tsm_type, name)
     else:
-        return "{} {}(1.0);".format(tsm_type, name)
+        ini = name_initree()
+        return "{} {}({}.get<double>(\"instat.theta\",1.0));".format(tsm_type, name, ini)
+
 
 
 def name_timesteppingmethod():

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

@@ -103,7 +103,11 @@ def define_boundary_lambda(name, boundary):
             from dune.perftool.ufl.visitor import UFL2LoopyVisitor
             from dune.perftool.pdelab import PDELabInterface
             visitor = UFL2LoopyVisitor(PDELabInterface(), "exterior_facet", {})
+            from pymbolic.mapper.c_code import CCodeMapper
+            ccm = CCodeMapper()
+            expr = visitor(boundary)
+
             from dune.perftool.loopy.target import numpy_to_cpp_dtype
             return "auto {} = [&](const auto& x){{ return ({}){}; }};".format(name,
                                                                               numpy_to_cpp_dtype("float64"),
-                                                                              visitor(boundary))
+                                                                              ccm(expr))

test/hyperbolic/linearacoustics.mini

@@ -18,7 +18,3 @@ extension = vtu
 [formcompiler]
 numerical_jacobian = 1, 0 | expand diff
 explicit_time_stepping = 0
-# This tests that all mass is transported out of the domain.
-# While this is not the best of tests, it is something easily checked for.
-#exact_solution_expression = g
-#compare_l2errorsquared = 1e-10

test/hyperbolic/linearacoustics.ufl

 cell = quadrilateral
 x = SpatialCoordinate(cell)
 
-#V = FiniteElement("DG", cell, 1)
-
 MV = VectorElement("DG", cell=cell, degree=1, dim=3)
 
 n = FacetNormal(cell)('+')

test/hyperbolic/lineartransport.mini

@@ -20,7 +20,4 @@ extension = vtu
 [formcompiler]
 numerical_jacobian = 1, 0 | expand diff
 explicit_time_stepping = 0, 1 | expand scheme
-# This tests that all mass is transported out of the domain.
-# While this is not the best of tests, it is something easily checked for.
-exact_solution_expression = g
 compare_l2errorsquared = 1e-10
\ No newline at end of file

test/hyperbolic/lineartransport.ufl

@@ -2,13 +2,12 @@ cell = triangle
 
 x = SpatialCoordinate(cell)
 initial = conditional(And(And(x[0] > 0.1, x[0] < 0.2), And(x[1] > 0.1, x[1] < 0.2)), 1, 0)
-g = 0.0
 dirichlet = conditional(x[0] + x[1] < 1.0, 1, 0)
 
 dsd = ds(subdomain_data=dirichlet, subdomain_id=1)
 dso = ds(subdomain_data=dirichlet, subdomain_id=0)
 
-V = FiniteElement("DG", cell, 1, dirichlet_constraints=dirichlet, dirichlet_expression=initial)
+V = FiniteElement("DG", cell, 1)
 
 u = TrialFunction(V)
 v = TestFunction(V)
@@ -26,4 +25,7 @@ r = -1.*u*inner(beta, grad(v))*dx \
   + inner(beta, n)*u*v*dso \
   + numerical_flux(n, g, u('-'))*v*dsd
 
-forms = [mass, r]
+forms = [mass, r]
\ No newline at end of file
+is_dirichlet = dirichlet
+dirichlet_expression = initial
+exact_solution = 0
\ No newline at end of file

test/sumfact/hyperbolic/linearacoustics.mini

@@ -19,7 +19,3 @@ extension = vtu
 numerical_jacobian = 1, 0 | expand diff
 explicit_time_stepping = 0
 sumfact = 1
-# This tests that all mass is transported out of the domain.
-# While this is not the best of tests, it is something easily checked for.
-#exact_solution_expression = g
-#compare_l2errorsquared = 1e-10

test/sumfact/hyperbolic/lineartransport.mini

@@ -19,5 +19,4 @@ sumfact = 1
 #fastdg = 1
 # This tests that all mass is transported out of the domain.
 # While this is not the best of tests, it is something easily checked for.
-exact_solution_expression = g
 compare_l2errorsquared = 1e-10

test/sumfact/hyperbolic/lineartransport.ufl

@@ -2,13 +2,12 @@ cell = quadrilateral
 
 x = SpatialCoordinate(cell)
 initial = conditional(And(And(x[0] > 0.1, x[0] < 0.2), And(x[1] > 0.1, x[1] < 0.2)), 1, 0)
-g = 0.0
 dirichlet = conditional(x[0] + x[1] < 1.0, 1, 0)
 
 dsd = ds(subdomain_data=dirichlet, subdomain_id=1)
 dso = ds(subdomain_data=dirichlet, subdomain_id=0)
 
-V = FiniteElement("DG", cell, 3, dirichlet_constraints=dirichlet, dirichlet_expression=initial)
+V = FiniteElement("DG", cell, 3)
 
 u = TrialFunction(V)
 v = TestFunction(V)
@@ -26,4 +25,7 @@ r = -1.*u*inner(beta, grad(v))*dx \
   + inner(beta, n)*u*v*dso \
   + numerical_flux(n, g, u('-'))*v*dsd
 
-forms = [mass, r]
+forms = [mass, r]
\ No newline at end of file
+exact_solution = 0.0
+is_dirichlet = dirichlet
+dirichlet_expression = initial
\ No newline at end of file