dss_factor_real export factorisation
Hello,I am working on a method where I need the factorisation of the coefficient matrix.Right now I call the routine dss_factor_real( handle, opt, rValues ) in C which creates a factorisation.Can...
View ArticleMKL for mac not available?
Hi, I downloaded l_mkl_2018.1.163.tgz but when launching it on my mac (macos 10.13.3) I get "The IA-32 architecture host installation is no longer supported."Google gave me an impression that I need...
View Articleinterface operator(s) for standard matrix math
Does the MKL software have available a Module with ready-to-use operators for typical matrix math? For example. Interface Operator (.x.) and Interface Operator (.tx.) ... A = B .x. C & A =...
View ArticleBounding threads to processors
Hi, This question is related to another post in this same forum (https://software.intel.com/en-us/forums/intel-math-kernel-library/topic/549769). We are using a code from there with minor...
View Article'mkl_dfti.f90' module error
I am trying to use intel mkl dft routines under windows/visual studio. Running intel mkl DGEMM routine presents no problem (hence automatic linking via seems to work). My problem is when I include...
View ArticleConvert dense matrix to sparse BSR
Hello,I am currently struggling to find a function that can convert a dense matrix to sparse BSR storage format. I know there is mkl_?dnscsr which does the job for CSR format but how about other sparse...
View ArticleError in mkl_sparse_z_set_value
Hello,I am struggling to change a single value of a matrix in internal Intel Sparse BLAS format using the mkl_sparse_z_set_value function. Most of the time this function returns a value of 3. However,...
View ArticleHow to multiply each row of matrix by a single column in CBLAS
I'd like to do the following in CBLAS:static float a[max_i][max_j]; static float b[max_i][max_j]; static float vec[max_i]; for (size_t i = 0; i < max_i; i++) { for (size_t j = 0; j < max_j; j++)...
View ArticleScalability of dense symmetric indefinite factorization
Hi,I am trying to speed up the factorization of a dense symmetric indefinite matrix, the size of my matrices is usually between 10 k and 20 k.I am using LAPACK (dsytrf) and MKL 2018 and I run it on a...
View ArticlevdSqr function seems not to return correct answer on Xeon Gold 6132
(Sorry for my poor English, I'm not a native.)Hello,vdSqrt function of MKL returns incorrect answer if I used this function under my developement environment.Here is my environment.OS:Windows 10 x64...
View ArticleMKL equivalent to IMSL lin_sol_gen
Is MKL/gesv an appropriate replacement for IMSL/lin_sol_gen? I believe both do an LU decomposition to get the solution. If I merely replace CALL LIN_SOL_GEN(A,B,X) with CALL GESVX(A,B,X), how close...
View ArticleMultiple jobs bug in makefile of FFTW3 interfaces
Dear MKL Team,I find that the current makefile of the FFTW interfaces, for example, /opt/intel/compilers_and_libraries_2017.4.196/linux/mkl/interfaces/fftw3xf/makefile does not allow multiple jobs...
View ArticleQuestion about mkl_?csrgemv
Hi, when I read reference manual of mkl, the mkl_?csrgemv routine calculates y=Ax (or y=A'x) where A is m-by-m matrix in 3-array CSR format.However, when I tested, I figured out that this routine also...
View Articlegetrf hangs when using intel threads in calling application
I am calling getrf from one thread in my application, the other threads are waiting at a taskwait. Task manager looks like all threads are maxed out and getrf never returns.The matrix size is quite...
View ArticlePossible Heap Corruption after calling DftiCreateDescriptor
Hello,Our application crashes some time after calling DftiCreateDescriptor; we believe because the heap has become corrupt. If we delete the call, the application no longer crashes, or if we reduce the...
View ArticleMemory usage during solve in function of number of cores
During the solve step of Pardiso of the system A*X=B where B is multiple right-hand-side we see that the memory consumed by Pardiso is highly dependent on the number of cores. So we observe mem(total)...
View ArticleSparse Matrix Multiplication Failed
Hello,I am trying to develop a function that multiplies matrix A and B which are in general format but sparse in nature. These matrices contain complex numbers. My issue is, when I do not use the...
View Articledss_factor_real export factorisation
Hello,I am working on a method where I need the factorisation of the coefficient matrix.Right now I call the routine dss_factor_real( handle, opt, rValues ) in C which creates a factorisation.Can...
View ArticleHow to multiply each row of matrix by a single column in CBLAS
I'd like to do the following in CBLAS:static float a[max_i][max_j]; static float b[max_i][max_j]; static float vec[max_i]; for (size_t i = 0; i < max_i; i++) { for (size_t j = 0; j < max_j; j++)...
View ArticleScalability of dense symmetric indefinite factorization
Hi,I am trying to speed up the factorization of a dense symmetric indefinite matrix, the size of my matrices is usually between 10 k and 20 k.I am using LAPACK (dsytrf) and MKL 2018 and I run it on a...
View Article