What is the memory management strategy in sparse_matrix_t ?
Hi all,I have two questions about memory management of sparse_matrix_t :1. When I create a matrix handle using routine mkl_sparse_?_create_csr, the handle will copy all data to new memory blocks? Or it...
View ArticleDoes MKL support real*16 datatype?
I am trying to convert some old Cray routines with double precision to the Windows Visual Fortran platform, but somehow am having difficulties in making the routines work. Just wondering if MKL have...
View Articledifference result of mkl_sparse_s_export_csr in two subroutines
Hi all,I am using mkl_sparse_s_export_csr function for reading data in sparse_matrix_t object, but the return result is error while in other function. Here is the example:#include "mkl_spblas.h"...
View Articledifferent results for LAPACKE_dgels between MKL and LAPACKE
Hi guys,I wrote some code to do dgels originally using LAPACKE. I also tried to use MKL to do the dgels. From my understanding, with the exact same code base but just recompile/relink with MKL library...
View ArticleFEAST Mathematica
Hi,I am running the FEAST algorithm on Mathematica/10.0. I have sparse matrices of dimension around 500,000, looking for around 50 eigenvalues. I have noticed the runtime depends massively on what...
View ArticleA strange result of sparse matrix addition with mol_sparse_s_add
Hi all，I have a new question about sparse matrix addition routine mkl_sparse_s_add, it return a strange result. And when I using double precision routine ml_sparse_d_add, everything is ok. Following is...
View ArticleWhat does the means of struct matrix_descr in mol_sparse_XXX routines?
Hi all, In many mkl_sparse_XXX routines (version 11.3), there is a a parameter of type matrix_descr, just like :sparse_status_t mkl_sparse_s_trsv (sparse_operation_t operation, float alpha, const...
View ArticleFEAST algorithm : feastinit input parameter setting problem
Hi,I have a problem with setting the input parameters at MKL FEAST algorithm library.(I currently use MKL 11.3 version and Intel parallel studio xe 2016 for c++ at Linux)I was trying to use the...
View ArticleBuild R with Intel MKL shared library
I build sucessfully R 3.2.2 with Intel MKL and ICC. Now I am wondering if library linking are indeed correct. Here is the output:% R CMD ldd BUILD/R-3.2.2/bin/exec/R linux-vdso.so.1...
View Articlepardiso error=-1, input inconsistent, what does it mean?
I am feeding pardiso a diagonal matrix with several zeros on the main diagonal. It is returning error=-1, "input inconsistent", when I was expecting -7 for "diagonal matrix is singular".What exactly...
View ArticleOpenMP MKL DGEMM Performance Issue
Hello,I am doing development on a 24-core machine (E5-2697-v2). When I launch a single DGEMM where the matrices are large (m=n=k=15,000), the performance improves as I increase the number of threads...
View ArticleUnexpected behavior of cluster_sparse_solver
Hello.I am trying to solve a system with `cluster_sparse_solver`, and am getting unexpected results. https://gist.github.com/ivan-krukov/157372d9a55db244c4b4In my case, given a sparse matrix `A`, we...
View Articledoes pardiso work for a 1 element matrix?
I am finding that pardiso does not work for a 1 element matrix. Calling pardiso with phase=11 returns error=0, which means success, but the pt(:) pointer array is all zero's. So when I subsequently...
View Articleneed help updating MKL software
The MKL version that I have has bugs, and I need to update it. I need help getting it updated.I am confused about what version of Intel Fortran I have installed. The Help About box in Visual Studio...
View ArticleWill Intel give me MKL 11.0.5 to replace 11.0.4 ?
I have encountered a confirmed bug with routine mkl_zcsrcoo in 11.0.4. It is my understanding that this bug has been fixed in 11.0.5. Would be it possible for me to get 11.0.5 ?
View ArticleData fitting of vector-valued function
Hello,I have a problem using the data fitting for a linear interpolation of a vector valued function. It seems that the format rhint is ignored by the function dfdInterpolate1D.#include...
View ArticleDSS Fortran Error
Hi all,Part of my Ph.D. research involves solving coupled reaction diffusion equations on arbitrary closed surfaces using Finite Element Method. In the past I had been using the ?GESV subroutine to...
View Articleneed help understanding error=-4 from pardiso
I am calling pardiso with mtype=11 (real nonsymmetric).I first call it with phase=11, iparm(1)=0, and the result is error=0. So all is good.I then call it with phase=22, and error=0. So all is still...
View ArticleCNR mode reporting incorrect SIMD version via KVM
HiWe had a problem with inconsistent results across some of our grid nodes, which I thought was worth sharing. After investigation we pinned this down to two different OS configurations returning...
View ArticleClik here to view.
Pardiso hangs in phase 33 when called from an OMP critical region in IVF 16.0
Hi all,I suspect this is a problem with my code, but I can't figure out what the issue is...I've been using similar code for a while now, and its worked fine until I installed 16.0, however it now...
View Article