dolfin team mailing list archive

[Marie Rognes <meg@xxxxxxxxxxx>]

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