Is DTRSM a sequential code in the case of a single right hand side?
Hello, DTRSM is the LAPACK function which I use to solve triangular linear system with a multiple right hand sides. I am perfectly happy with the performance and the parallel scalability of the...
View ArticleMKL cluster solver: number of threads always equal 1
Hello everyone,I'm using the MKL cluster sparse solver but I found a problem related to the number of threads. In fact the number of threads for my example is kept equal to 1 for every calculation and...
View ArticleSparse matrix complex number
We are going to solve linear systems with sparse matrix of complex numbers and we would like to know, which is the most suitable mkl routine to use. We read the documentation but we have still some...
View ArticleClik here to view.
I get different results from MKL and LAPACK.
Has anybody experienced getting different results from -mkl and -llapack?I've developed a fortran code that uses zgeev as inLWORK=-1 call...
View Articlenon-negative matrix factorization
Hello,does anyone know if the non-negative matrix factorization is available in mkl; or will it be implemented in mkl in the near future?Thanks!
View ArticleImpossible to use MKL in Visual Studio 2017
Hello, I know that this topic has been widely addressed but in my case I have tried with all the proposed solutions that I have found and no one worked for me. The testing code is the following://...
View ArticleCNR Support Functions in Intel MKLML small libraries
Hi,I'm using the Intel MKLML small libraries shipped with MKLDNN [1] to build a custom inference engine for neural networks. It's great to get MKL performance in such a small package.However it appears...
View ArticleMKL 2018 Update 3 (Windows) zgetri crash
MKL 2018 Update 3 has broken zgetri.Our QA process now crashes in a call to zgetri. Reverting to MKL Update 2 DLL's resolves the issue. It does not happen on every call to zgetriI will point out that...
View ArticleStrange partial pivoting of LAPACKE_dgetrf
Hello,I'm using LAPACKE_dgetrf to compute the LU factorization of square matrices in double precision. The matrix is in column major. Here is what I am doing. The environment is MKL 2018 Update 3 for...
View Articlelargest matrix size by the LAPACKE_dsyevr
Hi, I am trying to diagonalize a square matrix of size N=73789, with a MKL LAPACKE_dsyevr, and need to find all the eigenvalues and eigenvectors. The machine I am running on uses Intel(R) C++...
View ArticleClik here to view.
PDGEMM returns wrong result under specific conditions
Hello,I encountered a problem using Scalapack routine PDGEMM under specific conditions. After investigation, I managed to reproduce the problem with a simple test program. In my case, the error occurs...
View ArticleShared library using MKL crashes on runtime with EXC_BAD_ACCESS
Hi everyone!I'm trying to create an application for MacOS platforms that uses some MKL methods for its math computations, namely the `dfdInterpolate1D()`. I'm facing a difficult issue, with random...
View Articleiparm[35] = 1 or 2 question
Dear IntelI'm a user of MKL 2018 up1 parallel studio.Usually I handle the large matrix(>50Mx50M) using schur complement.'iparm[35] = 1' is the my general setting and it shows the schur complement...
View ArticleMust the Block in sparse BSR format be Square?
Dear all,I am new to the MKL sparse blas, and want to use BSR as my sparse matrix format.I notice that the block_size is specified by an integer. Does this mean the block must be a square matrix?Is...
View Articledyld: Library not loaded: @rpath/libmkl_intel_lp64.dylib when running scripts
I tried to run a series of simulations using GROMACS 4.6.5 on MAC PRO in which intel composer_xe_2013 was installed. If I manually submit an individual job using commands in a terminal, I do not have...
View ArticleClik here to view.
Potential Buffer overwrite in DGEEV when compiled in a debug mode.
Hi, We recently upgrade our MKL to version 18.3 and TBB to version 18.4. When running a debug version of our code we noticed that we were getting heap corruption when calling the DGEEV MKL function. It...
View ArticleMKL GEMM slower for larger matrices
For matrix mul A(m,k) * B(k,n):m=9, k=256, n=256 is faster than m=9, k=512, n=512 and all larger k and n.On my E5-2630v3 (16 cores, HT disabled), k,n=256 get 850 GFLOPS while k,n=512 only get 257...
View ArticleSegfault when using cblas_dgemm()
When attempting to use the cblas_dgemm() function, I am experiencing a segmentation fault. The segfault occurs only when running the code on a Linux machine (works on a Windows machine).The code is set...
View Articleside change for mkl_dcsrmm
Hi allthe current version of mkl_dcsrmm implements the operationsC=AB+C and C=A'B+C, where in both cases A is the csr sparse matrix.Is there any way/alternative such that B is the csr sparse matrix??In...
View ArticleProblems with environment variables and using code examples
I was following installation guide up to environmental variable step: https://software.intel.com/en-us/mkl-windows-developer-guide-setting-env.... I entered mklvars intel64 in the command shell. The...
View Article