|
From: | Yves Renard |
Subject: | Re: [Getfem-users] Assembing finite element matrices in GetFem++ |
Date: | Wed, 09 Oct 2013 10:25:29 +0200 |
User-agent: | Mozilla/5.0 (X11; Linux x86_64; rv:17.0) Gecko/20130330 Thunderbird/17.0.5 |
Dear Marco, if A(x), B(x) and C(x) are known matrix field described on a finite element method, you can use the generic assembly to compute your mass matrix, yes. It should be similar to the mass matrix assembly in getfem_assembling.h adding the matrix fields. This should give something like generic_assembly assem; assem.set("A=data(qdim(#1),qdim(#1),#2); B=data(qdim(#1),qdim(#1),#2); C=data(qdim(#1),qdim(#1),#2); M(#1,#1)+=comp(Base(#2).Base(#2).Base(#2).vBase(#1).vBase(#1))(k,l,m,:,i,:,j).A(j,o,k).B(o,p,l).C(p,q,m);"); assem.push_mi(mim); assem.push_mf(mf_u1); assem.push_mf(mf_data); assem.push_mat(M); assem.assembly(rg); assuming that mf_data is a scalar fem. Yves. Le 07/10/2013 16:37, Marco Pischedda a écrit :
-- Yves Renard (address@hidden) tel : (33) 04.72.43.87.08 Pole de Mathematiques, INSA-Lyon fax : (33) 04.72.43.85.29 20, rue Albert Einstein 69621 Villeurbanne Cedex, FRANCE http://math.univ-lyon1.fr/~renard --------- |
[Prev in Thread] | Current Thread | [Next in Thread] |