[Getfem-users] memory problem with gf_asm()

[Hicham Tber]

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
%%