dggev deterministic behaviour
Hello, I am using dggev (mkl versions 11.2.1 & 11.2.2) and i notice some non deterministic behaviour concerning eigenvalues near zero and infinite eigenvalues and eigenvectors. I use the parallel...
View ArticleDoes Intel MKL Pardiso come with selected Matrix Inversion
Hi, I am trying to compute the inverse of a sparse matrix. I saw in the pardiso website that manipulating certain parameters may result in computing the inverse of selected rows and columns of the...
View ArticleClik here to view.
Inplace permutation
Hello,I need to permute a vector according to an index array jpvt returned by ?geqp3. In the example lapack/source/dgeqpfx.f it is done using an auxiliary array. I would like to do the transformation...
View ArticleClik here to view.
ZGETRS memory corruption(?) with denormal numbers
We recently encountered a rather strange issue when passing denormal numbers in the right-hand-side matrix to ZGETRS. Attached is a small C++ file that reproduces this error.In this file, we set up a...
View Articlebug of sgemm, 2015 mkl gives wrong result, but 2013 gives correct
The compile option of 2015:/work1/soft/intel2015/composer_xe_2015.0.090/bin/intel64/ifort -i8 -openmp main.f90 -o i64_2015 -L /work1/soft/intel2015/composer_xe_2015.0.090/mkl/lib/intel64...
View Articleconversion from COO format to CSR
Dear all,I am using Intel MKL subroutine mkl_zcoocsr to convert a sparse matrix in COO form to CSR, but the results do not seem right. Here is my...
View ArticleClik here to view.
A bug in vsldConvExec
I use 64-bit AMD Linux machine. When I run the test file (attached) I get wrong results with the current MKL 11.2 update 2:EXAMPLE executing a convolution task result: z[0]=1343342.8 result:...
View ArticleClik here to view.
QR factorization without pivoting
Hello,I'm trying to compute the solution of a singular square system using a QR factorization. When I'm using a QR factorization with pivoting (i.e. ?geqp3), everything works OK. However, I'm concerned...
View ArticleClik here to view.
Strange (and seemingly wrong) result from Pardiso
I obtained strange (and seemingly wrong) result from Pardiso. I used Pardiso to solve a very ill-conditioned complex linear system arising from a circuital analysis of a three-dimensional conducting...
View Articleproblems setting MKL environment variables
Hello,I have searched the forums for a solution to my problem but haven't seen a similar topic.I have been unable to get the scripts that set the mkl environment varaibles (described here:...
View ArticleInput error in Pardiso error_num= 15
Hi! I am a new user. When I try to compile my fortran code in Pardiso. I got some errors. *** Error in Pardiso < sequence_ido,parameters> error_num= 15 *** Input check: matrix_type_new 0 <out...
View ArticleIntel® Math Kernel Library (Intel® MKL) 11.3 Beta and Intel® Parallel Studio...
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 ArticleHow to config the PARDISO OOC file and problems with SORTQQ
Hi,I have a question on how to config the OOC file for PARDISO in order to use the OOC mode.I'm trying to solve a very huge system, may take up to 500G RAM. I have a Xeon workstation with 250G RAM and...
View ArticleHow to set affinity of threads spawned by MKL?
I'm trying to do some optimization on a NUMA system at the moment and have run into trouble with regards to thread affinity.The essence of the code I'm using can be seen below: success =...
View ArticleSignificant Overhead if threaded MKL is called from OpenMP parallel region
Hello,my aim is to diagonalize quadratic matrices with different sizes dxd in parallel. To this end I wrote a for loop. In each iteration the aligned memory (dependent on the dimension d) is allocated...
View ArticlePARDISO vs. RCI CG
Hi, I am trying to solve a relatively large system (100.000 equations) using either Pardiso or CG.The system matrix is sparse, symmetric and converted to CSR format. For the matrix-vector...
View ArticleSlow rectangular matrix transposition ?
Hello, I'm working with MKL 11.2.0.090 on Gentoo. I have an "Intel(R) Core(TM) i7-2600 CPU @ 3.40GHz" processor. I'm trying to speed my inplace matrix transpositions and for that I thought that...
View ArticleIs there a Vector/Matrix Fill/Set?
Hello,Is there a function to set every element of a vector/matrix to a constant value? You can do it with mkl_calloc, but there a cases where you want to set an existing array and I was wondering if...
View ArticleComplex 1-D DFT not respecting DFTI_THREAD_LIMIT
Hi,I recently noticed that when using a threaded 1-dimensional DFT, a DFTI_COMPLEX domain DFT does not appear to respect the DFTI_THREAD_LIMIT and instead always uses the threading value set by...
View ArticleInstallation fails --Jave Class Not Found -- on MacBook Pro
I am trying to install the Math Kernel Library under: Product Subscription Information Download Latest Update Release PostedAcademic Research Performance Libraries from Intel (OS X*) On my MacBook...
View Article