Re: [Getfem-users] Transient analysis-dirichlet conditions

[Ronan Perrussel]

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

ronan_perrussel.vcf
Description: Vcard