Hi Ronan I've got the same problem although using an external subroutine calling push-pop.
Hicham > 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 >> >> > _______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
