Number of non-zeros less than nRows
Hello. I have two very related questions that I hope you can answer. I am writing code that uses the MKL to get the LDL' decomposition of a matrix that is real and symmetric but, in general,...
View ArticleUse Which Routine for Sparse Matrix-Sparse Vector Multiplication
I have a huge matrix with very few non-zero elements. Some of the columns and rows may also be completely zero. This matrix should be multiplied by a very long vector which has only few non-zero...
View ArticleQ: dtrnlsp_init for initializing the solver of a nonlinear least squares...
Hello All,We're replacing our current solver for a non-linear least squares problem without constraints with the MKL version.First tried implementing it in a simple problem: F(x) : n = 6 → m =...
View Articledtrnlspbc_get() returns 0xbaadf00d in "iter" argument
I'm using MKL 11 (tried on 11.04 and 11.0.5) and I'm getting a bad return value upon the return of this function. The "iter" argument gets filled with the magic number value 0xbaadf00d, which means...
View ArticleCan I pass a subset of a matrix into another function in MKL?
I am trying to optimize a lot of matrix calculations in MKL that requires me to allocate large blocks of memory using something like :double* test_matrix = (double*)mkl_malloc(n * sizeof(double),...
View ArticleClik here to view.
dcopy problem with very large arrays
I have run into a problem when using dcopy with very large arrays. Specifically, I am using dcopy to initialize a very large array to zero before doing other operations, and the dcopy is not setting...
View ArticleScalapack p*gemr2d return codes
Hi,I'm using the MKL Scalapack functions P*RGEM2D to distribute matrices between different contexts. Most MKL Scalapack functions take an 'info' argument to return various error information from the...
View ArticleMKL thread priorities
Hello All,In version 11 of the Linux-based MKL, is it possible to either:1) specify the pthread attributes to be used when creating any MKL threads,or2) have the MKL threads perform a callback to a...
View ArticleFor what reasons can mkl_free() fail?
I am working with Intel's MKL in C++ and I get a lot of runtime errors with themkl_free(some_matrix); call. Usually, it just hangs when it reaches the matrix. What are the reasons that can cause it to...
View Articlehow to make eigenvectors ( obtained from _geev ) orthorgonal?
I am looking for the eigenvectors & eigenvalues of a non-symmetric, complex-valued matrix. The matrix is a Hermitian matrix with some diagonal imaginary terms added, which is then non-hermitian....
View ArticleDoes mkl_dcsrmm require row-major matrices?
mkl_dcsrmm multiplies a CSR sparse matrix by a dense matrix, storing the result in a dense matrix.From what I can tell, this function expects the B and C dense matrices to be in row-major order, as...
View ArticleFFT 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 Article