FFT Speed up Generation
Hey thereI d like to use the MKL library to generate the 1D FFT(complex to complex), with the size of 256.To speed up the process i want to split the 256 FFT vector to 4 sub-vectors of 64 and run each...
View ArticleHelp using mkl_?imatcopy
The mkl_?imatcopy functions require an extra leading dimension parameter, called ldb. From the documentation, I am not sure what value it should be used for this parameter. As an answer, I would like...
View ArticleFortran intrinsic functions v.s. mkl functions or subroutines
Hi all,I am wondering which I should use in my code, for example if I do matrix multiplication A(100,100)*B(100,100), matmul(A,B) or gemm()?The same uncertainty for other functions, e.g. dot_product,...
View ArticleRenaming libiomp5md.dll for 64 bit
Hello. I want to use both the 32 and 64 bit versions from the same folder. How can I rename one of them to avoid conflict?Alternatively, I could statically link on one platform and dynamically on the...
View ArticleRenaming libiomp5md.dll for 64 bit
Hello. I want to use both the 32 and 64 bit versions from the same folder. How can I rename one of them to avoid conflict?Alternatively, I could statically link on one platform and dynamically on the...
View ArticleBLACS_EXIT doesn't call MPI_FINALIZE
I wanted to do some calculations with ScaLAPACK. To introduce myself to the API, I started with the"Hello World" example from Intel MKL Reference Guide.While running this code with a square number of...
View ArticleSkipping empty rows in mkl_dcsrmv()
In my application, I use mkl_dcsrmv(). I encounter matrices with empty rows. At most, one eighth of the rows are empty. Is there any solution to skip empty rows in CSR-based sparse matrix-vector...
View ArticleLink MKL 9.0 to VC++ 2008 express edition in Windows 8
Hi. I need to link the MK libraray to my project. On Windows Vista (32 bit) it was somehow easy:1) goto Tools->Options->VC++ Directories and add 'C:\Program Files\Intel\MKL\9.0\include' in...
View ArticleClik here to view.
Sparse Matrix Structure leads Zero Pivot Error
Hi all,I am using the DSS from MKL to solve the linear system of equations Ax=b, for which the resulting matrix A is sparse and not symmetric, but its structure is symmetric.Concerning the options, I...
View Articledss_solve_real takes more time to solve a linear system
I need to solve a system of linear equations with three righ-hand-side vectors. Initially, I was using the sequential version of MKL (compiling with "libmkl_sequential.a") and solving for each rhs...
View Articlefrom CSR to CSC for non-square (m x n) matrix
Hello everyone. I need to convert a non-square sparse matrix (m x n) from CSR format to CSC. But MKL has routines only for square matrix (m x m): mkl_?csrcsc. What can I try to make to complete this my...
View ArticleDoes MKL and PARDISO work properly with more than 64 cores?
This is a continuation of my old question "Does PARDISO have core access limit" with new information. I'll put a quick summary and the new question:I run a code on Windows based HP DL980 server with 80...
View ArticleIntel® Math Kernel Library 11.1 is now available
Intel® Math Kernel Library (Intel® MKL) is a highly optimized, extensively threaded, and thread-safe library of mathematical functions for engineering, scientific, and financial applications that...
View ArticleWhat performance I should expect from following code
Consider following two part of the codes:/* Perform LU factorization and store in DSS_handle */for(k = 0; k < N; k++){gettimeofday(&stTime, NULL);//DSS solver optionsMKL_INT solOpt =...
View ArticleClik here to view.
Segmentation fault trying to use mkl_dcscmm from a Matlab mex file
I'm trying to use mkl_dcscmm from inside a Matlab mex file (on 64-bit Linux, have tried Matlab versions 2009b and 2012a), but having trouble avoiding segfaults. I've tried static and dynamic linking to...
View ArticleMKL-DSS-DSS-Error, Out of memory
Hi all,I am using MKL DSS to solve a sequence of systems of equations in a finite element calculation.The structure of the systems is always the same, and so, DSS_CREATE, DSS_DEFINE_STRUCTURE...
View ArticleIndexing an array of size between 2^31 and 2^32-1 with LP64?
I am working with spblas. My matrix dimension is about 300k, nnz is between 2^31 and 2^32-1. To keep the memory consumption as small as possible. I would like to use 32bit unsigned integer to index my...
View ArticleHow to get LU sparse matrix decomposition?
Hello everyone,how could I get LU sparse matrix decomposition after PARDISO phase 22? May be, MKL has special routines to do it?Thanks.
View ArticleSparse FEAST
Hi,I ran into an issue with the sparse FEAST implementation of "dfeast_scsrev"& "dfeast_scsrgv".The actual matrix I need eigen-decomposed is of size 4096 x 4096. Here is a test case that presents...
View Article-mkl=parallel
I am using the Intel Fortran Composer with the MKL library. I am calling pardiso to solve a sparse system of equations. Everything was working well in OS X using the -mkl=parallel option with the...
View Article