Recent posts
https://software.intel.com/en-us/recent/541742
enLAPACKE_zsptri not threaded?
https://software.intel.com/en-us/forums/topic/338111
<p>Does anyone know if LAPACKE_zsptri is threaded in MKL 10.3? It doesn't appear to be, but I may be doing something wrong. It isn't listed in the user guide as being threaded, but it also has a list of functions and then says "and so on", so I was hoping it was part of the "so on". If not, is there any particular reason it isn't threaded?</p>
Thu, 15 Nov 12 10:18:52 -0800scossmann338111creating a MATLAB .mex file using dspsv
https://software.intel.com/en-us/forums/topic/280634
<p>I'm trying to use the MKL libraries to build a MATLAB routine to solve a matrix equation with a packed storage symmetric matrix (to save memory). I'm able to get it up and running fine using the following code, but for larger problems (say 5000x5000, which is still relatively small) it will run fine ~50% of the time, and crash MATLAB ~50% of the time. Does anybody see anything that I am obviously doing wrong? #include #include #include typedef size_t INT; #define MKL_INT INT#include "mkl_lapacke.h"void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]){ /* Locals */ lapack_int n, nr, lda, ldb, info; double *ar, *br; ar = mxGetPr(prhs[0]); /* pointer to first input matrix */ br = mxGetPr(prhs[1]); /* pointer to second input matrix */ // dimensions of input matrix A // lda = mxGetM(prhs[0]); n = mxGetM(prhs[1]); // dimensions of input matrix B // ldb = n; nr = mxGetN(prhs[1]); if (lda != (ldb*(ldb+1)/2)) mexErrMsgTxt("Error, size mismatch."); // Local arrays // lapack_int *ipiv; ipiv = (lapack_int *) mxMalloc(2*n*sizeof(lapack_int)); //lapack_int ipiv[N]; info = LAPACKE_dspsv( LAPACK_COL_MAJOR, 'U', n, nr, ar, ipiv, br, ldb); mxFree(ipiv);}</p>
Wed, 07 Dec 11 07:00:19 -0800scossmann280634