[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] Transient analysis-dirichlet conditions
From: |
Ronan Perrussel |
Subject: |
Re: [Getfem-users] Transient analysis-dirichlet conditions |
Date: |
Thu, 05 Jun 2008 08:32:53 +0200 |
User-agent: |
Thunderbird 2.0.0.14 (X11/20080505) |
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