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
