ffc team mailing list archive
-
ffc team
-
Mailing list archive
-
Message #01365
Re: Nedelec elements
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).
/Anders
Follow ups
References