diff --git a/test/sumfact/poisson/poisson_dg_3d.ufl b/test/sumfact/poisson/poisson_dg_3d.ufl index d5c0318228311f2823c36cd16daa62be35cb106d..b7e20157ef4f7015b856b572d0a1d84da1800f1e 100644 --- a/test/sumfact/poisson/poisson_dg_3d.ufl +++ b/test/sumfact/poisson/poisson_dg_3d.ufl @@ -1,12 +1,9 @@ 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)