Yes, this is a bug. Please report it as an issue and I will take a look at it. Looks like a rounding error issue.
Use the following simpler test program in the issue report: from dolfin import * mesh = UnitCubeMesh(14,14,14) V = FunctionSpace(mesh, "Lagrange", 1) v = Function(V) x = (0.75, 0.25, 0.125) print v(x) -- Anders On Thu, Dec 05, 2013 at 10:13:14AM +0100, Mikael Mortensen wrote: > Hi, > > I get the same error with development version of dolfin so this is likely a > bug. You can create an issue at https://bitbucket.org/fenics-project/dolfin. > > Mikael > > Den Dec 4, 2013 kl. 5:46 PM skrev Torsten Wendav: > > > from dolfin import * > import numpy > > # Mesh > mesh = UnitCubeMesh(14,14,14) > V = FunctionSpace(mesh, "Lagrange", 1) > > # BC > u0 = Expression("1 + x[0]*x[0] + 2*x[1]*x[1] + x[2]*x[2]") > bc = DirichletBC(V,u0,DomainBoundary()) > > # Define variational problem > u, v = TrialFunction(V), TestFunction(V) > f = Constant(-6.0) > a = inner(nabla_grad(u), nabla_grad(v))*dx > L = f*v*dx > > # Compute solution > u = Function(V) > set_log_level(PROGRESS) > solve(a==L,u,bc) > > x = numpy.array([0.75, 0.25, 0.125]) > print u(x) > > > _______________________________________________ > fenics mailing list > [email protected] > http://fenicsproject.org/mailman/listinfo/fenics _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
