Trust Region solver memory corruption

Michael_Sarahan — Thu, 05 Dec 2013 01:37:57 GMT

Hi there,

I'm trying to use MKL to replace existing Levenberg-Marquardt fitting using MINPACK or MPFit. As a test case, I'm fitting a Gaussian function. Things work pretty well, except that my program crashes at the final MKL_FreeBuffers call. I don't know if this is an MKL problem, as this crash occurs even when I comment out that MKL_FreeBuffers call.

I've pasted part of my code here: http://pastebin.com/SXuQ8egA - if more of the code would be helpful, I'll upload it. It's just a rather large codebase to ask you to crawl through.

I had a few questions, just to make sure I'm not doing anything dumb:

Is it OK to use vectors as I've done, instead of malloc? It seems to work (the algorithm fits my test Gaussian correctly), but memory corruption occurs at the end occurs in a vector's destructor. (Do you have a good way to identify which vector is causing the problem? I'm afraid I don't know how to debug that.)
Does MKL keep any handle on any of the temporary arrays that would cause that destructor problem?

For reference, I'm using VS2008, x64 and MKL 10.3.1.

Thanks for any tips...

PS: I also sometimes get

Michael_Sarahan — Thu, 05 Dec 2013 01:41:58 GMT

PS: I also sometimes get access violation errors in my djacobix call, but that's not every time.

Sorry, my error. I was not

Michael_Sarahan — Fri, 06 Dec 2013 00:56:07 GMT

Sorry, my error. I was not creating the fjac vector correctly - it was too small. I had only m elements.

Sorry for the noise!

topic Sorry, my error. I was not in Intel® oneAPI Math Kernel Library

Trust Region solver memory corruption

PS: I also sometimes get

Sorry, my error. I was not