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

Switch the poisson_dg_3d input to symbolic

parent 6ac8e4ff
No related branches found
No related tags found
No related merge requests found
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 = Expression("Dune::FieldVector<double,2> c(0.5); c-= x; return 4.*(1.-c.two_norm2())*std::exp(-1.*c.two_norm2());", cell=cell)
g = Expression("Dune::FieldVector<double,2> c(0.5); c-= x; return std::exp(-1.*c.two_norm2());", on_intersection=True, cell=cell)
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
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