Intel® Math Kernel Library

mkl_ddiamv multi-thread support

Does mkl_ddiamv support multi-thread calculation? 

In my project's properties I set "Use Intel MKL" to "Parallel". The majority of time of my programm is matrix-vector multiplication, which is performed by mkl_ddiamv. In my task manager I always see, that ONLY ONE thread is busy by my programm process. I use Intel C++ Compiler 14.0. I try to take many different sizes of multiplicated arrays. There is no difference. I unsuccessfully tried to find some info about the reason of one-thread calculation. 

mkl poisson solver with periodic and dirichlet boundary condition

Hello Sir,

 I am trying to solve electrostatic Poisson's equation using Poisson MKL library 'd_Poisson_2D_c.c' but my boundary conditions are periodic (X direction) and fixed in Y direction. How shall i incorporate periodic condition in the X direction ? 

I tried in this way but i donot think it is correct,

for( iy = 0; iy <= ny; iy++){

  bd_ax[iy] = bd_bx[iy];


  and  applied BCtype='PPDD' .

I am very new in this field. Please help me.


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/ intel64 and source /opt/intel/mkl/bin/ 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.

Subscribe to Intel® Math Kernel Library