cgemm3m, cgemm_compact AND cgemm give poor results for small problem 24*64
Hi all,I using sequential API and direct call to multilply matrices. C = 1*conj(A')*AA is 64*24 and C is 24*24 both are complex matrix (complex8).I have arrays of matrices: A_ARR (filled with random...
View Articlemkl_jit_create_cgemm on multi cores
Hi all,I want to use mkl_jit_create_cgemm on my setup where each thread is pinned to single core.In each thread I'll do the cgemm with the created Jitter.Do I need to create jit kernel specific for...
View ArticleSome subroutines of lapack in mkl gives Segmentation fault for matrix...
OS: Ubuntu 18.04.2 server, kernel 4.15.0-50-generic. MKL provided by intel parallel studio xe 2019.3.The C code reproducing the error#include <stdlib.h> #include <stdio.h> #include...
View ArticleClik here to view.
MKL FFT Error in Example Code in Linux
Hello!I was trying to learn about the MKL FFT libraries and wanted to check how the example programs ran, but when I check the output files, it says there is an error, status = 2 while trying to create...
View ArticleClik here to view.
MKL FFT Fortran in Test Code - Transformed Values are not correct
Hello!I was trying to learn about the MKL's FFT Function and wrote a small 1D program to show the forward transformed values. From my understanding, the arguments for the forward function - mentioned...
View ArticleCannot use MKL F95 eigensolver routines
I am using Intel Parallel Studio XE with Visual Studio. I would like to use the heev()/heevr() routines. I have written a simple code to test it, and I have enabled the "Use MKL libraries" option in...
View ArticleDifferent LU factorization result between libmkl_intel_ilp64.a and...
I am a novice in the field of Intel Math Kernel Library. However, when I tried to use the LU factorization function from LAPACK by compiling different libraries, I got different results as the...
View ArticleDifferent BLAS performance between CMake dynamic and manuly static link
I tried to compile my program by using CMake at first and my CMakeLists file is as the following. cmake_minimum_required(VERSION 3.11)project(LMMNET LANGUAGES...
View ArticleHow to get Cholesky diagonal?
HelloI need to extract diagonal from cholesky LDLt sparse factorization. Can I do this with Intel MKL?I can use pardiso to factor and solve sparse linear system. But cholesky factor matrix itself is...
View ArticleBug in mkl_sparse_d_mm
Hi All,I think that the result of "mkl_sparse_d_mm" is incorrect as evident for the following code. The code basically tries to multiply the Identity matrix stored in CSR format with a dense matrix and...
View ArticlePossible issue in MKL LAPACKE interface when using LAPACK_ROW_MAJOR
MKL provided by Intel parallel studio xe 2019.3.The minimal reproducing code is below:#include <stdlib.h> #include <stdio.h> #include "mkl_lapacke.h" #define N 2 int main(){ double m[N];...
View ArticleDiscontinued Xeon Phi support
Hello,Whether Intel will support Xeon in future versions of compilers(+MKL)?Malik.
View ArticleQuestions about iterative sparse solvers
Hello,I would like to use ISS from MKL for solving linear systems of equations, originating from discretization of partial differential equations used in CFD.I found out that CG and FGMRES are...
View ArticleMKL Packed Storage scheme for matrix with many zoroes
I have a very large matrix but many of its members are zero but it does not have an specific shape like being banded or etc.Is there any possiblity that I can reduce memory usage for storing this...
View ArticleMaintaining 15+ year old software package, MKL issue support
Hello, I've taken over an older software package my company wants to add some functionality to. We don't have funds or time to update all aspects of it, so as is, with new modules is the goal. I have...
View ArticleSparse Complex Matrix Multiplication
Hello, Let's say A and B are sparse matrices with complex values, which are stored in sparse handles (CSR format). I want to calculate D = A x conjugate(B), where conjugate is the complex conjugate of...
View ArticleError in one LAPACK example
Hi, I've located an error in one of your FORTRAN 77 interface to the LAPACK Single Value Decomposition routines (DGESVD). This is located in...
View ArticleHow to do an element-wise multiplication between two 3D-arrays?
Hi, everyone.I am trying to do some element-wise multiplication between two 3D arrays/matrices and then calculate the sum of all the elements returned. Something like this:res = 0 do k = kstart,kend do...
View ArticleClik here to view.
MKL_SPARSE_?_QR invalid input error (Fortran)
Hello,I am experimenting the Sparse QR library to solve a linear system of equations. When I run the attached file, I get the following output:0003The last output "3" indicates...
View ArticlePardiso 32-bit vs 64-bit performance
Using Pardiso in an iterative routine, where the Pardiso solver is called many times, I noticed the 64-bit version is significantly slower than the 32-bit version. I'm not sure yet if it depends by...
View Article