Hi,

I encountered a rather severe memory leak when using pygsl_rng.dirichlet_lnpdf. I'm using gsl 1.8 via pygsl 0.3.2.

The following code snippet reproduces the error:

import pygsl.rng
for i in range(800000):
    pygsl.rng.dirichlet_lnpdf([1.1,1.1,1.1,1.1],[[0.25,0.25,0.25,0.25]])

When running this you will see the memory consumption increase monotonously.

The problem did not occur when I implemented the density function myself using pygsl.sf.gamma, nor for a C test script that called gsl_ran_dirichlet_pdf directly. This seems to indicate that the problem resides somewhere in the pygsl interface to gsl_ran_dirichlet_pdf.

Best,
Benjamin



-------------------------------------------------------
Using Tomcat but need to do more? Need to support web services, security?
Get stuff done quickly with pre-integrated technology to make your job easier
Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
_______________________________________________
pygsl-discuss mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/pygsl-discuss

Reply via email to