MKL Pardiso full iterative solver not combined, and question about RCI ISS
Hi,We want to use Pardiso as full iterative solver, not combined direct - iterative, is there some way to do this.If not, please what's the best solver type to implement Parallel iterative solving, is...
View ArticleMKL DFT module not using AVX2 backend when the processor does support it
I have an application that uses MKL to perform DFTs. I'm using the latest release of MKL on CentOS 7. For a version reference, the MKL shared libraries are installed to...
View ArticleProblem in MKL example for CG with SSOR
I think there is something wrong with the sample "cg_ssor_precon_c.c" in the for loop that applies the SSOR iteration. The purpose of the loop is to apply SSOR repeatedly on the array &tmp[3*n]....
View ArticleVDEXP undefined on windows for mkl 2017 update 2
I have following compiling error using mkl 2017 update 2 on windows (It only happens on windows; it's ok on linux).error: identifier "VDEXP" is undefined VDEXP(&size, &x[0], &expx[0]);
View ArticlePardiso out of memory in phase 11
Hi,I have a problem with Pardiso running out of memory in the reordering/numerical factorization phase (11). It is not clear to me how the out of core option (iparm(60)) and environment variables...
View ArticleFGMRES preconditioner applied to?
I am using the MKL's preconditioned FGMRES solver and I am trying to understand what exactly is the vector that FGMRES is asking to apply the preconditioner to. From the reference for the solver,...
View ArticleClik here to view.
Poor scaling for real-to-real FFT with OpenMP
In the attached file I use MKL to compute a real-to-real FFT using OpenMP for multithreading.The code is compiled withicpc -o bench-fft -Wall -O3 -g -march=native -fopenmp bench-fft.cxx -mklThe machine...
View Articlenumroc returns incorrect value in scalapack
p.p1 {margin: 0.0px 0.0px 0.0px 0.0px; font: 11.0px Menlo; color: #ffe549; background-color: #2c67c8} span.s1 {font-variant-ligatures: no-common-ligatures}Thread Topic: Help Me
View ArticlePerforming split step fourier using MKL dfti
Hi,I am trying to perform split step fourier method with the following code. There seems to be a problem with the Fourier portion because when I compare my data with the values obtained with another...
View ArticleDFTI_REAL_REAL Speed?
I have some code that makes heavy use of 1-D DFTs using MKL (real-to-complex and complex-to-complex). I just realized that the library supports the DFTI_REAL_REAL layout, where the real and imaginary...
View ArticleTrust Region Random Results
We Implemented TR MKL Method for our purpose. It doesn't converge every time, and also produce different (random) solutions from the same initial condition and constraints.About Trust Region Size...
View ArticleEigenvalues not in ascending order
Hey there,is there a diagonalisation routine, that does not order the eigenvalues in an ascending way? Thanks in advance,sommerfeld
View ArticleStrange bad memory of MKL spline function
Dear all,I am using MKL compilers_and_libraries_2017.4.210 (update 3) and following the developer guide to construct a cubic spline interpolation workflow. However no matter how I configure my inputs...
View ArticleClik here to view.
MKL PARDISO: floating-point error with zero pivots ...
Hi,we have detected problems with the current version of the PARDISO solver.(MKL vers. 2017.4.210, Windows x64 architecture, static linking, Microsoft VS 2015).PARDISO produces a floating-point error...
View Articlelibmkl_blacs_openmpi compatibility
Hello everybody,I have a question regarding the libmkl_blacs_openmpi* libraries. Which Openmpi version is this library supposed to be compatible with ?I could not find this information in the usual MKL...
View Article"Trust Region Algorithm" Questions
Along this period, we have developed a calculation method that uses the Trust Region MKL API (with constraints).We found many difficulties, but after a lot of efforts we have obtained some quite good...
View ArticleAre LAPACKE_cgesdd and LAPACKE_cgesvd SVD calculations reliable?
I'm using LAPACKE_cgesdd and LAPACKE_cgesvd to compute the singular values of a matrix. Both the routines have the option to compute the singular values only. The problem I have is that, in the...
View ArticleSegmentation faults with sparse FEAST
Dear,I am using the interface dfeast_scsrev for computing eigenvalues and eigenvectors of a sparse matrix sorted using a CSR format (3-vector).It works fine with small sparse matrices with a size of...
View ArticleMKL library
SirI have to install a code. it requiyes linking of lapack n blas file. the code was written in 2009 using mkl 8 version. according to it for linking paths areLROOT = /opt/intel/mkl/lib/intel64/ LAPACK...
View ArticleDifference between cgesvd and LAPACKE_cgesvd or cgesdd and LAPACKE_cgesdd
I would like to know the difference between LAPACKE_cgesvd(see https://software.intel.com/en-us/node/521150) andcgesvd(see...
View Article