Hi all,
I have been reading up on the Nedelec element to use it
in an MHD problem, but I am having a little bit of trouble and
was hoping to ask a couple of questions. Specifically, I would
like to use the Nedelec element of first type of order k = 2.
I have dealii 7.0.0 installed. The questions I have are
1) In a scratch c++ file I have asked for the element by including the
file in " #include <fe/fe_nedelec.h> ". Using the command
FE_Nedelec<3> finite_element(1);
it looks like I get the lowest order Nedelec element of the first type
due to the number of degrees of freedom = 12 returned by
finite_element.dofs_per_cell;
But this seems to disagree with the detailed description on the webpage
http://www.dealii.org/developer/doxygen/deal.II/classFE__Nedelec.html#a8e6ff21bacb7ddf2e60075cff693b652
where the convention of Brezzi and Raviart is used and the command
FE_Nedelec<3> finite_element(0);
should give the lowest order element. However, this gives an error.
I was wondering why this might be?
2) I need the Nedelec element of first type and of order k=2, but when I ask
for this element with the command
FE_Nedelec<3> finite_element(2);
I get the error (the same error when 2 is replaced with 0) :
The violated condition was:
degree == 1
The name and call sequence of the exception was:
ExcNotImplemented()
I was wondering if this means that the Nedelec element of order 2
is not implemented? In an e-mail from Dr. Bangerth on deal.II version 7.0.0
it is said that Nedelec elements are now available to arbitrary order,
so I am confused.
Thanks for your time.
Daniel Brauss
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
