Clik here to view.
Different results getrf/getrs, dss and intel pardiso
Dear all,I have a small nonsymmetric linear system that is represented by a matrix in csr format (file fort.106). The task is to solve the system. To this end, I applied three different approaches. At...
View ArticleIntel MKL ERROR: Parameter 4 was incorrect on entry to DSTEIN2
Hi, Some users are seeing the following error messages when running VASP, built against Intel 2018.0.3.Intel MKL ERROR: Parameter 4 was incorrect on entry to DSTEIN2.Intel MKL ERROR: Parameter 6 was...
View ArticlePerformance characteristics of cblas_gemm_s16s16s32
Hi,I'm interested to get more details on the performance characteristics of the function cblas_gemm_s16s16s32. In my application, the performance gain over cblas_sgemm is lower than I would hope.Here...
View ArticleRandom Number Generator
Hi,I have developed a function to generate a vector of random numbers. However, there are two problems:1. When the number of random numbers to be generated is too large (e.g. 50000^2), it will output:...
View ArticlePackage script bug: link_install.sh and tr -s [:blank:]
Having installed the package intel-comp-l-all-vars-19.0.1-144 via apt on Ubuntu, I've run into a problem with the script at /opt/intel/compilers_and_libraries_2019.1.144/linux/bin/link_install.sh which...
View Articlemkl/blas routine for C=AA'B
Hi there,is there any mkl/blas function which performs the operation C=AA'B in on go. Currently I use an intermediate array T and dgemm: T=A'B;C=AT'. I am wondering whether there is a more efficient...
View Articleeigen value and eigen vector
I'm working on quantum structures, so I need to calculate eigenvalues and eigenvalues in a very precise drawing, I want to work with values such as n = 10000 and above. Mkl library that I use it is...
View Articlemkl(dgemm) performance problems on "superlarge" processors
Hi,I was running two subsequent dgemm operations: T=AB and C=A'T with A=(56,000x400,000), B=(400,000x30), T=(56,000x30) and C=B.Conditional on the CPU I measured these wall clock times (for the dgemm...
View ArticleEingenvalue solver dfeast_syev does not find values (info=1)
I am trying to use dfeast_syev to find eigenvectors of 4x4 matrix. I used FEAST, since I'be found that other methods give incorrect values. I call the routine with the following...
View ArticleClik here to view.
DGELSS Issue
I am attempting to use this routine but am getting an exception:Exception thrown at 0x01FA5428 (mkl_core.dll) in Console1.exe: 0xC0000005: Access violation reading location 0x00000000.I wonder if...
View ArticlePARDISO, How to escape the occurrence of zero pivot error
Hi everyone,I have zero pivot error -4, but I know the error occurs in the condftion of my calculation.So I want to have a solution using zero pivot replacement (the value i_parm(10)) and without...
View Articlemkl_set_num_threads() doesn't work
Hello everyone!I am using MxNet library with MKL support. I need my mxnet predictor use only one thread for calculations. For this purpose I used OMP_NUM_THREADS = 1 environment variable. So it works...
View ArticleRare crashes on MKL
I have implemented in C++ an algorithm in image processing using (among other things) fftw wrappers in MKL library (version 2018.3.210)I am working on a x64 machine with Intel Xeon E5-1650 v3 3.5 GHz...
View ArticleQuestion about the mkl_?omatcopy
I can not understand the manual about the following parameters of the function.rows The number of rows in matrix B (the destination matrix).cols The number of columns in matrix B (the...
View ArticleHow to access the number of non-zero elements in sparse_matrix_t?
I am writing a Python wrapper for calling the 'mkl_sparse_spmm' function. In order to export the result of matrix-matrix multiplication to a Python object, I need to know the size of the 'col_idx' or...
View ArticleIntel Pardiso error numerical factorization
Dear Pardiso users,I am using Intel Pardiso to solve a sparse system. For small matrices it works perfectly and very fast. However, with increasing system size I figured out that for some set of...
View ArticleFinding the eigenvalues (diagonalizing) of a block-diagonal matrix
I have to diagonalize a large matrix which take a lot of time. The matrix size 10,000 x 10,000.This matrix is Hamiltonian of a spin system which have some block structure. Is there way to diagonalize...
View ArticleUsing FMA in MKL routines
Hey everyone,I couldn't find any old topics that dealt with this question in detail, so here I am asking it again: is there a way to enable FMA math when using the MKL routines? Here is a sample...
View Articledfeast_sygv -4 error, BUT B IS POSITIVE DEFINITE!!!
Hello all,I am using the extended eigensolver routines specifically dfeast_sygv function, and I get the following error==>INFO code =: -4Intel MKL Extended Eigensolvers Error: Matrix B is not...
View Article