MKL source code
I'm evaluating MKL for internal use on both servers and potentially a Phi cluster. I've encountered a number of frustrating bugs, poor documentation, and performance problems. Does Intel ever make MKL...
View ArticleAvoiding unnecessary OpenMP synchronization
Say I wish to add a number of vectors:cblas_daxpy(n, 1.0, a, 1, b, 1);cblas_daxpy(n, 1.2, c, 1, d, 1);cblas_daxpy(n, 1.4, e, 1, f, 1);cblas_daxpy(n, 1.6, g, 1, h, 1);MKL will use OpenMP to parallelize...
View Articlemkl_dfti.h: No such file or directory
I installed ipp library successfully setting its environment paths as described in the manual and finally ran a program with:# gcc -o main main.c -I /opt/intel/ipp/include -L /opt/intel/ipp/lib/ia32...
View ArticleSolving non-linear Equations
Hi!I've recently discovered how great and optimal everything in the mkl library is but my knowledge of it is limited. I want to solve a set of 12 non-linear equations in c++ and was wondering how to do...
View ArticleHELP! Problems encountered when using PARDISO
Hi,I'm trying to solve a huge sparse linear system using the PARDISO. I intend to use a subroutine named psolver to include all information and solution processes required when using PARDISO. First, I...
View ArticleRegister for Intel® MKL Webinar
Dear Customer,We would like to invite you to participate in Intel Math Kernel Library (Intel MKL) webinar.This webinar includes the following Topics:Intel® Math Kernel Library (Intel® MKL) on the...
View ArticleRegister for webinar: Four steps to accelerate Monte Carlo calculations
Monte Carlo simulations are crucial to financial risk management. The challenge today is running more simulations for more complex instruments at a lower cost. In a 30-minute webinar (September 27th,...
View ArticleNumPy ImportError: libmkl_avx.so: undefined symbol: ownLastTriangle_64fc
I'm running Ubuntu 13.04, and have the latest versions of the MKL, icc, and ifort. I'm trying to build NumPy 1.8.0beta2 and have been running into an unusual error. The package builds and installs just...
View Articledtrnlspbc_solve yields NaN
Dear all,For my application I use the bounded nonlinear solver almost identical to the examples provided in MKL.I show sample of my code (unfortunately it's not possible to presend a more detailed...
View ArticleClik here to view.
SPARSE BLAS CSR FORMAT
There seems to be no consistency with the CSR format in mkl. Using the function "mkl_ddnscsr", I convert a matrix, A, to csr format. It works great, because I put that into "dfeast_scsrev" and have no...
View Article/iface option cause mkl_dcsrmm error
I tried test the MKL subroutine mkl_dcsrmm on my 64-bit windows 7 with visual studio 2010 and Intel Fortran composer 2013. If I set the /iface as default, the subroutine mkl_dcsrmm worked fine. But...
View ArticleClik here to view.
Scalapack Bug in PZUNGQR?
Hi,Attached is a test case for some Scalapack calls that are exibiting strange behavior that seems like a bug. I also attached screen shots of the output and of the compile output. I'm using Windows...
View ArticleEfficiency of Pardiso
Hi All,Has anybody compared the efficiency of Pardiso with some other solvers using ILU preconditioner?In my test, Pardiso takes much more time than the solver using ILU preconditioner. For example,...
View ArticleWhy there is no LAPACK DGBTRF in new release of MKL?
I am looking for a routine in MKL that is equivalent to DGBTRF in LAPACK but only found DGETRF. I think DGBTRF was provided in older MKL, as disucess in this post...
View ArticleMKL 11.1 and scalapack / blacs_intelmpi
Hi,I have just downloaded the new compiler package with ifort 14.0 but it seems to be missing the fileslibmkl_blacs_intelmpi_lp64.so and libmkl_scalapack_lp64.soas well as thier static counterpart and...
View Article8-Byte Blacs Library and 4/8 byte mpi library
In some of our work (using the Intel Fortran compiler), we are linking the intel blacs against the intel mpi library and seem to be having some 4/8 byte interface issues. The documentation is not very...
View ArticleFast Fourier transform, MPI and linking (and alignment)
Dear all, Adding -mkl to my version (for details, see below) of ifort produces problems. This is why I read through the docs and used the intel link advisor about doing an MPI DFT over a computing...
View ArticleCorrect way to call Pardiso for optimum memory use
In 32bit I do this and it seems to work OK. Never using Windows virtual memory and not using OOC mode unless the problem is very big:Set MKL_PARDISO_OOC_MAX_CORE_SIZE to 1000Call pardiso with...
View ArticleOptimization/Parallelization in MKL
I have a 16 core system with hyperthreading (32 threads). The matrix calculations are on a 4096x4096 matrix. Using "dfeast_syev", it takes a few mins. And using "cblas_dgemm" it takes a few seconds....
View ArticleTrouble converting routine to use parallel RNG
HelloI have a code where I need to initialise a massive array with random numbers. For increasing problem size, this part of the program is becoming a drag. I want to use Intel VSL to do this, but I'm...
View Article