From hg repository. pull&update 15 minutes ago.
BArtek
Garth N. Wells wrote:
> Which version of DOLFIN are you using?
>
> Garth
>
> Bartosz Sawicki wrote:
>> I realized strange problem in my code when Function is created using
>> FiniteElement and DofMap. PETSc have reported memory problems, but it
>> is difficult to repeat. I'm not able to locate the error, but it
>> seems that the bug is inside the FunctionSpace.
>>
>> Please consider following code. It compiles, it works works, but
>> Valgrind reports errors.
>>
>> ================
>>
>> #include <dolfin.h>
>>
>> using namespace dolfin;
>>
>> int main (int argc, char **argv) {
>>
>> UnitCube mesh(10,10,10);
>>
>> FiniteElement fe("FiniteElement('Lagrange', 'tetrahedron', 1)");
>> DofMap dofmap("FFC dof map for FiniteElement('Lagrange',
>> 'tetrahedron', 1)", mesh);
>>
>> Function j( FunctionSpace(mesh, fe, dofmap) );
>>
>> j.vector() = 1.0;
>>
>>
>> Point point( 0.5, 0.5, 0.5);
>> Point value;
>>
>> j.eval((double *)value.coordinates(), point.coordinates());
>>
>> std::cout << value[0] << std::endl;
>>
>> return 0;
>> }
>> ================
>>
>> $ ./demo
>> 1
>>
>> $ valgrind ./demo
>> ==20870== Memcheck, a memory error detector.
>> ==20870== Copyright (C) 2002-2007, and GNU GPL'd, by Julian Seward et al.
>> ==20870== Using LibVEX rev 1854, a library for dynamic binary
>> translation.
>> ==20870== Copyright (C) 2004-2007, and GNU GPL'd, by OpenWorks LLP.
>> ==20870== Using valgrind-3.3.1-Debian, a dynamic binary
>> instrumentation framework.
>> ==20870== Copyright (C) 2000-2007, and GNU GPL'd, by Julian Seward et al.
>> ==20870== For more details, rerun with: -v
>> ==20870==
>> ==20870== Invalid write of size 8
>> ==20870== at 0x4825B8A: VecGetValues_Seq (in
>> /usr/lib/libpetscvec.so.2.3.3)
>> ==20870== by 0x4851887: VecGetValues (in
>> /usr/lib/libpetscvec.so.2.3.3)
>> ==20870== by 0x41EAD22: dolfin::PETScVector::get(double*, unsigned,
>> unsigned const*) const (PETScVector.cpp:144)
>> ==20870== by 0x4188024: dolfin::Function::interpolate(double*,
>> dolfin::FunctionSpace const&, ufc::cell const&, unsigned, int) const
>> (Function.cpp:266)
>> ==20870== by 0x418041B: dolfin::FunctionSpace::eval(double*, double
>> const*, dolfin::Function const&, ufc::cell const&, unsigned) const
>> (FunctionSpace.cpp:149)
>> ==20870== by 0x41842D6: dolfin::FunctionSpace::eval(double*, double
>> const*, dolfin::Function const&) const (FunctionSpace.cpp:132)
>> ==20870== by 0x41886DA: dolfin::Function::eval(double*, double
>> const*) const (Function.cpp:206)
>> ==20870== by 0x8054EEF: main (main.cpp:38)
>> ==20870== Address 0x64f71d8 is 0 bytes inside a block of size 32 free'd
>> ==20870== at 0x40241AA: operator delete[](void*)
>> (vg_replace_malloc.c:364)
>> ==20870== by 0x4180107: dolfin::FunctionSpace::Scratch::~Scratch()
>> (FunctionSpace.cpp:302)
>> ==20870== by 0x4180A98: dolfin::FunctionSpace::~FunctionSpace()
>> (FunctionSpace.cpp:64)
>> ==20870== by 0x8054E82: main (main.cpp:22)
>> ==20870==
>> ==20870== Invalid write of size 8
>> ==20870== at 0x40ED491:
>> ffc_04_finite_element_0::evaluate_basis(unsigned, double*, double
>> const*, ufc::cell const&) const (ffc_04.h:173)
>> ==20870== by 0x4180496: dolfin::FunctionSpace::eval(double*, double
>> const*, dolfin::Function const&, ufc::cell const&, unsigned) const
>> (FiniteElement.h:60)
>> ==20870== by 0x41842D6: dolfin::FunctionSpace::eval(double*, double
>> const*, dolfin::Function const&) const (FunctionSpace.cpp:132)
>> ==20870== by 0x41886DA: dolfin::Function::eval(double*, double
>> const*) const (Function.cpp:206)
>> ==20870== by 0x8054EEF: main (main.cpp:38)
>> ==20870== Address 0x64f7228 is 0 bytes inside a block of size 8 free'd
>> ==20870== at 0x40241AA: operator delete[](void*)
>> (vg_replace_malloc.c:364)
>> ==20870== by 0x4180116: dolfin::FunctionSpace::Scratch::~Scratch()
>> (FunctionSpace.cpp:303)
>> ==20870== by 0x4180A98: dolfin::FunctionSpace::~FunctionSpace()
>> (FunctionSpace.cpp:64)
>> ==20870== by 0x8054E82: main (main.cpp:22)
>> ==20870==
>> ==20870== Invalid read of size 8
>> ==20870== at 0x41804B8: dolfin::FunctionSpace::eval(double*, double
>> const*, dolfin::Function const&, ufc::cell const&, unsigned) const
>> (FunctionSpace.cpp:158)
>> ==20870== by 0x41842D6: dolfin::FunctionSpace::eval(double*, double
>> const*, dolfin::Function const&) const (FunctionSpace.cpp:132)
>> ==20870== by 0x41886DA: dolfin::Function::eval(double*, double
>> const*) const (Function.cpp:206)
>> ==20870== by 0x8054EEF: main (main.cpp:38)
>> ==20870== Address 0x64f71d8 is 0 bytes inside a block of size 32 free'd
>> ==20870== at 0x40241AA: operator delete[](void*)
>> (vg_replace_malloc.c:364)
>> ==20870== by 0x4180107: dolfin::FunctionSpace::Scratch::~Scratch()
>> (FunctionSpace.cpp:302)
>> ==20870== by 0x4180A98: dolfin::FunctionSpace::~FunctionSpace()
>> (FunctionSpace.cpp:64)
>> ==20870== by 0x8054E82: main (main.cpp:22)
>> ==20870==
>> ==20870== Invalid read of size 8
>> ==20870== at 0x41804BB: dolfin::FunctionSpace::eval(double*, double
>> const*, dolfin::Function const&, ufc::cell const&, unsigned) const
>> (FunctionSpace.cpp:158)
>> ==20870== by 0x41842D6: dolfin::FunctionSpace::eval(double*, double
>> const*, dolfin::Function const&) const (FunctionSpace.cpp:132)
>> ==20870== by 0x41886DA: dolfin::Function::eval(double*, double
>> const*) const (Function.cpp:206)
>> ==20870== by 0x8054EEF: main (main.cpp:38)
>> ==20870== Address 0x64f7228 is 0 bytes inside a block of size 8 free'd
>> ==20870== at 0x40241AA: operator delete[](void*)
>> (vg_replace_malloc.c:364)
>> ==20870== by 0x4180116: dolfin::FunctionSpace::Scratch::~Scratch()
>> (FunctionSpace.cpp:303)
>> ==20870== by 0x4180A98: dolfin::FunctionSpace::~FunctionSpace()
>> (FunctionSpace.cpp:64)
>> ==20870== by 0x8054E82: main (main.cpp:22)
>> 1
>> ==20870==
>> ==20870== ERROR SUMMARY: 16 errors from 4 contexts (suppressed: 353
>> from 2)
>> ==20870== malloc/free: in use at exit: 16 bytes in 2 blocks.
>> ==20870== malloc/free: 463 allocs, 461 frees, 168,090 bytes allocated.
>> ==20870== For counts of detected errors, rerun with: -v
>> ==20870== searching for pointers to 2 not-freed blocks.
>> ==20870== checked 4,507,820 bytes.
>> ==20870==
>> ==20870== LEAK SUMMARY:
>> ==20870== definitely lost: 16 bytes in 2 blocks.
>> ==20870== possibly lost: 0 bytes in 0 blocks.
>> ==20870== still reachable: 0 bytes in 0 blocks.
>> ==20870== suppressed: 0 bytes in 0 blocks.
>> ==20870== Rerun with --leak-check=full to see details of leaked memory.
>>
>>
>> BArtek
>> _______________________________________________
>> DOLFIN-dev mailing list
>> [email protected]
>> http://www.fenics.org/mailman/listinfo/dolfin-dev
>
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
