Sébastien Janas a écrit :
Is there any systematic procedure to impose the dirichlet conditions
on my system A*X=B ?
I have tried with the penalty method, but the sparse matrix format
included in getfem set it very difficult (we can't make B[i] += A[i,i]
* BIG * val because A[i,i] give an 1x1 array and not a scalar...)
Dear Sébastien,
Did you try to use the instruction :
[spmat H,vec R] = gf_asm('dirichlet', int boundary_num, mesh_im mim,
mesh_fem mf_u, mesh_fem mf_d, mat Hd, vec Rd)?
explained here :
http://download.gna.org/getfem/doc/getfem_matlab/gfm_47.html
Ronan
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
