Calling LAPACK/MKL from parallel OpenMP region
Dear All,I often call some BLAS and LAPACK (MKL) routines from my Fortran programs. Typically, I try to place these calls outside of any parallel OpenMP regions while then making use of the parallelism...
View ArticleShould I place a barrier before calling pdpotri()?
I am using pdpotrf() in order to perform the Cholesky decomposition. Then I want to call pdpotri(), in order to invert the matrix. The function is called from every process, just after pdpotrf()....
View ArticleQuaternion Fourier Transform (QFT) in MKL
Hello,Would anyone know how to easily run Quaternion Fourier Transforms in MKL, or at least know of any other library that provides QFT on top of MKL's FFT? I have been using the Quaternion Toolbox for...
View ArticleBLACS broadcast 64-bit integer
Hi,I'm trying to broadcast a 64-bit integer with BLACS routines IGEBS2D() and IGEBR2D() --- Centos 6.5 Linux, ifort, composer_xe_2013_sp1.2.144, intel64, Intel MPI, ilp64 libs. Despite declaring all...
View ArticleIndices zero based or one based?
Are indices in the Intel MKL BLAS and LAPACK function zero-based or one-based. For example the function idamax I think returns a zero-based argument with MKL CBLAS. However, my experience with Netlib...
View Article[solved] no libmkl_tbb_thread.a after installation
Hello,I just finished to install a 2016 beta version of Intel parallel studio XE on centOS 7 (parallel_studio_xe_2016.0.035). I used the link advisor website to obtain the arguments list for linking...
View Articlemkl_?csrcoo duplicate COO values
Hi,I am using an old sparse linear solver in my code (of HSL) that receives the matrix in COO form and it accepts to have duplicate values. For example, (1,2,50) and (1,2,30) are two different entries...
View ArticleExploit symmetry in Cholesky
I am performing a Cholesky decomposition, like in this thread. Post #4 contains a minor example, that looks like as my code. The minor example of Ying is based on this example. As one can see, the...
View ArticleZYTRI sequential?
Hi I am testing the following code with composer_xe_2015lwork = 256*n call ZSYTRF( 'U', n, afull, n, ipiv, work, lwork, error ) call ZSYTRI( 'U', n, afull, n, ipiv, work, lwork, error )usinggfortran...
View Articlewhy defferent E5 v3 and E5 v2 in linpack
When I use intel mkl , intel mpi ,intel compiler to test linpack l in e5 v3 2660, I get RMAX/RPEAK IS 11% .BUT in e5 v2 is 93%. Why?
View Articlelinpack different in v3 and v2
When I use intel mkl , intel mpi ,intel compiler to test linpack l in e5 v3 2660, I get RMAX/RPEAK IS 11% .BUT in e5 v2 is 93%. Why?
View ArticleError -2 when running PARDISO with multi threads
Dear all,I am trying to solve a sparse symmetric indefinite real matrix using PARDISO. The matrix is of size n=55,650,566 with 164,833,721 nonzero entries. Everything is good when the number of enable...
View Articlegetting started on MKL - Error in opening the compiled module file.
I've freshly installed the Intel® Parallel Studio XE Composer Edition for Fortran OS X* (student version). It comes with the Math Kernel Library, which is why I bought it. I'm having a hard time...
View ArticlePARDISO Phases called in different routines,Mem leak?
HiI am calling pardiso with c/c++. for that I packed all phases of pardiso in a library , but splitted into different soubroutines using *pt and *iparm as global variables within that file/library...
View ArticleFail to preform DFTI.DftiComputeBackward on 64 bit platform (windows)
Hello.I've built an example application in C# using visual studio 2012. The program do FFT for the vector {1,2,3,4,5,6,7,8,9,10} and then reverse FFT and checks if the return values are the same as the...
View ArticleMath Kernel Library 10.0 Installation bug
Hello,Sorry if this is the wrong place to post this, but I needed some help. I'm attempting to install the Math Kernel Library on Mac (Several different OS versions. Yosemite being the primary test)....
View ArticleClik here to view.
cblas_dnrm2 much slower than cblas_ddot
Dear all,I run benchmarks on a sandy-bridge Intel processor (E5-4620) using Intel MKL 11.1. Here, I have found that cblas_dnrm2 is significantly slower (3.4 s) than the corresponding cblas_ddot call...
View Articlethreading for mkl_?bsrmv routines?
Linux, MKL 11.xAre the mkl_?bsrmv routines threaded? We're seeing threading if we use mkl_?csrmv routines. But when using mkl_?bsmrv routines we do not see threading.Link line is...
View ArticleBlocks of different sizes in ScaLAPACK?
I am performing a Cholesky factorization with Intel-MKL, which uses ScaLAPACK. I distributed the matrix, based on this example, where the matrix is distributed in blocks, which are of equal size (i.e....
View ArticleMultidimensional DFT and OpenMP
I'm working on a program that performs several 3 x 3d (N1xN2xN3) DFTs using the MKL DFT algorithm. I'm running most of the program in parallel using OpenMP and I'd like to get as much parallel...
View Article