Does not dealllocate memory
Hello everyone, I have an application where I have to create some worker threads. These threads call a function inside a dll, where mkl is used.I have read some articles about memory leaks, and they...
View ArticleMKL BLAS causes Intel i9 9900K to crash/reboot(?)
Running under Linux Ubuntu 18.04 with standard/default BLAS library I observe normal behavior. If I install a parallelized BLAS library such as MKL or openblas, any simple matrix operation causes the...
View Articlecblas_daxpy gives wrong results with multuple threads
Hello,I have installed the most recent version of MKL and tried to use in my application built using VS 2015 C++. Certain parts of my code use OMP. However, I never call MKL routines from within an OMP...
View Articleweighted sum of two matrices
How do I do this using Intel MKL?res(:,:,:,i+j+k) = var(:,:,:,1)*op(1)+var(:,:,:,2)*op(2)In this example, 'res' is a matrix with 4 dimensions, of which the last one is only two.'op' is a vector two...
View ArticleSGEMV and SGEMM performance degradation on Xeon E5-2609
Hi,We have noticed a performance problem on a particular CPU (Xeon(R) CPU E5-2609 0 @ 2.40GHz) that we have narrowed down to a change in either SGEMV or SGEMM (or both) that occurred somewhere between...
View ArticleError making custom blacs library on MacOS with openmpi
I was trying to make the custom blacs library (after successfully doing this last year) and just received this error:sudo make libintel64 interface=lp64 Password:...
View ArticleVersions older than 2017 available?
Hello,I would like to use mkl 2016.1.111. On https://registrationcenter.intel.com/en/products I can find older versions, but the latest is from August 2017. Is it possible to get the mkl 2016.1.111...
View Articlecgemm3m, cgemm_compact AND cgemm give poor results for small problem 24*64
Hi all,I using sequential API and direct call to multilply matrices. C = 1*conj(A')*AA is 64*24 and C is 24*24 both are complex matrix (complex8).I have arrays of matrices: A_ARR (filled with random...
View ArticlePardiso matrix size change over the simulation
Dear all, I have just started learning PARDISO and I face a question:How can we use PARDISO when our coefficient matrix size changes over the simulation?I know that we can simply solve this issue by...
View ArticleAny routines for sparse symmetric matrix-matrix production?
Hi,I notice that MKL has special support for symmetric dense matrix production. Also MKL inspector-executor sparse BLAS routine supports sparse matrix-matrix production. My task is to conduct matrix...
View Articlemkl_jit_create_cgemm on multi cores
Hi all,I want to use mkl_jit_create_cgemm on my setup where each thread is pinned to single core.In each thread I'll do the cgemm with the created Jitter.Do I need to create jit kernel specific for...
View ArticleSome subroutines of lapack in mkl gives Segmentation fault for matrix...
OS: Ubuntu 18.04.2 server, kernel 4.15.0-50-generic. MKL provided by intel parallel studio xe 2019.3.The C code reproducing the error#include <stdlib.h> #include <stdio.h> #include...
View ArticleClik here to view.
MKL FFT Error in Example Code in Linux
Hello!I was trying to learn about the MKL FFT libraries and wanted to check how the example programs ran, but when I check the output files, it says there is an error, status = 2 while trying to create...
View ArticleClik here to view.
MKL FFT Fortran in Test Code - Transformed Values are not correct
Hello!I was trying to learn about the MKL's FFT Function and wrote a small 1D program to show the forward transformed values. From my understanding, the arguments for the forward function - mentioned...
View ArticleCannot use MKL F95 eigensolver routines
I am using Intel Parallel Studio XE with Visual Studio. I would like to use the heev()/heevr() routines. I have written a simple code to test it, and I have enabled the "Use MKL libraries" option in...
View ArticleDifferent LU factorization result between libmkl_intel_ilp64.a and...
I am a novice in the field of Intel Math Kernel Library. However, when I tried to use the LU factorization function from LAPACK by compiling different libraries, I got different results as the...
View ArticleDifferent BLAS performance between CMake dynamic and manuly static link
I tried to compile my program by using CMake at first and my CMakeLists file is as the following. cmake_minimum_required(VERSION 3.11)project(LMMNET LANGUAGES...
View ArticleHow to get Cholesky diagonal?
HelloI need to extract diagonal from cholesky LDLt sparse factorization. Can I do this with Intel MKL?I can use pardiso to factor and solve sparse linear system. But cholesky factor matrix itself is...
View ArticleBug in mkl_sparse_d_mm
Hi All,I think that the result of "mkl_sparse_d_mm" is incorrect as evident for the following code. The code basically tries to multiply the Identity matrix stored in CSR format with a dense matrix and...
View ArticlePossible issue in MKL LAPACKE interface when using LAPACK_ROW_MAJOR
MKL provided by Intel parallel studio xe 2019.3.The minimal reproducing code is below:#include <stdlib.h> #include <stdio.h> #include "mkl_lapacke.h" #define N 2 int main(){ double m[N];...
View Article