Memory access issue #75
Description
@navjotk @pvelesko @vincepandolfo @mlange05
I have an issue with my new Branch test_gradient. I am running the simple adjoint test test_adjointA.py that worked perfectly. I however did some rewrite of Acoustic_codegen and fwi_operators. I would be gratefull if anyone could help with it.
Currently if I run the exact same test 10 times, about 8 out of 10 times I get NaN instead of the proper values, while the other 2 other times the test pass perfectly.
valgrind --tool=memcheck py.test -vs tests/test_adjointA.py &> logAdj
I ran valgrind on it and I get the following error :
==8642== Invalid read of size 16 ==8642== at 0x103E8E30: ForwardOperator (ae8148534a2554d9701a4f62b11b7b82a54c9a34.cpp:38) ==8642== by 0xBF25ADB: ffi_call_unix64 (in /usr/lib/x86_64-linux-gnu/libffi.so.6.0.1) ==8642== by 0xBF2540B: ffi_call (in /usr/lib/x86_64-linux-gnu/libffi.so.6.0.1) ==8642== by 0xBD135FD: _ctypes_callproc (in /usr/lib/python2.7/lib-dynload/_ctypes.x86_64-linux-gnu.so) ==8642== by 0xBD14F9D: ??? (in /usr/lib/python2.7/lib-dynload/_ctypes.x86_64-linux-gnu.so) ==8642== by 0x505F95: PyObject_Call (in /usr/bin/python2.7) ==8642== by 0x49B079: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x4A090B: PyEval_EvalCodeEx (in /usr/bin/python2.7) ==8642== by 0x499A51: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x4A090B: PyEval_EvalCodeEx (in /usr/bin/python2.7) ==8642== by 0x499A51: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x4A1C99: ??? (in /usr/bin/python2.7) ==8642== Address 0xdc926d8 is 33,656 bytes inside a block of size 33,664 alloc'd ==8642== at 0x4C2AB80: malloc (vg_replace_malloc.c:296) ==8642== by 0x8C8AC0B: ??? (in /usr/lib/python2.7/dist-packages/numpy/core/multiarray.so) ==8642== by 0x8C29321: ??? (in /usr/lib/python2.7/dist-packages/numpy/core/multiarray.so) ==8642== by 0x8C2BCD1: ??? (in /usr/lib/python2.7/dist-packages/numpy/core/multiarray.so) ==8642== by 0x8C8C647: ??? (in /usr/lib/python2.7/dist-packages/numpy/core/multiarray.so) ==8642== by 0x49EC75: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x4A090B: PyEval_EvalCodeEx (in /usr/bin/python2.7) ==8642== by 0x499A51: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x499EF1: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x4A1C99: ??? (in /usr/bin/python2.7) ==8642== by 0x505F95: PyObject_Call (in /usr/bin/python2.7) ==8642== by 0x536094: ??? (in /usr/bin/python2.7) ==8642== ==8642== Invalid read of size 16 ==8642== at 0x105EBDD8: AdjointOperator (b79c4970ce58dd574bc69b44d71769d497267caf.cpp:38) ==8642== by 0xBF25ADB: ffi_call_unix64 (in /usr/lib/x86_64-linux-gnu/libffi.so.6.0.1) ==8642== by 0xBF2540B: ffi_call (in /usr/lib/x86_64-linux-gnu/libffi.so.6.0.1) ==8642== by 0xBD135FD: _ctypes_callproc (in /usr/lib/python2.7/lib-dynload/_ctypes.x86_64-linux-gnu.so) ==8642== by 0xBD14F9D: ??? (in /usr/lib/python2.7/lib-dynload/_ctypes.x86_64-linux-gnu.so) ==8642== by 0x505F95: PyObject_Call (in /usr/bin/python2.7) ==8642== by 0x49B079: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x4A090B: PyEval_EvalCodeEx (in /usr/bin/python2.7) ==8642== by 0x499A51: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x4A090B: PyEval_EvalCodeEx (in /usr/bin/python2.7) ==8642== by 0x49AB44: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x4A1C99: ??? (in /usr/bin/python2.7) ==8642== Address 0xdc926d8 is 33,656 bytes inside a block of size 33,664 alloc'd ==8642== at 0x4C2AB80: malloc (vg_replace_malloc.c:296) ==8642== by 0x8C8AC0B: ??? (in /usr/lib/python2.7/dist-packages/numpy/core/multiarray.so) ==8642== by 0x8C29321: ??? (in /usr/lib/python2.7/dist-packages/numpy/core/multiarray.so) ==8642== by 0x8C2BCD1: ??? (in /usr/lib/python2.7/dist-packages/numpy/core/multiarray.so) ==8642== by 0x8C8C647: ??? (in /usr/lib/python2.7/dist-packages/numpy/core/multiarray.so) ==8642== by 0x49EC75: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x4A090B: PyEval_EvalCodeEx (in /usr/bin/python2.7) ==8642== by 0x499A51: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x499EF1: PyEval_EvalFrameEx (in /usr/bin/python2.7) ==8642== by 0x4A1C99: ??? (in /usr/bin/python2.7) ==8642== by 0x505F95: PyObject_Call (in /usr/bin/python2.7) ==8642== by 0x536094: ??? (in /usr/bin/python2.7)