Dear Iago, you should replace : "V(#1)+=comp(vBase(#2)(i,j).Normal()(j).Base(#1)(:)).d(i);" by : "V(#1)+=comp(vBase(#2).Normal().Base(#1))(i,j,j,:).d(i);" and it should work fine.
Best regards, Ronan Iago Barbeiro a écrit :
template<typename VECT1, typename VECT2>void asm_Fbc_source_term(VECT1 &F, const mesh_im &mim, const mesh_fem &mf,const mesh_fem &mf_data, const VECT2 &D, const mesh_region &rg) { GMM_ASSERT1(mf_data.get_qdim() == 2, "invalid data mesh_fem"); generic_assembly assem;assem.set("d=data(#2);" "V(#1)+=comp(vBase(#2)(i,j).Normal()(j).Base(#1)(:)).d(i);");assem.push_mi(mim); assem.push_mf(mf); assem.push_mf(mf_data); assem.push_data(D); assem.push_vec(F); assem.assembly(rg);}
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
