Error:
Traceback (most recent call last):
File "test.py", line 33, in ?
u0.assign(w)
File
"/home/meg/src/dolfin/dolfin/local/lib/python2.4/site-packages/dolfin/dolfin.py",
line 3172, in assign
return _dolfin.cpp_Function_assign(*args)
RuntimeError: *** Error: Unable to find finite element in library:
"Mixed finite element: [Crouzeix-Raviart finite element of degree 1
on a triangle, Crouzeix-Raviart finite element of degree 1 on a
triangle]".
Segmentation fault
Code:
from dolfin import *
# Geometry of the domain and mesh
shape = "triangle"
mesh = UnitSquare(16, 16)
# Setting up the elements
V = VectorElement("CR", shape, 1)
u0 = Function(V, mesh, 0.0)
v = TestFunction(V)
u = TrialFunction(V)
a = dot(grad(u), grad(v))*dx
L = dot(u0, v)*dx
pde = LinearPDE(a, L, mesh)
w = pde.solve()
u0.assign(w)
Easy fix?
--
Marie E. Rognes
Ph.D Fellow,
Centre of Mathematics for Applications,
University of Oslo
http://folk.uio.no/meg
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
