LAPACKE_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 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 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 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 Articleblacs for sgi mpt
I need to compile a software package which requires blacs, however since about composer_xe_2013.5.192 blacs seems to have disappeared from MKL (except blacs for intelmpi and MIC). Where do I get the...
View ArticleScalapack on Mac OS X
What are the chances of seeing scalapack in mkl for Mac OS X in near future?
View ArticleProblem with robust estimation of a covariance matrix.
Hi I have the matrix "x" and I want to compute the covariance matrix. The i column of the matrix stores the observations of the i variable.The matrix is 0.8147 0.9058 0.1270 0.9134...
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 Articleproblem in computation of 2D FFT using two seperate real and imaginary arrays
Hi,i want to compute FFT of a complex 2D array. So first i have tried it in 2 ways.1) Directly using complex 2D array(real and imaginary interleaved) .2) 2 seperate arrays(where real and imaginary are...
View Articlevsrnguniform
Dear all,Can vsrnguniform generate random multidimensional arrays or it just works for 1D arrays? I guess that one solution could be to generate a 1D array and after reshape it. However, I would like...
View Article