[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Getfem-users] Use of Nedelec finite elements.
From: |
Ronan Perrussel |
Subject: |
[Getfem-users] Use of Nedelec finite elements. |
Date: |
Tue, 30 Oct 2007 19:45:56 +0100 |
User-agent: |
Thunderbird 2.0.0.6 (X11/20071022) |
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
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Getfem-users] Use of Nedelec finite elements.,
Ronan Perrussel <=