Clik here to view.
DGELSS Issue
I am attempting to use this routine but am getting an exception:Exception thrown at 0x01FA5428 (mkl_core.dll) in Console1.exe: 0xC0000005: Access violation reading location 0x00000000.I wonder if...
View ArticlePARDISO, How to escape the occurrence of zero pivot error
Hi everyone,I have zero pivot error -4, but I know the error occurs in the condftion of my calculation.So I want to have a solution using zero pivot replacement (the value i_parm(10)) and without...
View Articlemkl_set_num_threads() doesn't work
Hello everyone!I am using MxNet library with MKL support. I need my mxnet predictor use only one thread for calculations. For this purpose I used OMP_NUM_THREADS = 1 environment variable. So it works...
View ArticleRare crashes on MKL
I have implemented in C++ an algorithm in image processing using (among other things) fftw wrappers in MKL library (version 2018.3.210)I am working on a x64 machine with Intel Xeon E5-1650 v3 3.5 GHz...
View ArticleQuestion about the mkl_?omatcopy
I can not understand the manual about the following parameters of the function.rows The number of rows in matrix B (the destination matrix).cols The number of columns in matrix B (the...
View ArticleHow to access the number of non-zero elements in sparse_matrix_t?
I am writing a Python wrapper for calling the 'mkl_sparse_spmm' function. In order to export the result of matrix-matrix multiplication to a Python object, I need to know the size of the 'col_idx' or...
View ArticleIntel Pardiso error numerical factorization
Dear Pardiso users,I am using Intel Pardiso to solve a sparse system. For small matrices it works perfectly and very fast. However, with increasing system size I figured out that for some set of...
View ArticleFinding the eigenvalues (diagonalizing) of a block-diagonal matrix
I have to diagonalize a large matrix which take a lot of time. The matrix size 10,000 x 10,000.This matrix is Hamiltonian of a spin system which have some block structure. Is there way to diagonalize...
View ArticleUsing FMA in MKL routines
Hey everyone,I couldn't find any old topics that dealt with this question in detail, so here I am asking it again: is there a way to enable FMA math when using the MKL routines? Here is a sample...
View Articledfeast_sygv -4 error, BUT B IS POSITIVE DEFINITE!!!
Hello all,I am using the extended eigensolver routines specifically dfeast_sygv function, and I get the following error==>INFO code =: -4Intel MKL Extended Eigensolvers Error: Matrix B is not...
View ArticleHow to compile and run Quantum Espresso with intel MKL and Open MPI on Mac OS X
I want to run Quantum Espresso faster on Mac Pro (multi-core). I think using intel MKL(+FTW3) and Open MPI with intel compiler is the best solution. It was a process of trial and error, and I finally...
View ArticleClik here to view.
MKL 5 times slower when called from a MEX function in Octave
Hi,I'm trying to use Intel MKL library in an Octave MEX function, but the performance that I achieve using some MKL functions such as cblas_cgemm is 5 time slower when called from Octave rather than a...
View ArticleMemory Leak in MKL
Our program analyzes several images every few seconds. Each image is processed in its own thread. The threads are ordinary Windows threads. MKL functions are used in the process. We used a memory leak...
View ArticleTypos in the mkl_spblas.h?
Some of the function declarations in mkl_spblas.h do not agree with the documentation. In particular, in the create and export routines, the parameter names for the CSC format coincide with those for...
View Articlebad interaction betwenn mkl_dynamic and potri in mkl 19.02???
Hithe program below implements the inversion of an autoregressive matrix.Program Test use blas95 use lapack95 USE IFPORT use mkl_service implicit none integer(kind=8) :: istat, n, c1, c2, ise...
View ArticleKeep some values of vector.
Hello, I want to keep some elements of a vector without using a for-loop. For example,x = [ 10 20 30 40 .... 100] and index = [1 5 10];In MATLAB you use x(1,index). Could you please tell me if there is...
View Articlepardiso out of memory (-2) in phase 33
Hi,pardiso stops with error message -2 in phase 33 when increasing the number of rhs from 600 to 700. 600 rhs runs fine with a memory use of 146.2GB ("RES") reported by "top". The system has 256GB of...
View Articlesome worry about FFT's half result
Hi,Intel masters, I am writing to ask for information about the FFT's half result . In matlab,i use fft to compute an real array ,and then use ifft ,i'll get the initial array.However , in MKL, i just...
View ArticleQuestion about matrix inverse
Recetnly, i was trying to replace the function for matrix inverse by that of MKL. However, after making it, i found that the time used for computing was not improved and be slower. I don't know the...
View Article