FEAST 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 Articlecluster_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 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 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 Article