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
begin:vcard
fn:Ronan Perrussel
n:Perrussel;Ronan
org;quoted-printable:Laboratoire Amp=C3=A8re
adr:;;36, avenue Guy de Collongue;Ecully;;69134;France
email;internet:[EMAIL PROTECTED]
title;quoted-printable:Charg=C3=A9 de recherche CNRS / Assistant scientist
tel;work:+33 4 72 18 61 09
x-mozilla-html:FALSE
url:http://perso.ec-lyon.fr/ronan.perrussel/
version:2.1
end:vcard
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
