MKL BLAS95 and LAPACK95 interface build
Hi, I'm trying to run the following MKL interface make files in RHEL MRG (3.10.33-rt32.34.el6rt.x86_64) in order to interface to MKL using the BLAS and LAPACK interface (using the mkl_ prefixed header...
View Articletranspose matrix stored in CSR format using MKL
All,I have some legacy parallel code that uses CSR format to store a very large, sparse matrix. In making some additions to the code, I have the need to transpose the matrix, storing the transpose in...
View Articlepardiso_getdiag
When I do acall pardiso_getdiag( pt, df, da, mnum, error)the arrays DF and DA are filled with 'factored' and 'original' diagonal pivots.The pivots are in the same 'order' of the 'original' matrix or...
View ArticleIntel® MKL Cookbook Recipes
Intel MKL Users,We would like to Introduce a new feature Intel® MKL Cookbook, an online Document with recipes for assembling Intel MKL routines for solving complex problems.Please give us your valuable...
View ArticleError in the documentation of feast
Hi,I would like to report an error in the documentation of the reference manual for the extended eigensolver (feast). In the rci interface description, the parameter ijob is detailed. When ijob = 30 or...
View ArticleEigenvalue Solver Error (dfeast_scsrgv )
I'm trying to solve eigenvalue problem for large sparse matrices using dfeast_scsrgv function. The function works fine for small problems (ex: 8*8 sparse matrix) but it gives...
View Articlemkl_?csradd
Dear alllooks as if I need to use mkl_?csradd (C++).However, my csr matrix is zero based; mkl_?csradd does only support one based sparse matrices.I would appreciate any suggestions for a pragmatic...
View ArticleFGMRES
All,I just started learning to use the mkl and fgmres in general and I'm having troubles implementing/understanding the fgmres. I've looked at the example that came in the library and I've been able to...
View Articlelibiomp5md.dll location (release build)
Hello, where would I find the dll without tracing code? Currently I use one from MKL 11.1.2 (64-bit, version 5.0.2013.1126, file size 1043kB, modified 2014-01-31 12:23) and the profiler shows this...
View ArticleClik here to view.
Error in simmetric sparse matrix analysis step
I am a new user of a MKL, so question will be silly. I am trying to use pardiso_64 function to solve a system with symmetric sparse matrix. IDE is MS VS 2013, language is C++. Matrix size in test case...
View ArticleStack Overflow Error in Eigenvalue Solver: dfeast_scsrgv
Hi,I'm trying to solve eigenvalue problem for large sparse matrices using dfeast_scsrgv function. The function works fine for small problems (ex: 8*8 sparse matrix) but it gives...
View ArticlePreconditioners for banded matrix (diagonal storage format)?
All,Is there a built-in routine for computing a preconditioner for a sparse matrix in diagonal storage format or what other alternatives exist? I plan on using it with the fmgres routine.ThanksM.
View ArticleWebinar announcement: A Tour of the Sparse Linear Algebra Functionality in...
Sparse matrix algorithms are encountered in a broad range of important scientific computing applications. Intel MKL offers a powerful set of functions that can be used to build a complete solution to...
View ArticleMKL FFT and FFTW
Does the MKL Library use FFTW internally?If there are other facilities as well, are there anmy limitations in the radixes?Thanks
View ArticleSpeed of solvers with CSR format
I have been using MKL in composer 2013. Especially, pardiso and preconditioned conjugate gradient solver with CSR format to solve symmetric matrices I wonder that using full element CSR format is...
View ArticleLinking problem
Hello, I am trying to solve a linear system of equations using LAPACK ROUTINES and I have tried to link but some of the codes seem to work whiles others do not. I get the error when use the CALL GESV....
View ArticlePardiso L and U factors
I have been trying to figure out a way to extract the L and U factors of a matrix using pardiso. I read somewhere that this is not possible, but I just wanted to check with people here to see if there...
View ArticleMemory doubled when calling zgetrf in Csharp,
Hi,I'm solving a large complex dense matrix (10k*10k) in Csharp by calling zgetrf first to LU-decompse the matrix. If I understand correctly, zgetrf will do a in-place LU-decomposition. So there should...
View ArticleVSL Summary Statistics VSL_SS_SUM error
Pls let me know why the code bellow shows invalid results(sum) when NR is greater than 1999.#define NR 2000 // MEAN=1 but invalid SUM ! //#define NR 1999 // MEAN=1 and SUM==1999 of course #define...
View ArticleSegfault in the dtpmqrt routine
Hello,While testing out MKL Lapack's dtpqrt and dtpmqrt routines, I've stumbled across a weird segfault. I replicated the error in this example (I should mention that I use Eigen just to make my life...
View Article