Use ?gbmv with a 3d array.
Let assume the we have a 3d array A and a matrix B of size 2x5. int A[3][2][2]= { {{1,2},{3,4} }, {{5,6},{7,8}}, {{9,10},{11,12}} }; I want to program the following:cs =0.0;for j=1 until 3Multiply the...
View ArticleExtract row/column of a matrix and subvector form a vector.
Hi! I would like to ask if there are routines for the following issues:1) Extract a row/column of a matrix.2) If theta is a vector extract a subvector. For example extract the elements 4 until 8 of...
View Article3D convolution
I am using MKL's VSL multi-dimensional convolution routines to perform a 3D convolution. As a test, I am using the input function as a 3D gaussian as follows:f(x,y,z) = exp(- (x-2)**2) exp(-(y+2)**2)...
View Articleselling R / Python pre-built with Intel MKL
Just a thought: instead of Intel selling MKL and then letting the clients build/extend R or Python with Intel MKL, it could be a lot easier if Intel could sell corresponding versions of R or Python...
View ArticleIs MKL plan to implement the parallel sparse blas?
I know the MKL can deal with the sparse blas. It is very useful for the finite element coding. However, for the parallel computing, it needs parallel sparse blas. So, Is MKL plans to implement the...
View ArticleProblem with MKL and MEX
I am trying to write a mex programm but it gives errors. When I use the VS2012 the code runs, but when I converted it to mex it didn't work. The code gives errors when it tries to use the mkl library...
View ArticleScalapack PXLAWRITE and PXLAREAD
Hi,Does MKL Scalapack include the PXLAWRITE and PXLAREAD (X=Z,C,D,S) subroutines? When I try to use them in my code I get 'undefined reference' link errors. I do not see them in the MKL...
View Articlefeast threading
hi guys, only getting 50% throughput CPU-wise when using FEAST's eigensolver. i am aware that hyperthreading is merely emulation, but isn't it possible to get the solver to still exploit all available...
View Articlefortran+fftw=nan (sometimes)
Hello,I run a fortran multithreaded code that uses ffts. I am now trying to switch tofttw available thru mkl. When I insert the code below the results moreoften than not are nans or infinity even...
View ArticleMKL from C
I am having trouble understanding C interface for MKL. In particular the const modifier. I need to use tridiagonal solver ?dtsvb which should have C interface:void ddtsvb (const MKL_INT * n, const...
View ArticleALLOCATE or MKL_MALLOC?
In Fortran, assuming I have a 64bit avx-able machine at my disposal, Is ALLOCATE statement guaranteed to allocate aligned arrays when the compiler is called with the "-align array64byte" flag? Or, do I...
View Articleprojection operation in MKL
Dear all,For a projection operation, I should implementF = ( I - Mf*Vfr*Vfr^T ) * FwhereMf is a sparse matrix of size sz_fXsz_fVfr is a vector of sz_fX1F is a block of vectors of size sz_fXP, and...
View ArticleMKL 11.0 performance problem on Opteron
Hello,I am observing a weird performance behavior when using MKL 11.0 on an AMD Opteron 6220 processor. I wrote a small benchmark that measures the Gflop rate of the _GEMM operation in S, D, C and Z...
View ArticleForum poll: Intel MKL and threading
Intel MKL users,We would like to hear from you how you are using Intel MKL with threading. Do you use the parallel or sequential MKL? How do your multithreaded applications use MKL? We would appreciate...
View ArticleCorrelation
Hi all, I have a lot of trouble understanding how the correlation/convolution functions are used. I have looked at the Information from the Manual but the description - involving finite functions - is...
View Articlepdlacpy does not work with matrices distributed differently.
I am using pdlacpy. Documentation says this function copies all or part of a distributed matrix (A) to other distributed matrix (B). It also says no communication is performed (e.g. the function...
View ArticlePerformance difference bewteen i7 and xeon CPU for intel mkl
Hi~ I think i might meet with some wired problem with Intel MKL.I install visual studio 2013 and intel xe composer 2013 sp1 on two machines, one is pc with windows 8.1 equped with intel core i5-3470...
View ArticleWebinar announcement: A Tour of the Sparse Linear Algebra Functionality in...
Sparse matrix algorithms are encountered in a broad range of important scientific computing applications. Intel MKL offers a powerful set of functions that can be used to build a complete solution to...
View Articlefeast threading
hi guys, only getting 50% throughput CPU-wise when using FEAST's eigensolver. i am aware that hyperthreading is merely emulation, but isn't it possible to get the solver to still exploit all available...
View ArticleIssues with executing PARDISO, dcsrilu0_example0 MKL routines on Xeon Phi...
Dear Mic people,I am facing problems executing Intel MKL routines, PARDISO and dcsrilu0, on the Xeon Phi co-processor.Please go through the report below. I would be very grateful if you could tell me...
View Article