Re: [Getfem-users] memory problem with gf_asm()

[Ronan Perrussel]

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
> address@hidden
> https://mail.gna.org/listinfo/getfem-users
>
>

ronan_perrussel.vcf
Description: Vcard