Библиотека Intel® Math Kernel Library

Pazzled with pdpotrf arguments

Recently I did a Cholesky decomposition with dpotrf(). My goal is to do that in a distributed environment, so I need to use pdpotrf(). However, I can not understand what values should ia, ja and desca. Can anyone please explain, since I did not find a C example about that?

Adding MKL to Path

I just installed Parallel Studio and need to make use of the C compiler, Fortran compiler, and math kernel library. After installing, I add source /opt/intel/bin/compilervars.sh intel64 and source /opt/intel/mkl/bin/mklvars.sh intel64 to my bashrc file (I'm running Fedora linux). When I go to check that everything has correctly added to my path I run "which ifort" and "which icc" and see the location of the ifort and icc executables. However, this method doesn't seem to work with MKL given that I don't believe there are executables for which I could run "which <blank>".

mkl_ddiamv example

I have some problems with using mkl_ddiamv (language C). Can anybody show me a code example of using mkl_ddiamv with initialization of all the input formal parametres? It is desirable to use the number of diagonals more than one. For example, 3-diagonals with offset=[-1 0 1] and size 5 by 5 with all elements are ones.Multiplicated vectors are ones too. In different sources the matrix in special diag format is treated as two-dimensional array or one-dimensional array with row-by-row storage.

dpotrf arguments

With the help of this great forum and old post, I was able to finally execute my first example on Cholesky factorization. However, I have a question. I call the function like this:

        &N, // A is NxN matrix

Since I am planning to do a distributed factorization, I have to be sure about the arguments. What is the latter 'N' used for? My first thought was about the dimension of the output, but that's the same array, so I am confused.

Simultaneous diagonalization of commuting matrices

Dear Intel MKL colleagues,

May I ask please to someone let me know, if there is any built-in function or easy to follow process for simultaneous diagonalization of multiple (or even two) commuting matrices in Intel MKL ?

Any help would be appreciated.

Best regards,



segment fault error for fortran code

hello, everyone, 


Newbie of the fortran code.   I was trying to use ilu0+gmres to solve Ax =b with A the csr format. But got the segment fault error when calling the subroutine. It seems that the arrays didn't pass successfully. Really confused about the error. Would anyone please help me about that?    Thanks!


 The variables are defined as : 


      integer nn, nnz




Example does not run

I am trying to run ScaLAPACK, through MKL. Since I have the student version, I had to install MPICH2, since what I want to do is to apply an inverse of a matrix, that can be done with Cholesky in a distributed environment.

I got this ex1.c from Intell 11. It compliles fine, but it fails to execute:

Kronecker product


I am new to Intel MKL and I am sorry if my question seems out of line.

Is there any function or subroutine in MKL that computes the kronecker product of two matrix (like kron in matlab).

Also is there a function in MKL that returns the identity matrix of a specified dimension (like eye(n) in matlab).

Thank you!

Подписаться на Библиотека Intel® Math Kernel Library