FFC

bubble elements not supported on manifolds

Bug #1154721 reported by Andrew McRae on 2013-03-13
6
This bug affects 1 person
Affects Status Importance Assigned to Milestone
FFC
Medium
Marie Rognes

Bug Description

Using trunk build. Mesh attached.

===========================
from dolfin import *
mesh = Mesh("sphere_ico3.xml")
E = FunctionSpace(mesh, "B", 3)
#E = FunctionSpace(mesh, "CG", 2)
===========================

Calling FFC just-in-time (JIT) compiler, this may take some time.
Traceback (most recent call last):
  File "bubblebug.py", line 3, in <module>
    E = FunctionSpace(mesh, "B", 3)
  File "/home/atm112/local/fenics/src/dolfin/current/local/lib/python2.7/site-packages/dolfin/functions/functionspace.py", line 390, in __init__
    FunctionSpaceBase.__init__(self, mesh, element, constrained_domain)
  File "/home/atm112/local/fenics/src/dolfin/current/local/lib/python2.7/site-packages/dolfin/functions/functionspace.py", line 84, in __init__
    ufc_element, ufc_dofmap = jit(self._ufl_element)
  File "/home/atm112/local/fenics/src/dolfin/current/local/lib/python2.7/site-packages/dolfin/compilemodules/jit.py", line 66, in mpi_jit
    return local_jit(*args, **kwargs)
  File "/home/atm112/local/fenics/src/dolfin/current/local/lib/python2.7/site-packages/dolfin/compilemodules/jit.py", line 154, in jit
    return jit_compile(form, parameters=p, common_cell=common_cell)
  File "/home/atm112/local/fenics/lib/python2.7/site-packages/ffc/jitcompiler.py", line 75, in jit
    return jit_element(ufl_object, parameters)
  File "/home/atm112/local/fenics/lib/python2.7/site-packages/ffc/jitcompiler.py", line 244, in jit_element
    compiled_form, module, form_data, prefix = jit_form(form, parameters)
  File "/home/atm112/local/fenics/lib/python2.7/site-packages/ffc/jitcompiler.py", line 197, in jit_form
    parameters=parameters)
  File "/home/atm112/local/fenics/lib/python2.7/site-packages/ffc/compiler.py", line 156, in compile_form
    ir = compute_ir(analysis, parameters)
  File "/home/atm112/local/fenics/lib/python2.7/site-packages/ffc/representation.py", line 86, in compute_ir
    for (i, e) in enumerate(elements)]
  File "/home/atm112/local/fenics/lib/python2.7/site-packages/ffc/representation.py", line 112, in _compute_element_ir
    element = create_element(ufl_element)
  File "/home/atm112/local/fenics/lib/python2.7/site-packages/ffc/fiatinterface.py", line 100, in create_element
    element = _create_fiat_element(ufl_element)
  File "/home/atm112/local/fenics/lib/python2.7/site-packages/ffc/fiatinterface.py", line 154, in _create_fiat_element
    return RestrictedElement(V, _indices(V, "interior", dim), None)
  File "/home/atm112/local/fenics/lib/python2.7/site-packages/ffc/fiatinterface.py", line 280, in _indices
    return element.entity_dofs()[dim][0]
KeyError: 3

Andrew McRae (andymc) wrote :
Marie Rognes (meg-simula) wrote :

Even easier to reproduce without custom mesh using:

from dolfin import *
mesh = BoundaryMesh(UnitCubeMesh(5, 5, 5), "exterior")
E = FunctionSpace(mesh, "B", 3)

Changed in ffc:
status: New → Confirmed
importance: Undecided → Medium
assignee: nobody → Marie Rognes (meg-simula)
To post a comment you must log in.
This report contains Public information  Edit
Everyone can see this information.

Other bug subscribers

Bug attachments