← Back to team overview

ffc team mailing list archive

Re: Nedelec elements

 

Anders Logg wrote:
On Fri, Nov 02, 2007 at 11:53:33PM +0100, Bartosz Sawicki wrote:
   element = FiniteElement("Nedelec", "tetrahedron", 1)
Yeeeap :), it was the problem. Now fcc compiles correctly.
It might compile, but I wouldn't count on getting anything sensible
(but you might be lucky). As Marie points out, the Nedelec elements
need some more work.
Usually I am lucky, but not this time.
ffc compiles, gcc compiles, but during the runtime:

   | Assembling over cells                                           |
   |-----------------------------------------------------------------| 0.0%
terminate called after throwing an instance of 'std::runtime_error'
   what():  evaluate_dof not implemented for this type of element
Aborted (core dumped)

... ups.
I'll be waiting for any news about Nedelec elements development and I'll be happy, if I can help in testing.

regards,
BArtosz

The first thing to try would be to *not* use a Nedelec element for the
right-hand side (the source term). Pick a standard Lagrange element,
say DG(1) for the right-hand. That should fix this problem, but I
still think something else will break (or simply give the wrong
values).


An error of this form would probably also occur, if you try to apply strong/essential boundary conditions for the Nedelec element space. Figuring out how to do this might take some time.

The "rest" should be working soon.

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



Follow ups

References