Clik 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 ArticleHow to compile the sample code in minGW complier?
Hi all,I am using minGW as the compiler and our code has a lot of FFT. In the past we use fftw and we want to try Intel's fft to improve it. I download the trial version of the MKL and install it on my...
View ArticlePARDISO gives wrong answers when given too many rhs's
When I call PARDISO with too many rhs's I get wrong answers.I first call PARDISO for phase 11 and phase 22. I keep the factored matrix in memory, and, subsequently leave the subroutine that calls...
View ArticleCompilation Interface issue with dfeast_scsrev
Hi,I recently upgrades to XE2013 to get access to the extended eigensolvers in the MKL11 libraries.I am trying to use dfeast_scsrev but am getting some compilation errors.Here is my code: SUBROUTINE...
View Articleadvice for getting rid of temporary creation and fortran95 interface for dsecnd
I am using intel fortran compiler and intel mkl for a performance check. I am passing some array sections to Fortran 77 interface with calls like<code>call dgemm( transa,transb,sz_s,P,P,& a,...
View Articlemixing fortran77 and fortran95 interfaces
Dear all,I programmed my code with the fortran 77 interfaces. However I was wondering if that is possible to mix fortran77 and fortran95 routines. For instance, looking at the sparse blas level 2...
View Articlefortran95 interface syevr return -1001 as info
Dear all,As the subject line suggests, fortran95 interface of dsyevr, namely, syevr returns -1001 as info parameter. I could not understand what this means, could you please help me with this...
View Article2D FFT
Hello,I have been assigned the task of converting a matlab script to C++, and am currently working on the FFT part.Given this input:octave:67> r57 r57 = 0.00000 0.20000 0.30000 0.30000...
View Articleinconsistent results from mkl_dcscmv and mkl_dcsrmv
Hello,I tried to calculate the (sparse) matrix-vector product using mkl_dcscmv and mkl_dcsrmv. However, sometimes they gave me different result. Here is the example I was using:Let M : = [ 1 1 0 0...
View Articlehow to tell mkl the underlying microarchitecture
Hi: I read the gcc document as well as intel compiler document, they both said that the default behavior will not detect the underlying microarchitecture in linux x86-64 (default is -xsse2). As a...
View Article