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).