Question: cycle count of 65536 MKL FFT DftiComputeForward(C++)
My code as followings: fft_mkl(int M,float * InputData,float * OutputData){MKL_LONG status;DFTI_DESCRIPTOR my_desc1_handle; DftiCreateDescriptor( &my_desc1_handle, DFTI_SINGLE,DFTI_COMPLEX, 1, M);...
View ArticleDistributed Cholesky
I am doing a distributed Cholesky. You can find the code I am using (almost the same) here. I am gathering the submatrices (after every node has executed Cholesky) into the master node exactly as shown...
View ArticleComputation of the Schur complement in MKL
Hello everyone,I recently came across a project in which I have to compute the Schur complement of a complex symmetrix matrix. I know that starting from Intel® MKL 11.2 update 1, MKL supports the...
View ArticleCalling 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 Article