Re: boundary region

[Edouard Oudet]

Dear all,

I think I solved my problem by a call like

getfem.asm_boundary(bregion_id, expr, mim, mf, mf, a)

where

 expr = 'a = data(#2); M(#1,#2)+= comp(vBase(#1).vBase(#2).Base(#1))(:, 
i,: ,i , j).a(j)'

I am not completely sure to understand the indexing  (:, i,: ,i , 
j).a(j) in previous expression..
Best,
Edouard.




On 6/12/21 8:39 PM, Edouard Oudet wrote:
> Dear Getfem users,
> It seems that the python function `asm_boundary_qu_term` does not 
> accept a (boundary faces) region argument.
> Is there a way to compute the qu matrix only considering a sublist of 
> boundary faces (/ segments) from python?
> Thanks,
> best,
> Edouard.

--
Edouard Oudet : http://www-ljk.imag.fr/membres/Edouard.Oudet/
IMAG - Bureau 164
700 avenue Centrale
38400 Saint Martin d'Hères
+33 (0)4 57 42 17 71  (office LJK)