Alter elements in diagonal of D-matrix during the LDL' factorization
Hi! I want to solve a system with mtype = -2, i.e. PAP' = LDL', using pardiso MKL. I have understood that the function pardiso_getdiag provides the diagonal of the initial and factorized system....
View Articleomp threading of sparse matrix operation does not yield expected seed up....
Hi all,I have implemented a kind of a gauss seidel algorithm for a sparse upper triangular coefficient matrix, where multiplication of a vector with a the sparse matrix (forward propagation) is a core...
View ArticleSparseQR v.2019 inital release problem
I'am trying to execute the SparseQR methods, but get an unhandled exception when staring the either mkl_sparse_qr_reorder() or mkl_sparse_d_qr()In both cases I construct the sparse_matrix_t object with...
View ArticleReturn mkl_sparse_spmm Results?
Hi all,I am trying to use mkl_sparse_spmm in a function, to calculate C = A*B, where A, B, C are in CSR format.I would like to return the information of C, i.e. rowIndex_C, column_C and values_C...
View ArticleClik here to view.
mp_linpack results problem
Good day!I have a weird problem with the Xeon Phi x200 (KNL) performance on the mp_linpack test (HPLinpack 2.1) from the set of MKL benchmarks of Intel Parallel Studio XE 2017.If I run the test on 1...
View ArticlePardiso Threadripper 2990wx versus Ryzen 1700
I have the same multi-physics finite element code generating a matrix. An old machine with a Ryzen 1700 (8 core) is faster than a threadripper 2990wx (32 core). Windows 10, intel64, mkl_rt.lib, and the...
View ArticleClik here to view.
MKL ScaLAPACK PDGEQRF routine on Stampede2 -- Hangs on 512 nodes with...
Hi, my name is Edward Hutter and I am a graduate student at the University of Illinois. I have been benchmarking my code against MKL's ScaLAPACK routine for Householder QR factorization, PDGEQRF, and...
View Articleoperate on diagonal of a sparse matrix in csr format
hi,I want to efficiently operate only on the diagonal of a large non-singular sparse_status_t matrix (in csr format).More precisely, if 'A' is my matrix then I want to do:A - alpha*Iwhere 'alpha' is a...
View Articlesegment error on mkl_sparse_?_ev when the matrix has degenerate eigenvalues
Hi everyone,I am trying to calculate eigenvalues of some real symmetric sparse matrices. What I notice is that the program crashes occasionally. I found a minimum working example:#include...
View ArticleMKL-DNN primitive examples
Hello people ,I'm trying to write a personal DNN primitive ( relu or maxpool ) and i need to compare performance with mkl-dnn primitives. the problem is that in the examples provided in github i could...
View ArticleMKL Windows 10 Installation Fails without error
Dear community, I am trying to install the latest standalone Intel MKL for windows (w_mkl_2019.0.117.exe). However when double clicking the file, I get asked for admin password and then... nothing...
View ArticleMKL > 2018.3 doesn't work with NVIDIA nvblas
Platform is Ubuntu 16.04.5. GPU NVIDIA1070 cuda-9.2 python3.5I use python with nvblas support by compiling numpy against intel MKL. It works with 2018.1 and 2018.2 but GPU is not used for MKL 2018.3...
View ArticleClik here to view.
Does MKL support sparse matrices in LAPACK?
From the LAPACK's FAQ, we know it does not support sparse matrices:Just to take a chance, does intel MKL break this restriction by introducing sparse matrices in LAPACK? If not, is there any other...
View ArticleFFT not working for N >= 2048
I'm attempting to compute an FFT on sizes N>=2048. What I can't explain is why the FFT output looks great for N<2048 but then goes to garbage when N>=2048. I've poured over other forum posts...
View ArticleDocumentation for MKL_DIRECT_JIT
Is there any documentation or explanation for MKL_DIRECT_JIT apart from a one-line mention in the release notes?
View ArticleZeros in the result of mkl_sparse_spmm
Hello,Per https://software.intel.com/en-us/mkl-developer-reference-c-sparse-blas-c... a matrix stored in CSR storage format should contain only non-zero values. We've ran into a case where the result...
View Articleifort/MKL vs Julia/OpenBLAS: ifort/MKL using less cores
I am new to Fortran and I am comparing ifort's speed to Julia. I would greatly appreciate any help or insight into my issues and taking the time to answer some of my questions. When running what...
View ArticleIntel® MKL version 2018 Update 4 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 ArticleClik here to view.
mkl_sparse_d_gv - Internal Error
Hi there,I am running some generalized eigenvalue buckling problems (Kg.v = (1/λ).Ke.v) in MKL and in most cases getting an internal error. I'm using this C# Library and coding in VB.NET. I have ruled...
View Articlemkl_sparse_?_syprd implementation in python
Hallo,I would like to use mkl_sparse_z_syprd in pythonhttps://software.intel.com/en-us/mkl-developer-reference-c-mkl-sparse-syprd In comparison with the older functions which require indices and data...
View Article