Symmetric Sparse Inversion
Hi,I am trying to find an inverse of a sparse matrix and then multiply it with another sparse matrix: A*inv(B). But I could not find any routine for making inverese of a sparse matrix. The only...
View ArticleWebinar recording: Powered by MKL - Accelerating NumPy and SciPy
The webinar recording is now available online: http://software.intel.com/en-us/articles/powered-by-mkl-accelerating-numpy-and-scipy-performance-with-intel-mkl-python NumPy/SciPy are scientific...
View ArticleUnable to run a mkl_malloc function from a package in Go
Hi, I am trying to use mkl_malloc() in Golang.I have created a package which wraps mkl_malloc(). I called the wrapper from other Go file.The problem comes when I am calling the wrapper. The package is...
View Articlemkl_ddnscsr problems
Hello everyoneI am trying to use the ddnscsr to build the arrays a, ia, ja which will then be used for the pardiso or other routines which need CSR formatting. Basically to make sure I understand...
View ArticleSSOR preconditioner using RCI
Hey Guys,I have a problem with impelementing SSOR precondtioner using RCI and hope you can help me with that.I have written a cg solver using RCI for a symmetric matrix that is stored in CSR format 3...
View ArticleHow to evaluate the incomplete gamma value with mkl?
Now i want to evaluate the incomplete gamma value.But i can not find the corresponding function in mkl. Is this function lacking in mkl?
View Articledgemm: slow performance
I am using dgemm function from Math Kernel Library in my FORTRAN source code to do the following matrix calculation: X = A - transpose(B) *CA is 400*400 dense matrix, B is 10000*400 sparse matrix and...
View ArticleScalapack BUG in PZUNGQR (Again)
Hi,Previously, we reported a possible Scalapack bug in the PZUNGQR function http://software.intel.com/en-us/forums/topic/473803That issue has still not been resolved, but it was stated that...
View ArticleInvitation to evaluate Intel® MKL Sparse Matrix Vector Multiply Format...
We are seeking interested parties to evaluate Intel® MKL SpMV Format Prototype Package for Intel® Xeon Phi™ coprocessors. Sparse Matrix Vector Multiply (SpMV) is an important operation in many...
View ArticleInvitation to evaluate Intel® MKL Sparse Matrix Vector Multiply Format...
We are seeking interested parties to evaluate Intel® MKL SpMV Format Prototype Package for Intel® Xeon Phi™ coprocessors. Sparse Matrix Vector Multiply (SpMV) is an important operation in many...
View ArticleSparse Eigensolver Question
I'm looking for a sparse eigenvalue solver that returns just the number of eigenvalues starting from the lowest one that I specify rather than all of them between emin and emax. For example, I would...
View Articlesegmentation fault
Hello everyone,I am new usint MKL and get an "segmentation fault" error with my first MKL-programm. Can anyone help me?My program computes:A = alpha * B + beta * CIt should be very easy. Segmentation...
View Articleoffload error: cannot load library to the device 0 (error code 20)
Hi, I'm trying to get a peice of compiler assisted offoad code working. The code uses MKL - it's just a test code though. At the moment, the code generates some random numbers on the host CPU,...
View ArticleLinking both ipp and MKL causes problems
I get linking issues of the form ippsmt.lib(pscopye9as_l9.obj) : error LNK2005: l9_ownsSet_32s_E9 already defined in mkl_core.lib(pscopye9as_l9_20121126.obj) with Intel C 14.0.1 on Windows x64. I am...
View ArticleSummary statistics and padding
When using matrices in MKL, I'm careful to align my columns as recommended by Intel for maximum throughput (this seems to be particularly important on the Intel Xeon Phi). MKL's BLAS interface makes it...
View ArticleScalapack linear solver (memory problem)
Hi all,I am trying to use Scalapack in order to solve a distributed linear system. The C++ source code is reported in attachment. The code compiles and runs without problems and it gives the right...
View ArticleReduce overhead for repeated offload convolutions on Xeon Phi
I have a series of arrays that I need to convolve. They are not all the same size, but their size does have an upper bound which is ~50. The offload looks something like this (simplified): #pragma...
View ArticleFFT Failure Within SideFX Houdini
Hi guys,This is a long shot but I'm hoping you can help, as I'm totally out of ideas and don't know where to go with this. I'm trying to run a plugin inside SideFX's Houdini that uses MKL's FFT. I had...
View ArticleProblem linking MKL with pthread and gomp
A minimal test program:#include "mkl_lapacke.h"int main(int argc, char** argv) { float a = 0.5f; LAPACKE_slasrt('I', 1, &a); return 0; }I used the "Library Link Line Advisor" to get the...
View ArticleLAPACKE_sgels example
I am trying to follow example available on http://www.nag.com/lapack-ex/node55.htmlwhere they have: DGGLSE Example Program Data 6 4 2 :Values of M, N and P -0.57 -1.28 -0.39 0.25 -1.93 1.08 -0.31...
View Article