Hi Hicham,

maybe you should follow the example on the documentation page concerning gf_workspace :
http://download.gna.org/getfem/doc/getfem_matlab/gfm_18.html

Nevertheless, even with these instructions I have also some problems for managing the memory used by getfem in matlab : I call an external subroutine, which follows the advices of the documentation concerning gf_workspace, but after several calls of the subroutine matlab finally crashes.
Does someone experiment the same behavior?

Best regards,
Ronan

Hicham Tber a écrit :
Hello,
I am using gf_asm(), from the Matlab interface, inside a loop. I realized that 
the used memory increases with iterations until my program runs out of memory. 
I would like to know if this is a bug or just something wrong with my program. 
The example below illustrate my problem.

Thank you in advance.

Hicham


%
% dx = 0.01;
Omega = gf_mesh('triangles grid',[0:dx:1],[0:dx:1]);
Vh = gf_mesh_fem(Omega,1);
gf_mesh_fem_set(Vh,'fem',gf_fem('FEM_PK(2,1)'));
im = gf_mesh_im(Omega, gf_integ('IM_EXACT_SIMPLEX(2)')) ;
VhNdof=gf_mesh_fem_get(Vh, 'nbdof');

e = ones(VhNdof,1); Nmax = 100;

for k = 1:Nmax
     Mat = gf_asm('laplacian', im, Vh, Vh,  ones(1,VhNdof));
% memory monitoring on Linux !free -m
end
%%





_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users


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

Reply via email to