dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #17575
Re: [Bug 521507] [NEW] Memory leak in Python Expression interface
Johan Hake wrote:
> On Saturday 13 February 2010 13:08:05 Garth Wells wrote:
>> Public bug reported:
>>
>> The below code leads to a memory leak (I see a leak over about 2MB over
>> 90s). The leak is related to the Expression. If I comment out the call
>> to
>>
>> eval(Array<double>&, const Array<double>&);
>>
>> from
>>
>> dolfin::GenericFunction::evaluate(double*, const double*, const
>> ufc::cell) const
>>
>> I don't detect a leak.
>
> I can confirm this bug. Not sure it is in the typemap, as it should also have
> shown up if we moved f = MyExpression(degree=1) outside of the loop.
>
I've add some scripts for testing memory to
lp:~dolfin-core/dolfin/sandbox
under 'memtest'.
I'm seeing a case where a call-back from C++ to a Python sublass of
NonlinearProblem is causing a big leak.
Garth
> Johan
>
--
Memory leak in Python Expression interface
https://bugs.launchpad.net/bugs/521507
You received this bug notification because you are a member of DOLFIN
Team, which is subscribed to DOLFIN.
Status in DOLFIN: Confirmed
Bug description:
The below code leads to a memory leak (I see a leak over about 2MB over 90s). The leak is related to the Expression. If I comment out the call to
eval(Array<double>&, const Array<double>&);
from
dolfin::GenericFunction::evaluate(double*, const double*, const ufc::cell) const
I don't detect a leak.
------
from dolfin import *
parameters["form_compiler"]["cpp_optimize"] = True
parameters["form_compiler"]["optimize"] = True
class MyExpression(Expression):
def eval(self, values, x):
if x[0] < DOLFIN_EPS:
values[0] = 1.0
mesh = UnitSquare(32, 32)
for t in xrange(20000):
for level in xrange(100):
V = FunctionSpace(mesh, "CG", 1)
v = TestFunction(V)
u = Function(V)
u0 = Function(V)
f = MyExpression(degree=1)
L = v*(u-u0)*dx + v*f*ds
b = assemble(L)
print "----"
References