Skip to content
Snippets Groups Projects
Commit 1755a04e authored by René Heß's avatar René Heß
Browse files

Fix poisson tests for dim=3 and dim=1

parent 97c8e54a
No related branches found
No related tags found
No related merge requests found
......@@ -3,7 +3,7 @@ cell = "interval"
x = SpatialCoordinate(cell)
c = (0.5-x[0])**2
g = exp(-1.*c)
f = 2*(1.-c)*g
f = 2*(1.-2*c)*g
V = FiniteElement("CG", cell, 1, dirichlet_expression=g)
u = TrialFunction(V)
......
......@@ -3,7 +3,7 @@ cell = "interval"
x = SpatialCoordinate(cell)
c = (0.5-x[0])**2
g = exp(-1.*c)
f = 4*(1.-c)*g
f = 2*(1.-2*c)*g
V = FiniteElement("DG", cell, 1)
......
......@@ -3,7 +3,7 @@ cell = "quadrilateral"
x = SpatialCoordinate(cell)
c = (0.5-x[0])**2 + (0.5-x[1])**2
g = exp(-1.*c)
f = 4*(1.-c)*g
f = 2*(2.-2*c)*g
V = FiniteElement("CG", cell, 1, dirichlet_expression=g)
u = TrialFunction(V)
......
......@@ -3,7 +3,7 @@ cell = "quadrilateral"
x = SpatialCoordinate(cell)
c = (0.5-x[0])**2 + (0.5-x[1])**2
g = exp(-1.*c)
f = 4*(1.-c)*g
f = 2*(2.-2*c)*g
V = FiniteElement("DG", cell, 1)
......
......@@ -3,7 +3,7 @@ cell = "hexahedron"
x = SpatialCoordinate(cell)
c = (0.5 - x[0])**2 + (0.5 - x[1])**2 + (0.5 - x[2])**2
g = exp(-1.*c)
f = 6*(1.-c)*g
f = 2*(3.-2*c)*g
V = FiniteElement("DG", cell, 1)
......
......@@ -3,7 +3,7 @@ cell = "quadrilateral"
x = SpatialCoordinate(cell)
c = (0.5-x[0])**2 + (0.5-x[1])**2
g = exp(-1.*c)
f = 4*(1.-c)*g
f = 2*(2.-2*c)*g
V = FiniteElement("CG", cell, degree, dirichlet_expression=g)
u = TrialFunction(V)
......
......@@ -3,7 +3,7 @@ cell = "hexahedron"
x = SpatialCoordinate(cell)
c = (0.5 - x[0])**2 + (0.5 - x[1])**2 + (0.5 - x[2])**2
g = exp(-1.*c)
f = 6*(1.-c)*g
f = 2*(3.-2*c)*g
V = FiniteElement("CG", cell, degree, dirichlet_expression=g)
u = TrialFunction(V)
......
......@@ -3,7 +3,7 @@ cell = "quadrilateral"
x = SpatialCoordinate(cell)
c = (0.5-x[0])**2 + (0.5-x[1])**2
g = exp(-1.*c)
f = 4*(1.-c)*g
f = 2*(2.-2*c)*g
V = FiniteElement("DG", cell, degree)
......
......@@ -3,7 +3,7 @@ cell = "hexahedron"
x = SpatialCoordinate(cell)
c = (0.5-x[0])**2 + (0.5-x[1])**2 + (0.5-x[2])**2
g = exp(-1.*c)
f = 6*(1.-c)*g
f = 2*(3.-2*c)*g
V = FiniteElement("DG", cell, degree)
......
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