Thank you for the quick reply. I created a bug report (Issue #168). Cheers, Torsten
> 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 > -- Torsten Wendav Graduate Student | Theoretical Optics & Photonics Humboldt Universitaet zu Berlin [email protected] +49-30-2093 7618 _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
