sparse pblas?
Hello everyone,I am trying to use MKL to solve sparse matrix multiplication (i.e. C=A*B) in parallel. I found there is sparse blas for sequential code but I am wondering if there is a parallel version...
View ArticleCGEMM performance strangeness on Haswell CPUs vs. Sandy Bridge
Hi AllI have investigated the performance of the CGEMM algorithm using both my own sandy bridge CPU and my colleagues newer computer with a Haswell cpu. The calculation is measured as the number of...
View ArticleCompilation error
HiI am trying to compile the C programme which uses the MKL subroutine.For example, mkl_scsrdia and mkl_scsrcoo,but when I am trying to compile it, I am getting error like this. It seems that I am not...
View ArticleAnnouncing new product: Intel® Data Analytics Acceleration Library 2016 Beta
We are pleased to announce the release of Intel® Data Analytics Acceleration Library 2016 Beta! Intel® Data Analytics Acceleration Library is a C++ and Java API library of optimized analytics building...
View ArticleReduced Row Echelon Form of a matrix (rref)
Hi,I'm newbie to the Intel MKL library and I'm trying to convert a code from Matlab to C using the C interface of Intel MKL routines.I haven't found a function to achieve the RREF of a matrix in Intel...
View Articlewhat is ARR_HDR_SIZE in IPP samples
Hello All,I am implementing the JPEG encoding and everything looks fine. But, I want to know the ARR_HDR_SIZE in uic samples #include <stdlib.h> #include <etxt.h> #include...
View Articledoes MKL FFT use pre-calculated twiddle factor table or calculates the...
Hi, I'm evaluating the MKL FFT performance, I have a quick question, does MKL FFT use pre-calculated twiddle factor table or calculate the twiddle factor parallely when doing FFT? Thank youJohn
View ArticleOpenMP not using all processors
I am trying to use MKL libraries and OpenMP in a MSVS C++ application on Windows7. The application shows affinity for all 24 processors (2 nodes, 6 processors, HyperThreaded). omp_get_num_procs() also...
View ArticleOpenMP not using all processors
I am trying to use MKL libraries and OpenMP in a MSVS C++ application on Windows7. The application shows affinity for all 24 processors (2 nodes, 6 processors, HyperThreaded). omp_get_num_procs() also...
View ArticleClik here to view.
Modified Cholesky factorisation
Hello,I'm using MKL to calculate Cholesky factorisation of a covariance matrix. MKL (?POTRF function) is of course much faster than my own naiive implementation (input: 6500x6500 matrix), however there...
View ArticleClik here to view.
MKL's (distributed) FFT library fails with a floating-point error
When repeatedly calling MKL's distributed (cluster) DFT library via the FFTW3 interface library, it will fail with a floating-point error with certain combination's of grid sizes and MPI processes...
View ArticleWhat is the calling convetion for i_malloc and friends on 32bit Windows
Hello,I am trying to redefine memory allocation functions for MKL by using function pointer i_malloc and friends as described in https://software.intel.com/en-us/node/528395. On 32bit Windows however I...
View ArticleCalls to dcscmv lock
Calling dcsmv is causing my program to lock -- periodically when I invoke dcsmv, the call never returns and the process appears with a Sl+ status when I run ps. When I launch the program in a debugger...
View Articlehow do I know MKL FFT being called used AVX-512
Hi I am trying to look at MKL FFT performance by calling 5 lib functions as below, I configured the optimization and enabled AVX-512 etc in the project properties setting(vs2013 integrated with...
View ArticleIs it save to call ?feast_syev/?feast_heev from multiple threads?
Hi,I'm developing an application that needs to compute various eigenvalue decompositions. Is it possible to call zfeast_heev from multiple threads in parallel? Ofcourse, each thread has it's own...
View ArticleParallel Direct Sparse Solver for Clusters in Windows OS
Is the Parallel Direct Sparse Solver for Clusters supported on Windows OS?
View Articleoptimizing dgemm on NUMA systems
HelloI have been trying to optimize matrix multiplication on NUMA systems but so far without much luck.I have played around with the dgemm routine and first touch.A snippet of my code looks like this:...
View ArticleNumerical integration package in MKL
Hi, I have "old" codes that call Fortran numerical_libraries routine GQRUL and DGQRUL for calculating Gauss-Legendre quadrature rule to perform numerical integration. I used to be able to just put a...
View ArticleDSS symmetric solve
We have been using the following code compiled with Intell Fortran on Linux clusters linked with MKL 10.2.2.025 for a number of years. We are trying to update to Intel 15 compilers and MKL 11.2 update...
View Articlemkl_pardiso time consumption during solution phase
Hello, I am working with a sparse matrix A and trying to solve <Ax =b>. b is an array of 1 million and it is in double precision. The matrix is sufficiently sparse (number of non-zeros in A (%):...
View Article