cluster_sparse_solver cause segmentation fault in mkl 11.3
Hi: My environment: linux64, mpicxx for MVAPICH2 version 2.0b, icpc version 13.1.3 (gcc version 4.7.0 compatibility). In order not to confuse with the mkl library in icpc version 13.1.3, I put the...
View ArticleGetting the content of *pt in Pardiso
Hiis there a way to see the Content of Workspace for pardiso ( moslty call void PT[64] ), It seems that the Permutation Matrix of PAx=Py preserves quite well the sparsity in LU, or LL* ,...
View ArticleInverse of very small matrix
Dear all, I have a piece of code in Fortran90 in which I have to solve both a non-linear (for which I have to invert the Jacobian matrix) and a linear system of equations. When I say very small I mean...
View Articleunexpected ?potrf subroutine failure
AHi,I need to generate a gaussian correlated noise based on a covariance matrix with fortran. I am using spotrf subroutine to do the cholesky decomposition. However, the decomposition is always failed...
View Articleundefined symbol mkl_blas_avx2_cgemm_copyb_ext
I have tried to embed MKL in a java application by embedding the static mkl libraries in a shared lib of mine and then opening that shared lib from the jni. However, I got a link error when I created...
View Articlemkl_ddiamm bug ?
I am trying to multiply two matrices using mkl_ddiamm method:C = A * Bwhere A is diagonal matrix 3x3 and B is general matrix 3x3. No matter what I try, i get as a result no A*B, but B*A. This is my...
View ArticleDifference between C++ code with Blas/Lapack and Matlab
Hi All!I am trying to implement one Linear programming algorithm in C++. For the matrix multiplication, I use Blas and Lapack. However, I find C++ code performs worse than Matlab as the size of problem...
View ArticleExtra precision
I am evaluating MKL 11.3 for Windows. Could you tell me if MKL provides some functionalities like XBLAS (which supports the extra precision or up to twice the working precision)? I notice that MKL has...
View Articlemkl_zcsrcoo faster computation on subsequent calls?
Hi,I have a sparse matrix in coordinate format (row, col, A) and I transform it to CSR to be used for PARDISO. The sparsity pattern never changes (that is, row and col are always the same). Vector A...
View ArticleError in running VASP 5.3.5 with mkl, ifort and mpif90
I did installed VASP executable successfully, only I changed FC=mpif90 (openmpi compiled using Intel compiler) whatever you mentioned in the following...
View ArticleProblem with PARDISO in Windows 10
I have used PARDISO in Fortran successfully in Windows 7. When I moved to Windows 10 my code suddenly did not want to run and it crashes in the call to PARDISO. It seems something changed in Windows...
View ArticleCall dfsNewTask1D from C#
Hey,I'm trying to get the spline functionality to work in C#. I don't seem to know what to put in place for the DFTaskPtr. I can see it's a void pointer in the mkl_df_types.h file, but I don't know how...
View ArticleQuestion with df?EditPPSpline1D in Data Fitting Functions
Hi all,I have a question about parameter bc of function df?EditPPSpline1D. In the documantation ( https://software.intel.com/en-us/node/522222#D69BA4E7-E8BD-4413-A2A9-769D5002F832 ) says :bcconst...
View ArticlePARADISO Compile Error
Hi, I am trying to use PARADISO from the MKL library, however, I cannot compile my program because of the following C2059 syntax error '(' on Line 71 of mkl_paradiso.h Has anyone else had this...
View Articleslow out-of-place matrix transposition using mkl_domatcopy
I created test program to test speed of matrix transposition for various scenarios, in particular:1. in-place transposition of square matrix using mkl_dimatcopy method 2. out-of-place transposition of...
View ArticleXerbla and CBLAS/BLAS
Hello,I have a problem using xerbla and cblas/blas functions. The following code ends up with a segmentation fault:#include "mkl.h" #include <iostream> void XERBLA(const char * Name, const int *...
View ArticleNo PZLAHQR in MKL
Hello,I am using ScaLapapck from MKL.Is there a PZLAHRQ routine in MLK?The reference manual says it has only a real version https://software.intel.com/en-us/node/521530If no, then how can I do a Schur...
View ArticleSparse BLAS - sparse_dense matrix multiplication with zero indexing and...
Hello,I would like to use MKL sparse blas in C++ for computing multithreaded sparse-dense matrix multiplication. I would like to compute C = S * B, with S a sparse matrix in format coo or csc in...
View ArticleMatrix Inversion LAPACKE_zsytri
hello everyone,I'm testing the inversion of a symmetric matrix, in order to do that, I wrote this code:Complex[,] c = (Complex[,])a.Clone(); //I need to use it later Complex[,] d =...
View ArticleCompiling in a mingw64/msys2 environment
I'm trying to build my C program in a mingw64 (msys2) environment using cmake.I successfully compiled my example code with inkl mkl_malloc() call, but once including vdMul() function I get the...
View Article