cblas_dasum is not adding up negative numbers
Hello,I am trying to get the sum of a vector using MKL's cblas_dasum function in C using Intel C++ compiler but the negative numbers are not subtracted. I would expect 40.0 + 2.0 - 42.0 == 0.0 but I am...
View ArticleClik here to view.
*** glibc detected *** ./test.out: free(): invalid next size (fast)
Sir,i am trying to multiply 2 block sparse matrices ie. block elemnts along the diagnol. in my case two 4*4 blocks.a_mat is the nonzero array for A. since A is 8*8 a_mat(columns array as per...
View ArticleOutput of 2DFFT is not matching when dimensions of size 128*128, 256*256,512*512
Hi,I have written a program to check the correctness of 2D FFT using two different menthods.First one has input data in complex format(real and imaginary are interlevaed).2nd one has input data in...
View ArticleWebinar Registration: What's new in Intel® MKL 11.2 Beta
Dear Customers,Please register and join us for the Intel® MKL 11.2 Beta Webinar on May 7th 9.00 AM to 10.00 AM PDT. More Information on new features introduced in MKL 11.2 Beta will be presented in the...
View ArticleManipulation of matrices with sparse storage
Hi,recently I got involved with memory optimization of our in house FE program. I decided to use csr format storage with dss solver, which worked fine so far. However in our implementation, we have so...
View ArticleMKL spurns only minimal number of threads for SGEMM
HiI have a Xeon E5-2620 processor and benchmarking with SGEMM. Why does MKL spurn only 6 threads (hardware threads) instead of the expected 12 threads (hardware plus software threads)?The same code on...
View ArticleExample of submatrix multiplication
Hi,Are there any example of using dgemm for multiplying submatrices of two matrices? Thanks,Pawan
View ArticleHow to get a response from Intel?
Hi,I previously posted topic "strnlspbc_solve divide by zero" that illustrates a bug in the MKL. It has been several weeks now and I haven't received any response from Intel. Is there another method...
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 ArticleMKL compile failure with scalapack/pblas
Hi, When I try to compile some c++ code using scalapack and pblas on Windows, I get the compile errors shown below. A simple test case is #include "mkl_pblas.h" #include "mkl_scalapack.h" int main() {...
View ArticleIntel MKL downloading problem
Dear all,I entered this link : https://software.intel.com/en-us/non-commercial-software-development, and made registration in order to download the Intel MKL library for Linux OS.After finishing...
View ArticleFEAST sparse with parallelism error
Intel MKL 11.1 update 3In some cases, using feast with a sparse matrix, dfeast_scsrgv(), the program can't find any eigenvalue. I've run the same problem with same data and parameters in many computers...
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 ArticleIntel® 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 ArticleCode Examples for CG Solvers on Xeon , Xeon Phi ?
Are there any code examples which runs on the CPU and also on the Xeon Phi,we already have the MKL libaray but i wasn't able to find any big example, this code just should read a mtx file and should...
View ArticleConfused about some eps values for ?trnlspbc_init
The documentation is not quite clear about some of the eps array containing the stop criteria for the TR solver, especially eps(5) and eps(6). Page 4164 of the MKL 11.1 documentation.1. eps(5) : The...
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 Article