zgetrf/zgetri cannot invert matrix
Using code that worked with previous versions (Intel 12.1 OpenMPI 1.4.4 MKL 10.2.5.035 on CentOS 6) of MKL, zgetrf/zgetri cannot invert a matrix that it previously inverted.Could the newer intel...
View ArticleNNZ coefficients in dss/paridso LU factor
Hi,I was wondering whether there is any way to obtain the number of non-zero coefficients of the matrix factor generated by the MKL function mkl_dss_real? I checked the MKL manual but there seems to be...
View Articlepotential bug in mkl configuration in version 2018 update 2
I installed Intel® Parallel Studio XE with VS 2017 integration on a Win 10 platform.The parallel MKL option in the Properties page works fine as expected. But it seems like the sequential option...
View ArticleBandwidth of matrix
Hi,I'm using the DSS interface for solving a large linear set of equations. In order to do some analytical work on the performance of my implementation I would like to know the bandwidth of the system...
View ArticleScaLAPACK: pzheev error 706: (eigenvalue computation) needs MB=NB
UPDATE: sorry, my mistake, it IS in the documentation:* Alignment requirements * ====================== * * The distributed submatrices A(IA:*, JA:*) and C(IC:IC+M-1,JC:JC+N-1) * must verify some...
View ArticleIntel ode solvers
Здравствуйте. Можно ли узнать почему intel отказался от поддержки и развития библиотеки intel ode solvers ?
View Articleapt repository broken
Dear all,I get the following errors in Ubuntu 18.04 when runningwget https://apt.repos.intel.com/intel-gpg-keys/GPG-PUB-KEY-INTEL-SW-PRODUCTS... apt-key add GPG-PUB-KEY-INTEL-SW-PRODUCTS-2019.PUB wget...
View Articlepossible bug in 18.0.2 mkl
Hi there,the following code yield as segfault at run time:Program Test use lapack95 Implicit none Real*8, allocatable :: x(:,:) Integer*8 :: ISError if(allocated(x)) Then call...
View ArticleIntel® MKL version 2018 Update 3 is now available
Intel® Math Kernel Library (Intel® MKL) is a highly optimized, extensively threaded, and thread-safe library of mathematical functions for engineering, scientific, and financial applications that...
View ArticlePardiso - storing factorized array
Hello, I am using pardiso through a subroutine in a program (made with Intel Visual Fortran). I call pardiso with phase=12, then I call with phase=33 (to obtain the solution). There are some cases...
View ArticleNumpy and Scipy with MKL: installation troubles
Hi,I am currently trying to build Numpy and Scipy with MKL on Ubuntu 18.04 following your tutorial (but not using the Intel compilers). It seems that everything works for Numpy. Here is the...
View ArticleDefault value for mkl_set_threading_layer() in SDL?
I am using MKL with SDL (mkl_rt.lib). I know I need to use mkl_set_threading_layer() to set the threading layer in my code.My question is, if I don't call mkl_set_threading_layer(), is there a default...
View ArticlePardiso Low rank Update does not accelerate the decomposition
Hello, I have a problem: Pardiso Low rank Update does not accelerate the decomposition. The usual pardiso decomposition is faster than decomposition with low rank update. I have a complex symmetric...
View ArticleFGMRES ILU, Diagonal and RCI iterative solvers
I am trying to put together an FGMRES iterative sparse solver routine, that uses one of the examples given by INTEL, (I have tried many technics; the ILU, Diagonal preconditioner and the RCI; Reverse...
View ArticleClik here to view.
mkl_sparse_z_export_csr of Inspector-executor Sparse BLAS return bizarre...
Dear MKL experts, I'm testing the function,mkl_sparse_z_export_csr, in Inspector-executor Sparse BLAS, to convert COO format into CSR format and output the CSR format. The function returns...
View Articleconvert dense matrix to sparse CSR form
previously, I have used mkl_sdnscsr to convert a dense matrix to a sparse matrix with CSR format. However, with current update(2018 update 3), this function is depreciated, and the user manual instruct...
View ArticleMKL_DNN convolution has the wrong output order on Intel(R) Xeon(R) CPU...
Hello all,I recently implemented the convolution of the intel mkl library as described in the example included with the library. Everything is fine and dandy on my Laptop with a i5-3210M. However when...
View ArticleIs 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 Article