question about the complexity to calculate just one eigenvalue
As we know that the complexity of diagonalization of a matrix should be N^3, where N is the dimension of the matrix. In MKL, there are a lot of functions that can calculate just some of the...
View ArticleClik here to view.
Segmentation fault in mkl_pds_lp64_assemble_csr_full
Hello,I'm trying to use the following file with 3 MPI processes but I end up with the following trace:(gdb) bt #0 0x00007ffff52ac38f in mkl_pds_lp64_assemble_csr_full () from...
View ArticleVSL for Convolution
I'm getting started with using MKL VSL for convolutions. For learning, I first I want to conduct a 1D convolution. I'm confused with the usage example given here:...
View ArticleWrong results from cblas_sgemm
Here is the code:float max_val(const float * vec, size_t sz) { float result = 0.0f; for (size_t i = 0; i < sz; ++i) { float val = abs(vec[i]); if (val > result) result = val; } return result; }...
View ArticleIntel MKL dll redistribution
I have a question about the license and usage possibilities of the intel mkl .dll's. Is it possible to include the (unaltered) mkl .dll's* in a (Windows) product for our customers? I understand that...
View Articlequestion about the manual of scalapack
The parameter "lwork" of the function pdsyevr can be determined by some formula.I don't understand some of the formula:why the subscript is decreased by 1,e.g., mb_ -1?what is csrc_nb_? I can not find...
View Articlepardiso error with a particular set of data
Today in running a simulation test case I have found the following error:See attached txt file since by copying the test I received 'Your submission contains invalid characters ....".The compiler used...
View ArticleWarning signs when compiling Fortran 2008 code that uses...
Hello, I'm compiling a code that is linked to the MKL cluster sparse solver. This code is written in Fortran 2008 and compilation uses the -stand:f08 intel Fortran compiler flag. I use the include...
View Articlequestion about the blacs library in MKL
I can find a lot of examples use the syntax like Cblacs_get(...), Cblacs_pinfo(...), but I can not find the c-type declaration in the "mkl_cblacs.h". Is it legal to use this kind of c-wrapper function...
View Articleprogram sleep when set OMP_NUM_THREADS more than one
Dear all, I have written some codes with MKL, but recently I found this codes did not work if the set OMP_NUM_THREADS more than one. Here is an example code. I compile it with icc test.c -mkl. When I...
View ArticleChanges to the routine CGETRS between MKL 11.1.3 and MKL 11.3.3
Hi,After an upgrade from MKL 11.1.3 to MKL 11.3.3 I suspect that the mkl routine CGETRS gives slightly different results. This sounds very strange (since I believe it just forwards the calculation to...
View ArticleGuidance on integrating cluster_sparse_solver into my application
I am trying to integrate cluster_sparse_solver into my application, however, I am confused by this in the documentation:NoteMost of the input parameters (except for the pt, phase, and comm parameters...
View Articlein-core and out-of-core mode in pardiso
Hi, i am a beginner of MKL and am recently using the pardiso package as slover. I am hoping to run parallel cases by turn on the OpenMP in the solver. However, i do not know whether i should use the...
View ArticleIntel® Math Kernel Library (Intel® MKL) on CentOS 7, CPU IA-64 E5-2683v3
Hi all,I registered free license for MKL and installed on CentOS 7.2, CPU E5-2683v3 with errors:Error: A license for MKernL could not be obtained Your license restricts the platforms on which our...
View ArticleProblem compiling MKL with OpenMP
I'm trying to compile and load a code that uses MKL libraries to solve an eigenvalue problem. I use the following commands based on what the Intel MKL link line advisor tells me: ifort -c -i8 -qopenmp...
View ArticleCan't set number of threads >32 for Pardiso solver in MKL library
I'm running the pardiso solver on a machine with CentOS 6.5 using Fortran. I'm using MKL v 11.3When I set the number of threads using either environment variable MKL_NUM_THREADS or inside the Fortran...
View ArticleRe-use cluster pardiso factorization crashed
I'm trying to use cluster pardiso to solve large equations with complex hermitian matrix. I have multiple (>1000) right hand side but matrix A is the same, so I do phase 11 and 22 factorization...
View ArticleThe use of FFTW3 on my Xeon phi card (MIC)
hello , When i use the fftw3 in the MKL on the Intel Xeon Phi, i use the followed compile command : icpc -g -O0 -tbb -mkl=parallel...
View ArticleCNR setting does not work as expected
I am working on different machines (notebook and cloud, all Win7) and thus want to set CNR manually to get reproducible results (e.g. when working with Matlab).Thus, and using the latest MKL version, I...
View Articledf?dSearchCells1D give different result between Composer 16 and 16 Update3
Hello,I'm using the MKL in the Composer 16update 3 release and I notice a weird behaviour of the df?dSearchCells1D routine when a sitePoint is equal to the last breakPoint.I did a little test : const...
View Article