Pardiso scaling inversely with number of threads for openmp
Hello. After a number of failed debugging attempts and tests, I'm hoping to get some input on using pardiso in parallel with openmp. The software in question uses the intel fortran compiler (ifort)...
View ArticlePardiso scaling inversely with number of threads for openmp
Hello. After a number of failed debugging attempts and tests, I'm hoping to get some input on using pardiso in parallel with openmp. The software in question uses the intel fortran compiler (ifort)...
View ArticlePardiso scaling inversely (openmp)
Hello. After a number of failed debugging attempts and tests, I'm hoping to get some input on using pardiso in parallel with openmp. The software in question uses the intel fortran compiler (ifort)...
View ArticleProblem using pdhseqr
Hi,I have some trouble using pdhseqr from ScaLAPACK (it is supposed to compute the eigen value and Schur factorization of a matrix already in Hessenberg form). It happens for both pshseqr and pdhseqr....
View ArticleProblems in running MKL fortran example
Dear all,After the example fortran code vdrnggaussian.f successfully compiled (in visual studio 2010+Intel Parallel Studio XE 2013), the resulting executable file returns the following error...
View Articlepardiso multiple right hand side with sparse matrix
Hi All,Pardiso can handle multiple right side, AX=B, meaning B is a matrix. Every column of B is a right hand side. If B is very sparse, is there any way to exploit this feature, instead of storing B...
View ArticleLinking BLA95 and LAPACK 95 to a shared library in Linux
I'm in the process of porting an application from Windows to Linux, and one component of my application is a shared library (a DLL on Windows) compiled with Intel Fortran and with heavy use of the BLAS...
View ArticleUsing Intel® C++ Composer XE for Multiple Simple Random Sampling without...
<!--break-->IntroductionRandom sampling is often used when pre- or post-processing of all records of the entire data set is expensive, as in the following examples. When the file of records or...
View ArticleLarge Scale Weighted Least Squares
Hello,I would like to use Intel MKL to solve a Large Scale Weighted Least Squares problem.The Matrices are sparse yet their elements number might be huge. This is a result of working on images and...
View ArticleThe PZGESV subroutine
Hello everybody, Please I need your help to understand the scalapck subroutine PZGESV. I am trying to use this subroutine to solve a linear equation system A*X = B where A is a (N by N) matrix...
View ArticleMKL BLAS95 and LAPACK95 interface build
Hi, I'm trying to run the following MKL interface make files in RHEL MRG (3.10.33-rt32.34.el6rt.x86_64) in order to interface to MKL using the BLAS and LAPACK interface (using the mkl_ prefixed header...
View Articletranspose matrix stored in CSR format using MKL
All,I have some legacy parallel code that uses CSR format to store a very large, sparse matrix. In making some additions to the code, I have the need to transpose the matrix, storing the transpose in...
View ArticleOctave 3.8.1 with MKL
I am trying to rebuild GNU Octave with the Intel MKL on Ubuntu 14.04uname -aLinux zoli-Precision-WorkStation-T3500 3.13.0-27-generic #50-Ubuntu SMP Thu May 15 18:06:16 UTC 2014 x86_64 x86_64 x86_64...
View Articlepardiso_getdiag
When I do acall pardiso_getdiag( pt, df, da, mnum, error)the arrays DF and DA are filled with 'factored' and 'original' diagonal pivots.The pivots are in the same 'order' of the 'original' matrix or...
View ArticleIntel® MKL Cookbook Recipes
Intel MKL Users,We would like to Introduce a new feature Intel® MKL Cookbook, an online Document with recipes for assembling Intel MKL routines for solving complex problems.Please give us your valuable...
View ArticleError in the documentation of feast
Hi,I would like to report an error in the documentation of the reference manual for the extended eigensolver (feast). In the rci interface description, the parameter ijob is detailed. When ijob = 30 or...
View ArticleEigenvalue Solver Error (dfeast_scsrgv )
I'm trying to solve eigenvalue problem for large sparse matrices using dfeast_scsrgv function. The function works fine for small problems (ex: 8*8 sparse matrix) but it gives...
View Articlemkl_?csradd
Dear alllooks as if I need to use mkl_?csradd (C++).However, my csr matrix is zero based; mkl_?csradd does only support one based sparse matrices.I would appreciate any suggestions for a pragmatic...
View ArticleFGMRES
All,I just started learning to use the mkl and fgmres in general and I'm having troubles implementing/understanding the fgmres. I've looked at the example that came in the library and I've been able to...
View Articlelibiomp5md.dll location (release build)
Hello, where would I find the dll without tracing code? Currently I use one from MKL 11.1.2 (64-bit, version 5.0.2013.1126, file size 1043kB, modified 2014-01-31 12:23) and the profiler shows this...
View Article