FGMRES 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 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 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 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 ArticleMKL FFT: fftw_mpi_plan_many_transpose
Hello,I have been trying to use the MKL FFT through the FFTW3 interface and now am stuck with plan creation. The following line of code produces a NULL value, which is subsequently caught and kills the...
View ArticleProblems about how to use multithreaded intel MKL
Hello， I have some problems about multithreaded intel MKL1. I use the function MKL_SET_NUM_THREADS(2) to change the thread number. Then I want to check the value of MKL_NUM_THREADS, by using...
View ArticleBasic linking problem with MKL library
Hi everybody,Long ago I used to use mkl library 9.1.027, and used it with great joy to run some programs that included a lot of operations on matrices. And it was easy to include in the program, just...
View ArticleAccess violation in dss_solve_real
Hello,We are calling MKL from C# code by using general mechanism described in article Using Intel® MKL in your C# program. The code is using DSS (Direct Sparse Solver) for a large system of equations...
View ArticleFortran95 Interface causes segfault
I am trying to use the zheevd or zheev with their Fortran95 interface:program main use lapack95 implicit none complex(8) :: H(2,2), w(2) H = 1d0 call zheevd(H,w, 'V') write (*,*) H write (*,*) "######"...
View Article