[Getfem-users] Use of Nedelec finite elements.

[Ronan Perrussel]

Dear getfem++ users,

I would like to use Nedelec finite elements for some computations. I 
made some small two dimensional experiments for evaluating the use of 
these finite element family but I do not manage to make all what I wish.
For assembling a term like "\int_\Omega \curl(E) \curl(E') dx", it seems 
to be fine with the following piece of code :
""
   getfem::generic_assembly assem;
   assem.push_mi(mi);
   assem.push_mf(mf);
   assem.push_mat(K);
   assem.set("M(#1,#1)+=comp(vGrad(#1).vGrad(#1))(:,2,1,:,2,1)"
         "+comp(vGrad(#1).vGrad(#1))(:,1,2,:,1,2)"
         "-comp(vGrad(#1).vGrad(#1))(:,2,1,:,1,2)"
         "-comp(vGrad(#1).vGrad(#1))(:,1,2,:,2,1)");
   assem.assembly();
""
but for assembling a term like "\int_\Omega E . E' dx", the simple piece 
of code for the mass matrix does not give the good answer :
""
  getfem::asm_mass_matrix(M, mi, mf);
""
and I do not understand exactly why.

Does somebody have a suggestion for assembling the term " \int_\Omega E 
. E' dx"?

Thank you in advance for the answer,
Best regards,
Ronan Perrussel

ronan_perrussel.vcf
Description: Vcard