My version of the adaptivity prototype code from the 2010 Rognes and Logg paper
"Automated goal-oriented error control I" results in a segmentation fault error:
Generating forms required for error control, this may take some time...
Solving variational problem adaptively
Adaptive iteration 0
Solving linear variational problem.
Segmentation fault: 11
Python quit unexpectedly while using the libdolfin.1.2.dylib plug-in.
Is that a coding issue or is my installation not working properly?
I am running fenics 1.2.0 on OS X 10.8.5 with Xcode 5.0.2
Below is the code.
Best, Achim
from dolfin import *
mesh = UnitSquareMesh(4, 4)
V = FunctionSpace(mesh , "CG", 1)
u = TrialFunction(V)
v = TestFunction(V)
u_h = Function(V)
a = dot(grad (u), grad (v))*dx
L = v*dx
def boundary (x):
return x[0] < 1E-14
bc = DirichletBC(V, Constant (0.0), boundary)
M = u*dx
solve(a == L, u_h, bcs=bc, tol=1e-3, M=M)
_______________________________________________
fenics-support mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics-support
