ffc team mailing list archive

[Bartosz Sawicki <sawickib@xxxxxxxxxxxxx>]

Hello everybody,

I'm trying to use Nedelec element for curl-curl equation. I can't find
any example of using Nedelec edge elements in ffc and dolfin, so I'm not
sure how to use it.

My first, the simplest try was:

element = VectorElement("Nedelec", "tetrahedron", 1, 3)
v = TestFunction(element)
t = TrialFunction(element)
c = VectorConstant("tetrahedron")
a = dot(curl(v),curl(t))*dx
L = c*dx

but compilation fails:

$ ffc -d1 Nedelec.form
This is FFC, the FEniCS Form Compiler, version 0.4.3.
For further information, go to http://www/fenics.org/ffc/.
Preprocessing form file: Nedelec.form --> Nedelec.py
Traceback (most recent call last):
  File "/usr/bin/ffc", line 180, in ?
    sys.exit(main(sys.argv[1:]))
  File "/usr/bin/ffc", line 107, in main
    execfile(outname, ns)
  File "Nedelec.py", line 15, in ?
    a = dot(curl(v),curl(t))*dx
  File
"/usr/lib/python2.4/site-packages/ffc/compiler/language/operators.py",
line 181, in curl
    return rot(v)
  File
"/usr/lib/python2.4/site-packages/ffc/compiler/language/operators.py",
line 175, in rot
    raise FormError, (v, "Rotation only defined for v : R^3 --> R^3")

I really don't understand what I'm doing wrong. Could you give me any
example of proper use of Nedelec elements in ffc?
It would be great if such example will be in dolfin/src/demo directory :)

I'm using the newest FIAT and FFC just from hg repository.

Regards,
BArtosz Sawicki