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)


Reply via email to