misunderstanding of dss_create Direct Sparse Solver (DSS) Interface Routines
Dear sir or madam:I need to use dss interface in c++ program.From Intel® Math Kernel Library Reference Manual MKL 11.1, Page 2973:The default value of the parameter opt is MKL_DSS_MSG_LVL_WARNING +...
View ArticleIs there a problem with hemm() ?
If I multiply two complex diagonal 2x2 matrices, with chemm(), I don't get a diagonal matrix ! Here is the code for C=A.B Matrix A is stored as an upper triangle hermitian matrix. #include...
View ArticleCompiling R 3.0.2 with MKL
I'm not sure if R compiled correctly with MKL:When I do ldd on R I don't see any reference to MKLbenjamin@benjamin-Lenovo-IdeaPad-Y510P:/opt/intel/mkl/tools$ ldd /usr/lib/R/bin/exec/R...
View ArticleFunction to solve the linear "less-equality"-constrained least squares (LSE)
Hi, I know there is ?gglse that solves the linear equality-constrained least squares problem using a generalized RQ factorization minimize ||c - A*x||2 subject to B*x = d However I would like to solve...
View ArticleSymmetric sparse matrix - dense matrix multiplication
I need to multiply a symmetric sparse matrix A with a dense matrix X (Y = A*X) using multi-thread/core. The matrices I'm using are the adjacency matrix of graphs, with large number of nodes (up to 2...
View ArticleLAPACKE_sgels example
I am trying to follow example available on http://www.nag.com/lapack-ex/node55.htmlwhere they have: DGGLSE Example Program Data 6 4 2 :Values of M, N and P -0.57 -1.28 -0.39 0.25 -1.93 1.08 -0.31...
View ArticleClik here to view.
Problems with DSS routine
I am beginner of MKL. I just want to use the DSS routine to solve a large sparse matrix equation. I revised the DSS example according to the manual. But it could not work with the error message...
View ArticleHow do we select the column of a matrix?
Hello!I am new to MKL and to C in general. Lets say that we have two matrices A, B of size (n,k). We would like for j =1 until k to subtrack the j column of A from B. In Matlab language the code is:...
View ArticleI can't run MKL with Visual studio.
Hello!I installed on my computer Intel® C++ Studio XE for Windows 2013 SP1 and Microsoft Visual Studio Pro 2012. Then I followed the instructions that I watched on this video...
View ArticleClik here to view.
PARDISO segmentation fault
idbc wrote after 80% of LL' factorization:Program received signal SIGSEGV mkl_blas_mc_sgem2vu_odd () in /mnt/storage/opt/intel/composer_xe_2013_sp1.0.080/mkl/lib/intel64/libmkl_mc.soin the attachment...
View ArticleUse ?gbmv with a 3d array.
Let assume the we have a 3d array A and a matrix B of size 2x5. int A[3][2][2]= { {{1,2},{3,4} }, {{5,6},{7,8}}, {{9,10},{11,12}} }; I want to program the following:cs =0.0;for j=1 until 3Multiply the...
View ArticleExtract row/column of a matrix and subvector form a vector.
Hi! I would like to ask if there are routines for the following issues:1) Extract a row/column of a matrix.2) If theta is a vector extract a subvector. For example extract the elements 4 until 8 of...
View Article3D convolution
I am using MKL's VSL multi-dimensional convolution routines to perform a 3D convolution. As a test, I am using the input function as a 3D gaussian as follows:f(x,y,z) = exp(- (x-2)**2) exp(-(y+2)**2)...
View ArticleSparse Matrix storage to compressed sparse row format
Hi :-)I am new here and I have a question about using "PARDISO - Parallel Direct Sparse Solver Interface". As you may know in order to use this subroutine for solving a set sparse linear equations as...
View Articleselling R / Python pre-built with Intel MKL
Just a thought: instead of Intel selling MKL and then letting the clients build/extend R or Python with Intel MKL, it could be a lot easier if Intel could sell corresponding versions of R or Python...
View ArticleMKL LAPACK SBGVX example
SBGVX computes a few eigenvalues/eigenvectors of a real symmetric matrix so it should be very useful for FEA where we want only a few eigenvalues/vectors of a very large matrix. So, I implemented a...
View ArticleLinking Scalapack with Microsoft MPI in 32 bit
Hi Everyone, I have a project in Visual Studio 2010 - C++ that I am trying to link with MKL Scalapack and Microsoft MPI HPC 2012. I can create a 64bit binary, but I am stuck creating a working 32bit...
View ArticleECCN for MKL
Hi all,Could you please tell me what is the ECCN for MKL ? Is it the same for all versions ? Are there specific libraries inside MKL with specific ECCN (I am thinking about crypto related libs if such...
View ArticleEigenvalue decomposition with CSR sparse matrix
Hi, at the moment I use dsyevd to compute the eigenvalues and eigenvectors of a large matrix A (n = 22000). This takes about half an hour. I know that they are a lot of zeros in matrix A (90% are...
View ArticleClik here to view.
mkl_dcsrmv slower than openMP implementation
Hi,I'm trying to find the fastest way to do a multithreaded sparse matrix-vector multiply. I've written some benchmarking code to form a large random sparse matrix in CSR format, and then time 3...
View Article