Intel® Math Kernel Library 11.1 Update 3 is now available
Intel® Math Kernel Library (Intel® MKL) is a highly optimized, extensively threaded, and thread-safe library of mathematical functions for engineering, scientific, and financial applications that...
View Articlerandom numbers with vdrnguniform
Hello,I want to generate random number with the SFMT19937 generator. Everything works fine as long as I store the random numbers to a vector, e.g.... brng=VSL_BRNG_SFMT19937...
View ArticleClik here to view.
Help! Negative output of IPARM (18) : Number nonzeros in factors.
Hi,Recently I solved a sparse non-symmetric real matrix equation with 2,785,003 unknowns,The question is that I found Negative output of IPARM (18), which is the number of nonzeros in factors... ,...
View Articleprojection operation in MKL
Dear all,For a projection operation, I should implementF = ( I - Mf*Vfr*Vfr^T ) * FwhereMf is a sparse matrix of size sz_fXsz_fVfr is a vector of sz_fX1F is a block of vectors of size sz_fXP, and...
View ArticleParameters of the objective function for the nonlinear fitting routine
I want to conduct a nonlinear least square fitting using MKL. However, I don't know how I can introduce some additional parameters for the objective function. Does anybody have some ideas?This is the...
View ArticleMKL 11.0 performance problem on Opteron
Hello,I am observing a weird performance behavior when using MKL 11.0 on an AMD Opteron 6220 processor. I wrote a small benchmark that measures the Gflop rate of the _GEMM operation in S, D, C and Z...
View Articlestrnlsq related routines for nonlinear fitting
I want to solve a nonlinear least square fitting problem to estimate a and b. g(x)=a(x^2)+bx. I have two vectors working as the training data. fx1={1.0, 2.0, 3.0, 4.0}; fx2={2.0, 4.0, 6.0, 8.0}. I...
View ArticleGemm BLAS routine
Has anyone found an issue on using a BLAS/LAPACK routine after putting the computer in a power save mode? I have found that the SGEMM routine does not work after the computer awakens from sleep or...
View Article?trnlsq related routines for nonlinear fitting
I want to solve a nonlinear least square fitting problem to estimate a and b. g(x)=a(x^2)+bx. I have two vectors working as the training data. fx1={1.0, 2.0, 3.0, 4.0}; fx2={2.0, 4.0, 6.0, 8.0}. I...
View ArticleForum poll: Intel MKL and threading
Intel MKL users,We would like to hear from you how you are using Intel MKL with threading. Do you use the parallel or sequential MKL? How do your multithreaded applications use MKL? We would appreciate...
View ArticleClik here to view.
Low Accuracy with CG Solver ?
I am just wondering why , the CG Solver is so inaccurate and also in some cases the solver doesnt covert. So how to make the solution of x more accurate and how to set the maximum number of itereations...
View ArticleCorrelation
Hi all, I have a lot of trouble understanding how the correlation/convolution functions are used. I have looked at the Information from the Manual but the description - involving finite functions - is...
View Articlepdlacpy does not work with matrices distributed differently.
I am using pdlacpy. Documentation says this function copies all or part of a distributed matrix (A) to other distributed matrix (B). It also says no communication is performed (e.g. the function...
View ArticlePerformance difference bewteen i7 and xeon CPU for intel mkl
Hi~ I think i might meet with some wired problem with Intel MKL.I install visual studio 2013 and intel xe composer 2013 sp1 on two machines, one is pc with windows 8.1 equped with intel core i5-3470...
View ArticleIs it possible to update Schur complement quickly with MKL/Pardiso?
I have a matrixA = [A11 A12;A21 A22];and I'd like to compute the Schur complement of A11 in A using pardiso: S = A22 - A21 * A11^-1 * A12. This seems well supported by the PARDISO-SCHUR routine...
View ArticleIntel MKL 11.1 java example build issue
Hello,I am trying to build Intel MKL 11.1 java example as discripted here (https://software.intel.com/en-us/articles/performance-tools-for-software...). I got the following error while running nmake....
View ArticleIs there any way to use VML on 2-dimensional array?
Suppose I have a 2-dimensional array a[m][n],and I want to get the result array of one row of a[m][n] adding another row of a[m][n].Is there any way to use VML in this case？I‘ve tried to directly use...
View Articleerror for R calling MKL with multithreading
Specifically, when using mkl functions (such as lapacke_dgesvd) in C++ with multi-threading, it works fine. But with the same function and R binding to call it under multi-threading, the result of SVD...
View ArticleMKL_malloc interface missing from Fortran module
Hi,It looks like the interface for MKL_malloc is missing from the mkl_service Fortran 90 module, at least in the versions I have. Is it supposed to be in there? The documentation seems to indicate so,...
View ArticleMKL ILU for BSR
I am the technical manager in a reservoir simulation project and working with various self made and existing math libraries for my speed comparisons. I use ILU (Incomplete LU) factorization on a...
View Article