diff --git a/test/sumfact/poisson/poisson_dg_3d.ufl b/test/sumfact/poisson/poisson_dg_3d.ufl index b7e20157ef4f7015b856b572d0a1d84da1800f1e..d5c0318228311f2823c36cd16daa62be35cb106d 100644 --- a/test/sumfact/poisson/poisson_dg_3d.ufl +++ b/test/sumfact/poisson/poisson_dg_3d.ufl @@ -1,9 +1,12 @@ 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 +#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) V = FiniteElement("DG", cell, degree)