[Top][All Lists]

[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 (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 :


Attachment: ronan_perrussel.vcf
Description: Vcard

reply via email to

[Prev in Thread] Current Thread [Next in Thread]