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

Port all scalar non-sumfact tests

parent 38689d42
No related branches found
No related tags found
No related merge requests found
...@@ -13,7 +13,6 @@ extension = vtu ...@@ -13,7 +13,6 @@ extension = vtu
[formcompiler] [formcompiler]
explicit_time_stepping = 0, 1 | expand scheme explicit_time_stepping = 0, 1 | expand scheme
exact_solution_expression = g
compare_l2errorsquared = 1e-7 compare_l2errorsquared = 1e-7
# Disable explicit tests for now # Disable explicit tests for now
......
f = Expression("Dune::FieldVector<double,2> c(0.5); c-= x; return 4.*(1.-c.two_norm2())*std::exp(-1.*c.two_norm2());") f = Expression("Dune::FieldVector<double,2> c(0.5); c-= x; return 4.*(1.-c.two_norm2())*std::exp(-1.*c.two_norm2());")
g = Expression("Dune::FieldVector<double,2> c(0.5); c-= x; return std::exp(-1.*c.two_norm2());") g = Expression("Dune::FieldVector<double,2> c(0.5); c-= x; return std::exp(-1.*c.two_norm2());")
V = FiniteElement("CG", "triangle", 1, dirichlet_expression=g) V = FiniteElement("CG", "triangle", 1)
u = TrialFunction(V) u = TrialFunction(V)
v = TestFunction(V) v = TestFunction(V)
mass = (u*v)*dx mass = (u*v)*dx
poisson = (inner(grad(u), grad(v)) - f*v)*dx poisson = (inner(grad(u), grad(v)) - f*v)*dx
forms = [mass, poisson] forms = [mass, poisson]
\ No newline at end of file dirichlet_expression = g
is_dirichlet = 1
exact_solution = g
\ No newline at end of file
...@@ -13,7 +13,6 @@ extension = vtu ...@@ -13,7 +13,6 @@ extension = vtu
[formcompiler] [formcompiler]
explicit_time_stepping = 0, 1 | expand scheme explicit_time_stepping = 0, 1 | expand scheme
exact_solution_expression = g
compare_l2errorsquared = 1e-7 compare_l2errorsquared = 1e-7
# Disable explicit tests for now # Disable explicit tests for now
......
...@@ -29,3 +29,6 @@ poisson = inner(grad(u), grad(v))*dx \ ...@@ -29,3 +29,6 @@ poisson = inner(grad(u), grad(v))*dx \
mass = (u*v)*dx mass = (u*v)*dx
forms = [mass, poisson] forms = [mass, poisson]
dirichlet_expression = g
is_dirichlet = 1
exact_solution = g
\ No newline at end of file
...@@ -11,5 +11,4 @@ reference = nonlinear_ref ...@@ -11,5 +11,4 @@ reference = nonlinear_ref
extension = vtu extension = vtu
[formcompiler] [formcompiler]
exact_solution_expression = g
compare_l2errorsquared = 6e-4 compare_l2errorsquared = 6e-4
\ No newline at end of file
f = Expression("return -4;") f = Expression("return -4;")
g = Expression("return x.two_norm2();") g = Expression("return x.two_norm2();")
V = FiniteElement("CG", "triangle", 1, dirichlet_expression=g) V = FiniteElement("CG", "triangle", 1)
u = TrialFunction(V) u = TrialFunction(V)
v = TestFunction(V) v = TestFunction(V)
r = (inner(grad(u), grad(v)) + u*u*v - f*v)*dx r = (inner(grad(u), grad(v)) + u*u*v - f*v)*dx
forms = [r] forms = [r]
exact_solution = g
dirichlet_expression = g
is_dirichlet = 1
\ No newline at end of file
...@@ -11,6 +11,4 @@ reference = nonlinear_dg ...@@ -11,6 +11,4 @@ reference = nonlinear_dg
extension = vtu extension = vtu
[formcompiler] [formcompiler]
exact_solution_expression = g
compare_dofs = 1e-1
compare_l2errorsquared = 6e-4 compare_l2errorsquared = 6e-4
...@@ -31,3 +31,4 @@ r = inner(grad(u), grad(v))*dx \ ...@@ -31,3 +31,4 @@ r = inner(grad(u), grad(v))*dx \
- gamma*g*v*ds - gamma*g*v*ds
forms = [r] forms = [r]
exact_solution = g
...@@ -13,7 +13,6 @@ extension = vtu ...@@ -13,7 +13,6 @@ extension = vtu
[formcompiler] [formcompiler]
numerical_jacobian = 1, 0 | expand num numerical_jacobian = 1, 0 | expand num
matrix_free = 1 matrix_free = 1
exact_solution_expression = g
compare_l2errorsquared = 5e-3 compare_l2errorsquared = 5e-3
# Disable numerical methods (not working in PDELab?) # Disable numerical methods (not working in PDELab?)
......
...@@ -13,5 +13,4 @@ extension = vtu ...@@ -13,5 +13,4 @@ extension = vtu
[formcompiler] [formcompiler]
numerical_jacobian = 1, 0 | expand num numerical_jacobian = 1, 0 | expand num
matrix_free = 1 matrix_free = 1
exact_solution_expression = g
compare_l2errorsquared = 6e-4 compare_l2errorsquared = 6e-4
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