ffc team mailing list archive

[Marie Rognes <meg@xxxxxxxxxxx>]

Bartosz Sawicki wrote:
> 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.
>
>   

The reason for why are there are no examples using the Nedelec elements 
is probably that we have not finished implementing and testing them yet. 
In other words, they are not working properly yet.


> My first, the simplest try was:
>
> element = VectorElement("Nedelec", "tetrahedron", 1, 3)
>   

Further, the Nedelec elements are vector-elements by  default, so you want

   element = FiniteElement("Nedelec", "tetrahedron", 1)

here. (With the notation you used above, you get a vector of 3 Nedelec 
vector elements.)

> v = TestFunction(element)
> t = TrialFunction(element)
> c = VectorConstant("tetrahedron")
> a = dot(curl(v),curl(t))*dx
> L = c*dx
>   

In addition L should be a linear functional. Thus c*dx where c is a 
vector does not make sense in the ffc-regime. Maybe you want something like

   L = dot(curl(v), c)*dx

instead?


--
Marie E. Rognes
Ph.D Fellow, 
Centre of Mathematics for Applications, 
University of Oslo
http://folk.uio.no/meg