dolfin team mailing list archive

[Ted Kord <question98201@xxxxxxxxxxxxxxxxxxxxx>]

New question #98201 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/98201

Hi

I get the following error withh ffc when trying to generate the C++ header for a Q1 element.

In the ufl file, I have: 

element = FiniteElement("Lagrange", "quadrilateral", 1)
v      = TestFunction(element)
u      = TrialFunction(element)

a = inner(grad(u), grad(v))*dx
L = v*f*dx

The output of ffc -l dolfin ?.ufl is:

Compiler stage 1: Analyzing form
--------------------------------
  Name:                               a
  Rank:                               2
  Cell:                               <quadrilateral of degree 1>
  Topological dimension:              2
  Geometric dimension:                2
  Number of facets:                   4
  Number of coefficients:             0
  Number of cell integrals:           1
  Number of exterior facet integrals: 0
  Number of interior facet integrals: 0
  Number of macro cell integrals:     0
  Number of cell domains:             1
  Number of exterior facet domains:   0
  Number or interior facet domains:   0
  Number of macro cell domains:       0
  Arguments:                          [v_0, v_1]
  Coefficients:                       []
  Argument names:                     [v, u]
  Coefficient names:                  []
  Unique elements:                    CG1
  Unique sub elements:                CG1
  
  Quadrature rule: default
  Quadrature order: auto
  Automatic selection of representation not implemented, defaulting to quadrature.
  Integral quadrature degree is 0.
Traceback (most recent call last):
  File "/home/tedkord/FEniCS/build/bin/ffc", line 167, in <module>
    sys.exit(main(sys.argv[1:]))
  File "/home/tedkord/FEniCS/build/bin/ffc", line 155, in main
    compile_form(ufd.forms, ufd.object_names, prefix, options)
  File "/home/tedkord/FEniCS/build/lib/python2.6/site-packages/ffc/compiler.py", line 87, in compile_form
    form, form_data = analyze_form(form, object_names, options)
  File "/home/tedkord/FEniCS/build/lib/python2.6/site-packages/ffc/compiler.py", line 155, in analyze_form
    form_data.ffc_elements = [create_element(element) for element in form_data.elements]
  File "/home/tedkord/FEniCS/build/lib/python2.6/site-packages/ffc/createelement.py", line 52, in create_element
    ffc_element = FFCFiniteElement(ufl_element, domain)
  File "/home/tedkord/FEniCS/build/lib/python2.6/site-packages/ffc/finiteelement.py", line 82, in __init__
    ufl_element.cell().domain(), ufl_element.degree())
  File "/home/tedkord/FEniCS/build/lib/python2.6/site-packages/ffc/finiteelement.py", line 237, in __choose_element
    fiat_shape = ufl_domain2fiat_domain[shape]
KeyError: 'quadrilateral'

What's the problem?

Ted



-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.