strnlspbc_solve divide by zero
I'm running MKL 11 Update 3 and am getting a floating point divide by zero exception at _mkl_trs_dmintr_ls_lc() + 0x10a1 from a call to strnlspbc_solve. I'm using LS to fit a Gaussian, and the code...
View Articlecompiling with GCC using ilp64 interface, keep getting seg-fault
Here is my source file(gm.c) in C, multiplying two n by n matrix, A and B, using dgemm() within Inel MKL:#include <stdio.h> #include <stdlib.h>#include "mkl.h"int main(){ int n = 500;...
View ArticleSparse direct solve on Xeon phi?
HiI have gone through the MKL documentation but find it very unclear which MKL routines can actually offload to Xeon Phi. Does the MKL include any sparse direct solver capable of offloading to Xeon...
View ArticleALLOCATE or MKL_MALLOC?
In Fortran, assuming I have a 64bit avx-able machine at my disposal, Is ALLOCATE statement guaranteed to allocate aligned arrays when the compiler is called with the "-align array64byte" flag? Or, do I...
View ArticleFatal error in MPI_Recv: Invalid tag
Hello,I was developing a new software using scaLAPACK/BLACS MKL implementation. In my small testing cluster everything worked ok. However, when I have moved to the 'big' cluster to start computations,...
View ArticleFatal error in MPI_Recv: Invalid tag
Hello,I was developing a new software using scaLAPACK/BLACS MKL implementation. In my small testing cluster everything worked ok. However, when I have moved to the 'big' cluster to start computations,...
View ArticleClik here to view.
FFT forward and backward test
Hi,I am testing the FFT routines in MKL library. Following the example code, I was able to do 1D forward FFT for an array input supplied by a text file. Then a backward FFT is performed and the result...
View ArticleMaths Apps using Microsoft Excel
Hello AllI stumbled across this community, so not sure if this is the correct forum to discuss my app; apologies if it is not.I am a software developer with 10 years professional experience. I am also...
View Articlecblas_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 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 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 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 Article