Re: [Getfem-users] boundary source term

[Ronan Perrussel]

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);
>       
>     }

ronan_perrussel.vcf
Description: Vcard