[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] memory problem with gf_asm()
From: |
Hicham Tber |
Subject: |
Re: [Getfem-users] memory problem with gf_asm() |
Date: |
Wed, 08 Oct 2008 14:08:29 +0200 |
User-agent: |
Thunderbird 2.0.0.6 (X11/20070801) |
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
>> address@hidden
>> https://mail.gna.org/listinfo/getfem-users
>>
>>
>