Call lapack routine from inside SImulink S function on 64bit

Question

Matthew 2013-7-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/82957-call-lapack-routine-from-inside-simulink-s-function-on-64bit

I am trying to call a LAPACK routine (sgetrf) from inside a Simulink S function on a 64bit machine using libmwlapack. However, I get the error message: MKL ERROR: Parameter 1 was incorrect on entry to SGETRF

Here is the function prototype: extern void sgetrf_(mwSignedIndex *, mwSignedIndex *, float (*)[3], mwSignedIndex *, mwSignedIndex *, mwSignedIndex *);

And the function call: sgetrf_(&nvars,&nvars,(float (*)[3])Gmat,&nvars,ipiv,&info);

I also tried using int instead of mwSignedIndex, but I get the same error message.

The mwlapack library is not compiled with -g so I can't see my variables inside sgetrf.

Any suggestions?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Matthew 2013-7-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/82957-call-lapack-routine-from-inside-simulink-s-function-on-64bit#answer_92580

Here is a useful resource for understanding how to compile C code using mex on 64 bit devices: http://www.mathworks.com/help/matlab/matlab_external/upgrading-mex-files-to-use-64-bit-api.html?s_tid=doc_12b#bsgpkjp-3

The main pitfall is that the mex command uses the -compatibleArrayDims flag (32-bit API) by default. This results in the addition of the -DMX_COMPAT_32 flag in the call to mex. As mentioned above, this resulted in an error message stating that the array dimension I was passing in was incorrect: MKL ERROR: Parameter 1 was incorrect on entry to SGETRF. Merely changing the integers to a 64bit compatible type was not sufficient.

The fix was to include the -largeArrayDims flag in the call to mex and declare array indices as type mwSize. Since the S function is recompiled by the call to sim, it is necessary to set Simulation>Model Configuration Parameters>Device type> to the relevant 64-bit setting.

I don't understand why the -compatibleArrayDims flag is set automatically on a 64-bit machine. Perhaps someone can enlighten me.